CA1067576A - Gyromagnetic resonance spectroscopy employing spin echo spin-spin decoupling and two-dimensional spreading - Google Patents
Gyromagnetic resonance spectroscopy employing spin echo spin-spin decoupling and two-dimensional spreadingInfo
- Publication number
- CA1067576A CA1067576A CA278,232A CA278232A CA1067576A CA 1067576 A CA1067576 A CA 1067576A CA 278232 A CA278232 A CA 278232A CA 1067576 A CA1067576 A CA 1067576A
- Authority
- CA
- Canada
- Prior art keywords
- spin
- omega
- resonators
- resonance
- spin echo
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired
Links
- 238000004611 spectroscopical analysis Methods 0.000 title claims abstract description 9
- 230000007480 spreading Effects 0.000 title abstract description 3
- 238000000034 method Methods 0.000 claims abstract description 31
- 230000003595 spectral effect Effects 0.000 claims abstract description 28
- 238000001514 detection method Methods 0.000 claims abstract description 6
- 230000005415 magnetization Effects 0.000 claims description 17
- 239000013598 vector Substances 0.000 claims description 17
- 230000000694 effects Effects 0.000 claims description 14
- 230000003993 interaction Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 3
- 230000001427 coherent effect Effects 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 claims 4
- 206010013710 Drug interaction Diseases 0.000 claims 1
- 238000001228 spectrum Methods 0.000 abstract description 32
- 230000008878 coupling Effects 0.000 abstract description 12
- 238000010168 coupling process Methods 0.000 abstract description 12
- 238000005859 coupling reaction Methods 0.000 abstract description 12
- 238000002592 echocardiography Methods 0.000 abstract description 7
- 230000001808 coupling effect Effects 0.000 abstract description 6
- 238000005259 measurement Methods 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 10
- 230000009466 transformation Effects 0.000 description 6
- 238000004590 computer program Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000006698 induction Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000000655 nuclear magnetic resonance spectrum Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 102000004190 Enzymes Human genes 0.000 description 1
- 108090000790 Enzymes Proteins 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 102000004196 processed proteins & peptides Human genes 0.000 description 1
- 108090000765 processed proteins & peptides Proteins 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 238000009774 resonance method Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 229920001059 synthetic polymer Polymers 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/46—NMR spectroscopy
- G01R33/4633—Sequences for multi-dimensional NMR
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S505/00—Superconductor technology: apparatus, material, process
- Y10S505/825—Apparatus per se, device per se, or process of making or operating same
- Y10S505/842—Measuring and testing
- Y10S505/843—Electrical
- Y10S505/844—Nuclear magnetic resonance, NMR, system or device
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Lasers (AREA)
Abstract
Application for Patent of RICHARD R. ERNST
for GYROMAGNETIC RESONANCE SPECTROSCOPY EMPLOYING SPIN ECHO
SPIN SPIN DECOUPLING AND TWO-DIMENSIONAL SPREADING
ABSTRACT OF THE DISCLOSURE
A method of gyromagnetic resonance spectroscopy is disclosed which simplifies the resultant spectra by eliminating spin-spin coupling effects. The method is particularly useful for eliminating homonuclear spin-spin coupling effects. In the method, spin echo resonance of the resonators is obtained and the time elapsing between tipping of the resonators and detection of the resultant echo t1 is changed from one measure-ment to the next. The resultant echos are measured at equal intervals of time t2. Echo resonance data is Fourier trans-formed from the time domain into the frequency domain to obtain spectral data in-the frequency domain as a function of t1. A second Fourier transform of the resonance data is then effected to transform the spectral data into data which is a function of both .omega.1 and .omega.2 where .omega.1 and .omega.2 are each related to ? and ? respectively. The spectral data in the .omega.1 - .omega.2 plane is then projected at 45° to the .omega.2 and .omega.1 axes to derive spectral data free of spin-spin coupling effects. The data may also be projected perpendicular to and onto the .omega.1 axis to obtain a J coupling spectrum.
for GYROMAGNETIC RESONANCE SPECTROSCOPY EMPLOYING SPIN ECHO
SPIN SPIN DECOUPLING AND TWO-DIMENSIONAL SPREADING
ABSTRACT OF THE DISCLOSURE
A method of gyromagnetic resonance spectroscopy is disclosed which simplifies the resultant spectra by eliminating spin-spin coupling effects. The method is particularly useful for eliminating homonuclear spin-spin coupling effects. In the method, spin echo resonance of the resonators is obtained and the time elapsing between tipping of the resonators and detection of the resultant echo t1 is changed from one measure-ment to the next. The resultant echos are measured at equal intervals of time t2. Echo resonance data is Fourier trans-formed from the time domain into the frequency domain to obtain spectral data in-the frequency domain as a function of t1. A second Fourier transform of the resonance data is then effected to transform the spectral data into data which is a function of both .omega.1 and .omega.2 where .omega.1 and .omega.2 are each related to ? and ? respectively. The spectral data in the .omega.1 - .omega.2 plane is then projected at 45° to the .omega.2 and .omega.1 axes to derive spectral data free of spin-spin coupling effects. The data may also be projected perpendicular to and onto the .omega.1 axis to obtain a J coupling spectrum.
Description
1~6'7~7f~
B~CKG~OUND OF T~E INVENTION
The present invention relates in general to gyromagnetic resonance spectroscopy and more particularly to such spectroscopy employing spin echo and two-dimensional spreading techniques to obtain simplified spectra.
DESCRIPTION OF THE PRIOR ART
Heretofore, spin echo gyromagnetic resonance has been employed for obtaining a spectrum of spin-spin coupling constants which is free of magnetic field inhomogeneity effects and chemical shifts. Such a technique is disclosed and claimed in U.S.
Patent 3,753,0~1 issued August 14, 1973 and assigned to the same assignee as the present invention.
It is also known to resolve the multiplet spectral structure produced by coupled gyromagnetic resonators, such as hetero-nuclear coupling, by inducing a train of transient free induction decay resonances and detecting, during a time t2, the free induction decay resonance. A decoupling RF magnetic field is applied to excite one of the groups of resonators which is coupled to the other during the free induction decay resonance of the first group for decoupling the spins of the first and second groups~ The duration tl of the dccotlpling effect is changed from ono free induction decay resonance to the next.
The detected resonance data which is a function of two time lntervals, i.e., tl and t2, is then double Fourier transformed lnto the frequency domain and displayed as a two-dimensional plot for resolving the multiplet structure of the spectra of the first group of gyromagnetic bodies. Such a resonance technique is disclosed in U.S. Patent No. 4,045,723 issued August 30, 1977 and assigned to the same assignee as the present invention.
.
B~CKG~OUND OF T~E INVENTION
The present invention relates in general to gyromagnetic resonance spectroscopy and more particularly to such spectroscopy employing spin echo and two-dimensional spreading techniques to obtain simplified spectra.
DESCRIPTION OF THE PRIOR ART
Heretofore, spin echo gyromagnetic resonance has been employed for obtaining a spectrum of spin-spin coupling constants which is free of magnetic field inhomogeneity effects and chemical shifts. Such a technique is disclosed and claimed in U.S.
Patent 3,753,0~1 issued August 14, 1973 and assigned to the same assignee as the present invention.
It is also known to resolve the multiplet spectral structure produced by coupled gyromagnetic resonators, such as hetero-nuclear coupling, by inducing a train of transient free induction decay resonances and detecting, during a time t2, the free induction decay resonance. A decoupling RF magnetic field is applied to excite one of the groups of resonators which is coupled to the other during the free induction decay resonance of the first group for decoupling the spins of the first and second groups~ The duration tl of the dccotlpling effect is changed from ono free induction decay resonance to the next.
The detected resonance data which is a function of two time lntervals, i.e., tl and t2, is then double Fourier transformed lnto the frequency domain and displayed as a two-dimensional plot for resolving the multiplet structure of the spectra of the first group of gyromagnetic bodies. Such a resonance technique is disclosed in U.S. Patent No. 4,045,723 issued August 30, 1977 and assigned to the same assignee as the present invention.
.
-2-.' q~ '' ~
' ~
.~ . ... .. :
- . .. . .
1~67576 SUMMARY OF T~E PRESENT INVENTION
The principal object of the present invention is the provision of an improved method and apparatus for performing gyromagnetic resonance spectroscopy and particularly to an improved method and apparatus for simplifying spectra by eliminating undesired homonuclear spin-spin coupling effects.
According to the present invention,there is provided in a method of gyromagnetic resonance spectroscopy, operative upon an assembly of first and second groups of gyromagnetic resonators immersed in a polarizing magnetic field ~herein the magnetization vectors of said gyromagnetic resonators precess about said polarizing magnetic field, said first and second groups being coupled through spin-spin interaction, said spin-spin interaction causing a defocusing effect whereby . the phase coherents of precession among said resonators is progressively lost, the steps of:
:: ¦ periodically t.ipping the magnetization vectors of said spin-spin coupled first and second groups of gyromagnetic resonators with respect to said polarizing magnetic field;
periodically flipping the magnetization vectors of said ¦ precessing resonators relative to the direction of the polarizing magnetic field for reversing the defocusing effect I of the precession effect of the precession of said resonators bl about the polarizing magnetic field to obtain spin echo 1 resonance of æaid resonators;
, periodically detecting the spin echo resonance of said ;i resonators; and changing the period of time t between tipping of said magnetization vectors and detection of the spin echo resonance i 30 for a series of said periodically detected spin echo reson-ances to obtain detected gyromagnetic resonance data about said first and second spin-spin coupled groups of resonators
' ~
.~ . ... .. :
- . .. . .
1~67576 SUMMARY OF T~E PRESENT INVENTION
The principal object of the present invention is the provision of an improved method and apparatus for performing gyromagnetic resonance spectroscopy and particularly to an improved method and apparatus for simplifying spectra by eliminating undesired homonuclear spin-spin coupling effects.
According to the present invention,there is provided in a method of gyromagnetic resonance spectroscopy, operative upon an assembly of first and second groups of gyromagnetic resonators immersed in a polarizing magnetic field ~herein the magnetization vectors of said gyromagnetic resonators precess about said polarizing magnetic field, said first and second groups being coupled through spin-spin interaction, said spin-spin interaction causing a defocusing effect whereby . the phase coherents of precession among said resonators is progressively lost, the steps of:
:: ¦ periodically t.ipping the magnetization vectors of said spin-spin coupled first and second groups of gyromagnetic resonators with respect to said polarizing magnetic field;
periodically flipping the magnetization vectors of said ¦ precessing resonators relative to the direction of the polarizing magnetic field for reversing the defocusing effect I of the precession effect of the precession of said resonators bl about the polarizing magnetic field to obtain spin echo 1 resonance of æaid resonators;
, periodically detecting the spin echo resonance of said ;i resonators; and changing the period of time t between tipping of said magnetization vectors and detection of the spin echo resonance i 30 for a series of said periodically detected spin echo reson-ances to obtain detected gyromagnetic resonance data about said first and second spin-spin coupled groups of resonators
3 --from which to derive simplified gyromagnetic resonance spectral data.
In the described embodiment spin echo resonance is obtained of gyromagentic resonance samples having homonuclear spin-spin coupling. The spin echos are detected during a time inter~al t~ as a function of the change in the time tl between tipping of the resonators and detecting of the resonance.
This results in obtaining detected resonance data which is - the function of both tl and t2 and from which simplified spectral data may be derived.
I ~ In another feature of the described embodiment, the detected resonance data which is a function of both tl and t2 i5 double Fourier transformed into a ~1 and ~2 plane to obtain two-dimensional resonance data to facilitate resolving the multiplet structure of the spectra under analysis.
In another feature of the described embodiment, the double Fourier transformed resonance data in the frequency ¦ domain is projected at an angle to either the ~I or ~2 planes in such a manner as to obtain a simplified resonance spectrum of the sample under analysis.
In another feature of the described embodiment, the two-dimensional resonance data, in th~ frequency domain and ~ 2 plane, is projected at 45 to either the ~1 or ~2 ., axes to derive spectral data free of spin-spin coupling effects.
In another feature of the described embodiment, the two-dimensional resonance data in the frequency domain is projected in the ~ plane parallel to the ~2 axis onto the ~1 axis to deri~e a J spectrum of the sample under , .
.
' '~ .
! ~ - 3a -10~7576 ana].ysis.
Embodiments of the present invention will now be de.s.cribed, by way of example, with reference to the accompanying drawings in which:-Fig~ 1 is a schematic drawing, partly in block diagram form, of a gyromagnetic resonance spectrometer, Fig. 2 is a plot of RF nagnetic field intensity versus time depicting a method for exciting spin echo resonance of the sample under analysis, Fig. 3 is a plot similar to that of Fig. 2 depicting an 1 alternative method for exciting spin echo resonance of the sample under analysis, Fig. 4 is a schematic line diagram depicting the storage of resonance spectral data in the memory of the computer of ' Fig. 1, Fig. 5 is a diagram similar to that of Fig. 4 depicting a block of resonance data after Fourier transformation thereon in accordance l~ith the equation under Figs. 4 and 5, Fig. 6 is a plot depicting a block of resonance data corresponding to that o~ Fig. S a~tcr invcrting the matrix thereof in accordance with the equations under the block of Fig. 6, Fig. 7 is a schematic diagram depicting a block of spectral data stored in the memory of the computer of Fig. 1 and corresponding to a second Fourier transform thereof in accordance with thc formulas depicted below Fig. 7, ;' ' , ' .
~ 30 ' 1, ~
~ 4 -106~7~i I Fig. 8 is a diagram similar to that of Fig. 7 depicting ¦ inversion of the matrix of the data of Fig. 7 according to the ¦ equations depicted below Fig. 8, ¦ Fig. 9 is a schematic diagram of a block of resonance data S ¦ derived from the data of Fig. 8 in accordance with the equation ¦ depicted immediately below Fig. 9, ¦ Fig. 10 is a diagram similar to that of Fig. 9 depicting ¦ a block of resonance data stored in the memory of the computer - ¦ of Fig. 1 and derived from the data of Fig. 8 in accordance ¦ with the equation depicted below Fig. 10, ¦ Fig. 11 is a diagram similar to those of Figs. 9 and 10 ¦ depicting storage of the data of Figs. 9 and 10 and inter-polation thereof so as to equalize the frequency intervals in both the ~ and ~ directions, 15 I Fig. 12 is a schematic diagram similar to that of Fig. 11 ¦ depicting the step of projecting the spectral data of Fig. 11 I at an angle of 45 to the ~ axis, ¦ Fig. 13 is a two-dimensional display of the resonance ¦ spectral data and depicting, at waveform c, the projocted , 20 ¦ sp~ctrum derived ~rom per~orming th~ ~tcp o Fig. 12, ¦ Pig. 14 is a computer program flow diagram depicting operation of the computer of Fig. 1 for carrying out the resonance method of the described embodiment, Fig. lS is a computer program flow chart depicting the ' 251 echo train subroutlne portion of the flow chart of Fig. 14, and Fig. 16 lS a computer program flow chart depicting , operation of the computer program for projecting the spectral data, such program being executed at the termination of the program of Fig. 14.
:
` ~0~7576 DI~SCRIPTION OF THE PREFERRED_EMBODIMENTS
Referring now to Fig. 1 there is shown a gyromagnetie resonance spectrometer ll. sriefly, the speetrometer ll includes a container 12 for containing gyromagnetic resonators such as atomic nuelei or unpaired electrons to be analyzed.
In a typical example, the sample container 12 may contain relatively complex molecules such as biomolecules, enzymes, peptides, proteins or complicated organic moleeules in general.
A common transmitter/receiver coil 13 is disposed coaxially surrounding the container 12, such coil being wound in axial alignment with the Y axis of the Cartesian coordinate system indicated in Fig. 1. The single transmitter/receiver coil 13 ' is connected to a single coil gyromagnetic resonance spectro-i meter 14, such as a Varian Model CFT-20 or a Bruker Model , ,' SPX ~100.
`~ The sample under analysis is disposed in a relatively ,,~
intense unidirectional polarizing magnetic field Ho produced between the pole faces 15 and 16 of a relatively large electro-magnet, such as a 15" diameter pole face electromagnet or, in 1 20 a preferred embodiment, in the field o a supercollclucting I magnet having a magnetic ~ield intensity corresponding to a ¦ Larmor resonance frequency of the gyromagnetic resonators in the range of 220-360~-1z.
The spectrometer 14 is interfaced with a digltal computer 17, such as a Varian 620/L-100 having a 12K bit memory, via the intermediary of an analog-to-digital converter 18. One ; output of the computer is fed to a display print-out 19 for obtaining two-dimensional(2D)spectral displays of the resonance , i iO6'75'76 spectra of the sample under analysis. A typical 2D display is shown in Fig. 13. A synchronize and execute line 21 feeds signals from the computer 17 to the spectrometer 14 for placing the spectrometer under the control of computer 17.
S In operation, the spectrometer 14 is controlled by the ; computer 17 in such a manner as to excite spin echo resonance of the gyromagnetic resonators within the sample. In a typical example, the gyromagnetic resonators may comprise the protons of a relatively complex hydrocarbon molecule. The spectrometer includes an internal RF transmitter which applies a pulse 22 of radio frequency magnetic field to the sample under analysis . with the polarization vector of the RF magnetic field being at right angles to the direction of the polarizing magnetic field. The frequency of the radio frequency magnetic field lS is selected at the Larmor resonance frequency of the gyro-magnetic resonators under analysis and the intensity and the duration of the applied RF magnetic field are selected so as to tip the magnetization vectors of the gyromagnetic resonators at right angles to the direction of the polarizing magnetic field. This is indicated by pulse 22 of Fig. 2. When the resonators have been tippe~ by 9~ the RP pulse is terminated and the resonators begin to precess about the polarizing magnetic field Ho~ After a period of time corresponding to tl a second pulse of the radio frequency magnetic f;eld is applied with an intensity and a duration to flip the magnetizatio vectors of the precessing gyromagnetic resonators by 180 thereby reversing the defocusing effect of their precession about the direction of the polarizing magnetic field Ho~ This is indicated by pulse 23. In a time tl after the center of .~ , . .,",...
;' ".,, pulse 23 a spin echo resonance signal, resulting from application ¦ of pulses 22 and 23, will reach a maximum amplitude in the transmitter/receiver coil 13. The induced resonance signal is l picked up by the transmitter/receiver coil 13 and detected by 5 ¦ the spectrometer by sampling same at a number of equal intervals of time t2 commencing after tl. Thus, it is seen that at a ¦ time tl after application of the first pulse 22 the detection I of the spin echo resonance signal 24 commences and the spin ¦ echo resonance signal is detected at a number of equal intervals 10 ¦ of time t2. The detected resonance data is stored in the ¦ memory of the computer 17.
¦ As an alternative to application of a single 180 pulse 23 ¦ a succession of such pulses may be employed, as indicated in ¦ the method of Fig. 3, during the time interval tl. In the spin 15 ¦ echo method of Fig. 3, the first 180 pulse follows the 90 ¦ pulse by a period T and successive 180 pulses follow the ¦ first pulse by intervals of time 2T. The time tl between I tipping of the magnetization vectors and detection of the ¦ spin echo 24 remains at tl which is equal to 2nT.
20 I The method of Fig. 3 has the advalltage of providing better ; ¦ resolution because diffusion effects are compensated within ¦ the sample. However, it is slightly more complicated than the ¦ method of Fig. 2 and an additional complication in the method ¦ of Fig. 3 is that the peak amplitude of the detected spin echo 25 ¦ resonance changes sign from a negative at the first or every odd echo to positive for every even echo. Thus, when the I method of Fig. 3 is utilized it is necessary to change the sign ¦ of every second recorded echo or to change the phase of the RF
I transmitter pulse.
I
~ ~6'7576 In accordance with the present invention, the detected spin echo resonance is measured at a plurality of successive time displaced intervals such as 64 intervals of t2, for each value of tl. Each sampled value of the detected spin echo resonance signal, for a given spin echo signal, is recorded in a corresponding location in the memory of the computer 17, as depicted in the storage data block of ~ig. 4. There are then m number of spin echo resonance signals recorded, there being one for each different value of tl. This is indicated by each different row in Fig. 4. In a simple example consider only eight resonance data values Mjk for each spin echo resonance signal and assume there are only eight echo signals, each corresponding to a different value of tl. The resultant 64 resonance data signal values are then stored in a storage block as indicated in Fig. 4.
In the next step the resonance data, indicated by Fig. 4, - is Fourier transformed from the time domain into the frequency domain with the resultant data stored in a memory block in the manner as indicated in Fig. 5. In the process of the Fourier transform of the data of the block of Fig. 4 into the block of Fig. 5, eight zero values are added to each row of the data o Fig. 4 prior to effecting the Fourier transform into the frequency domain according to the equation depicted below Figs. 4 and 5. This permits the real and imaginary portions o the resonance spectral data to be obtained and stored in the memory block as indicated in Fig. 5.
` Next, the spectral datà of the memory block of Fig. 5 is inverted according to the equations at the bottom of Fig. 6 and stored a memory block in the manner as depicted in Fi~. 6.
g .
' l l '~......... . I ~: I
,: .
106~57~j ~ Next, an additional equal number of zero values are added ¦ to each rol~ of the memory block of Fig. 6. Then the data of ¦ the expanded memory block of Fig. 6 is again Fourier transformed and stored in a memory block as indicated in Fig. 7. The 5 ¦ second Fourier transformation of ~he data of the expanded memory ¦ block of Fig. 6 into the data of block 7 is done in accordance ¦ with the Fourier transform equations depicted below Fig. 7.
¦ Next, the spectral data in the frequency domain of Fig. 7 ¦ is inverted and storedin the inverted form as shown in the 10 ¦ memory block of Fig. 8. The inversions are in accordance ¦ with the equations depicted below Fig. 8.
¦ Next, the real and imaginary spectral data of the memory ¦ block of Fig. 8 is converted into absolute value spectral data ¦ for positive frequency values of ~ and stored in a memory 15 ¦ block as depicted in Fig. 9. The transformation of the positive ¦ frequency data from the memory block of Fig. 8 into the data block of Fig. 9 is in accordance with the equation below Fig. 9.
In this step of the process, the memory data block is reduced to a matrix 8x8. Similarly, the negative frequency ~ data is converted from the data of the memory block of Pig. 8 into the data of the memory block of Fig. 10 in accordance with the equation below Fig. 10.
Next, the data from the memory blocks of Figs. 9 and 10 -is combined into one memory data block with equal increments of frequency along both the ~ and ~ axes as shown in Fig. 11.
In Fig. 11 it is assumed that the ~ increment ~ is four times the magnitude of the ~ increment A~ and therefore interpolated values Sl/4, Sl/2, S3/4 are inserted into the rows extending In he ~ direction so that equal increments of ~ -10-., I
~ ~ 106~576 frequency are obtained along both the ~ and ~ axes. Inter-polation is merely a linear interpolation between the two adja-cent previously recorded signal data values. In addition, the negative ~ frequency data is interpolated and stored in the memory block of Fig. 11. Equalizing the frequency increments along the orthogonal ~ and ~ axes facilitates a 45 projection of the resonance data of the block of Fig. 11.
The projection, which is the next step, is shown in Fig. 12 and is readily accomplished by shifting the data of successive rows in the ~ direction to the right by one increment. The projection is obtained by summing in the ~ direction the various signal values for a given column. The projected sums are thell displayed along the ~ axis to obtain a greatly simp]ified resonance spectrum as shown by waveform c of Fig. 13.
Waveform c is a spectrum free of spin-spin multiplet structure.
Waveform a of Fig. 13 shows the resonance data spectrum as obtained with all of the spin-spin coupling effects present and corresponding to a projection of the spectral data of the memory block of Fig. 11 onto the ~ axis. The coupled spectrum ~a) shows six chemically shifted groups with each group split into ovcrlapping multiplot lines due to homonuclear spin-spin coupling. The plot o Fig. 13 is an isometric projection and display of the resonance data in the block of ~ig. 11.
In a next step, the resonance data of the block of Fig. 11 is projected onto the ~ axis by summing the values of each row onto the ~ axis to obtain a J coupling spectrum of the sample under analysis. The J coupling spectrum yields important and valuable data as previously disclosed in the aforecited U.S, Patent 3,753,081.
.,`, . , .
,~, . .
~.` ~ , . . I
. .
In the method as above described, it may be advantageous to suppress weak background nuclear magnetic resonant signals.
This can be done before the aforedescribed projections by using a rather coarse digitization of the data which suppresses weak signals. In addition, time averaging obtained during the projection smooths the digitization steps almost completely.
In addition, it is often not desirable to record complete 2D spectra of complicated molecules. This may require too much memory. Therefore, in such cases a 2D spread of only a selected spectral region is employed. This is conveniently accomplished by recording complete echos of the entire spectrum and trans-. forming the same in the first Fourier transform step of Fig. 5.After the first Fourier transorm step, the interesting spectral range is selected, stored and used in the second Fourier transformation of Fig. 7. As another alternative, an analog filter is employed to select the resonances which are of particular interest in the response. Also, the number of echos which have to be recorded to obtain sufficient resolution in the ~ direction is rather limited due to the small range of coupling constants in the ~ direction.
The advantage of the present invclltion is that it permits great simplificatiGn ln the recorded spectra for purposes of analysis without losing information. The technique permits a two-dimensional spread of complicated nuclear magnetic resonance spectra, i.e., of biomolecules or synthetic polymers. At the same time, the techn.ique permits complete homonuclear decoupling of the NMR spectra. This is particularly useful for proton spectra and has a special utility in connection with blochemical application hich are ~enerally limited to proton spectroscopy -i2~
.' I .
and which spectra are particu~arly complicated by the numerous ¦ spin-spin couplings. The method of the present invention permits, ¦ for the first time, the operator to untangle extremely compli-¦ cated spectra by means of the two dimensional spread and 5 ¦ complete decoupling. The method of the presént invention is ¦ particularly useful with regard to relatively weakly coupled ¦ spin-spin systems. Most biological applications use high-field ¦ spectrometers, i.e., superconducting high-field spectrometers ¦ in the Larmor resonance range of 220-360 ~egahertz. At these 10 ¦ high-fields, most spectra are weakly coupled.
¦ The theory behind the present invention is that echo ¦ amplitudes in spin echo experiments are not effected by the ¦ chemical shift and reflect exclusively the effects of nuclear ¦ spin-spin coupling constants and of relaxation phenomena as 15 I long as the couplings are sufficiently weak. The free decay ¦ of the individual echos, on the other hand, is governed by ¦ the complete nuclear Hamaltonian. The free decay of an echo ¦ is composed of the various magnetization vectors Mjk~tl,t2).
They describe thc observable transverse magnetization of the ¦ resonance line k in the multiplet belonging to a set j oE
magnetically equivalent nuclei with thc Zeeman requency ~j, I Mjk~tl,t2) ~ Mjk~O,O) cos (~L jktl ~ jkt2) ~Xp(-tl/T2jk-t2/T2jk) I with ~jk ~j + ~ jk. The multiplet splitting is denoted by i~ l ~ jk = 2~ Jjlmlk with the coupling constants Jil and the 25 1 magnetic quantum nuMbers mlk of nucleus 1 . T2jk is the trans-verse relaxation time of resonance line jk and T2jk includes ¦ additionally the effects of magnetic field inhomogeneity. The two time parameters tl and t2 are defined in Fig. 1.
'.~ . ' .
~ 1~- .~
. . ~ . . . ~
.
,-- ~ - .. ..
I 1~67576 l To obtain a ~D J resolved spectrum~ a complete set o~
¦ echos for various tl values is recorded. A two-dimensional I Fourier transformation of ~Ittl,t2) produces the 2D spectrum ¦ S(~ ,~ ). The contribution of the magnetization component - 5 1 Mjk(tl,t2) to the absolute value ¦S¦~ "~ ) is given for ~0 by lS I jk(~ '~ )=2Mjk(0,0)[1/T2jk+(~ k) ] [1/T2~k ¦ 2 jk A 2D J-resolved spectrum of a composite sample is shown I in Fig. 13 at (c). Each pea~ of the original spectrum (Fig. 13 10 ¦ at ~a) is represented by a peak with the proper intensity in ¦ the 2D plot. The selectivity along ~ is given by the unper-¦ turbed resonance ~requencies ~jk~ whereas the spread in the ~
¦ direction is determined exclusively by the multiplet splitting ~~ jk~ A comparison with the original spectrum shows that the 15 ¦ multiplet resolution has been significantly enhanced, making I the analysis even of very complicated patterns rat'ner easy.
¦ It is crucial to notice that the peaks k of each multiplet ¦ j lie on a straight line passing through the point ~j on the I; ¦ ~ axis. By mea]ls of a projection of the 2D spec~rum along '!~ 20 1 tllis direction onto the ~ axis it is now possible to obtain ¦ a completely decoupled spectrum. This is demonstrated by Fig.
13 at c which clearly shows the six peaks corresponding to the six sets of nonequivalent protons in the sample. The obtained I resolution is severely limited by the 64 samples used to 25 ¦ represent the spectrum. There is no principle limitation to obt~in much better resolution by using a larger data array.
1.
I
~ ,' .
.
:: .~, .
_ 1 1067S'76 It is by no means necessary to record a complete 2D
spectrum to obtain information about a particular frequency region. It is easily possible to select a narrow portion of ¦ the spectrum after the first Fourier transformation of the 5 1 various echos with respect to t2.
A computer flow chart for the computer program is shown in Figs. 14-16 and the suitable computer sourcc program in Varian assembly language for use with a DAS assembler and a 620 series ¦ Varian data machine computer, commercially available from 10 ¦ Varian Data Machines of Irvine, California.
I r~
In the described embodiment spin echo resonance is obtained of gyromagentic resonance samples having homonuclear spin-spin coupling. The spin echos are detected during a time inter~al t~ as a function of the change in the time tl between tipping of the resonators and detecting of the resonance.
This results in obtaining detected resonance data which is - the function of both tl and t2 and from which simplified spectral data may be derived.
I ~ In another feature of the described embodiment, the detected resonance data which is a function of both tl and t2 i5 double Fourier transformed into a ~1 and ~2 plane to obtain two-dimensional resonance data to facilitate resolving the multiplet structure of the spectra under analysis.
In another feature of the described embodiment, the double Fourier transformed resonance data in the frequency ¦ domain is projected at an angle to either the ~I or ~2 planes in such a manner as to obtain a simplified resonance spectrum of the sample under analysis.
In another feature of the described embodiment, the two-dimensional resonance data, in th~ frequency domain and ~ 2 plane, is projected at 45 to either the ~1 or ~2 ., axes to derive spectral data free of spin-spin coupling effects.
In another feature of the described embodiment, the two-dimensional resonance data in the frequency domain is projected in the ~ plane parallel to the ~2 axis onto the ~1 axis to deri~e a J spectrum of the sample under , .
.
' '~ .
! ~ - 3a -10~7576 ana].ysis.
Embodiments of the present invention will now be de.s.cribed, by way of example, with reference to the accompanying drawings in which:-Fig~ 1 is a schematic drawing, partly in block diagram form, of a gyromagnetic resonance spectrometer, Fig. 2 is a plot of RF nagnetic field intensity versus time depicting a method for exciting spin echo resonance of the sample under analysis, Fig. 3 is a plot similar to that of Fig. 2 depicting an 1 alternative method for exciting spin echo resonance of the sample under analysis, Fig. 4 is a schematic line diagram depicting the storage of resonance spectral data in the memory of the computer of ' Fig. 1, Fig. 5 is a diagram similar to that of Fig. 4 depicting a block of resonance data after Fourier transformation thereon in accordance l~ith the equation under Figs. 4 and 5, Fig. 6 is a plot depicting a block of resonance data corresponding to that o~ Fig. S a~tcr invcrting the matrix thereof in accordance with the equations under the block of Fig. 6, Fig. 7 is a schematic diagram depicting a block of spectral data stored in the memory of the computer of Fig. 1 and corresponding to a second Fourier transform thereof in accordance with thc formulas depicted below Fig. 7, ;' ' , ' .
~ 30 ' 1, ~
~ 4 -106~7~i I Fig. 8 is a diagram similar to that of Fig. 7 depicting ¦ inversion of the matrix of the data of Fig. 7 according to the ¦ equations depicted below Fig. 8, ¦ Fig. 9 is a schematic diagram of a block of resonance data S ¦ derived from the data of Fig. 8 in accordance with the equation ¦ depicted immediately below Fig. 9, ¦ Fig. 10 is a diagram similar to that of Fig. 9 depicting ¦ a block of resonance data stored in the memory of the computer - ¦ of Fig. 1 and derived from the data of Fig. 8 in accordance ¦ with the equation depicted below Fig. 10, ¦ Fig. 11 is a diagram similar to those of Figs. 9 and 10 ¦ depicting storage of the data of Figs. 9 and 10 and inter-polation thereof so as to equalize the frequency intervals in both the ~ and ~ directions, 15 I Fig. 12 is a schematic diagram similar to that of Fig. 11 ¦ depicting the step of projecting the spectral data of Fig. 11 I at an angle of 45 to the ~ axis, ¦ Fig. 13 is a two-dimensional display of the resonance ¦ spectral data and depicting, at waveform c, the projocted , 20 ¦ sp~ctrum derived ~rom per~orming th~ ~tcp o Fig. 12, ¦ Pig. 14 is a computer program flow diagram depicting operation of the computer of Fig. 1 for carrying out the resonance method of the described embodiment, Fig. lS is a computer program flow chart depicting the ' 251 echo train subroutlne portion of the flow chart of Fig. 14, and Fig. 16 lS a computer program flow chart depicting , operation of the computer program for projecting the spectral data, such program being executed at the termination of the program of Fig. 14.
:
` ~0~7576 DI~SCRIPTION OF THE PREFERRED_EMBODIMENTS
Referring now to Fig. 1 there is shown a gyromagnetie resonance spectrometer ll. sriefly, the speetrometer ll includes a container 12 for containing gyromagnetic resonators such as atomic nuelei or unpaired electrons to be analyzed.
In a typical example, the sample container 12 may contain relatively complex molecules such as biomolecules, enzymes, peptides, proteins or complicated organic moleeules in general.
A common transmitter/receiver coil 13 is disposed coaxially surrounding the container 12, such coil being wound in axial alignment with the Y axis of the Cartesian coordinate system indicated in Fig. 1. The single transmitter/receiver coil 13 ' is connected to a single coil gyromagnetic resonance spectro-i meter 14, such as a Varian Model CFT-20 or a Bruker Model , ,' SPX ~100.
`~ The sample under analysis is disposed in a relatively ,,~
intense unidirectional polarizing magnetic field Ho produced between the pole faces 15 and 16 of a relatively large electro-magnet, such as a 15" diameter pole face electromagnet or, in 1 20 a preferred embodiment, in the field o a supercollclucting I magnet having a magnetic ~ield intensity corresponding to a ¦ Larmor resonance frequency of the gyromagnetic resonators in the range of 220-360~-1z.
The spectrometer 14 is interfaced with a digltal computer 17, such as a Varian 620/L-100 having a 12K bit memory, via the intermediary of an analog-to-digital converter 18. One ; output of the computer is fed to a display print-out 19 for obtaining two-dimensional(2D)spectral displays of the resonance , i iO6'75'76 spectra of the sample under analysis. A typical 2D display is shown in Fig. 13. A synchronize and execute line 21 feeds signals from the computer 17 to the spectrometer 14 for placing the spectrometer under the control of computer 17.
S In operation, the spectrometer 14 is controlled by the ; computer 17 in such a manner as to excite spin echo resonance of the gyromagnetic resonators within the sample. In a typical example, the gyromagnetic resonators may comprise the protons of a relatively complex hydrocarbon molecule. The spectrometer includes an internal RF transmitter which applies a pulse 22 of radio frequency magnetic field to the sample under analysis . with the polarization vector of the RF magnetic field being at right angles to the direction of the polarizing magnetic field. The frequency of the radio frequency magnetic field lS is selected at the Larmor resonance frequency of the gyro-magnetic resonators under analysis and the intensity and the duration of the applied RF magnetic field are selected so as to tip the magnetization vectors of the gyromagnetic resonators at right angles to the direction of the polarizing magnetic field. This is indicated by pulse 22 of Fig. 2. When the resonators have been tippe~ by 9~ the RP pulse is terminated and the resonators begin to precess about the polarizing magnetic field Ho~ After a period of time corresponding to tl a second pulse of the radio frequency magnetic f;eld is applied with an intensity and a duration to flip the magnetizatio vectors of the precessing gyromagnetic resonators by 180 thereby reversing the defocusing effect of their precession about the direction of the polarizing magnetic field Ho~ This is indicated by pulse 23. In a time tl after the center of .~ , . .,",...
;' ".,, pulse 23 a spin echo resonance signal, resulting from application ¦ of pulses 22 and 23, will reach a maximum amplitude in the transmitter/receiver coil 13. The induced resonance signal is l picked up by the transmitter/receiver coil 13 and detected by 5 ¦ the spectrometer by sampling same at a number of equal intervals of time t2 commencing after tl. Thus, it is seen that at a ¦ time tl after application of the first pulse 22 the detection I of the spin echo resonance signal 24 commences and the spin ¦ echo resonance signal is detected at a number of equal intervals 10 ¦ of time t2. The detected resonance data is stored in the ¦ memory of the computer 17.
¦ As an alternative to application of a single 180 pulse 23 ¦ a succession of such pulses may be employed, as indicated in ¦ the method of Fig. 3, during the time interval tl. In the spin 15 ¦ echo method of Fig. 3, the first 180 pulse follows the 90 ¦ pulse by a period T and successive 180 pulses follow the ¦ first pulse by intervals of time 2T. The time tl between I tipping of the magnetization vectors and detection of the ¦ spin echo 24 remains at tl which is equal to 2nT.
20 I The method of Fig. 3 has the advalltage of providing better ; ¦ resolution because diffusion effects are compensated within ¦ the sample. However, it is slightly more complicated than the ¦ method of Fig. 2 and an additional complication in the method ¦ of Fig. 3 is that the peak amplitude of the detected spin echo 25 ¦ resonance changes sign from a negative at the first or every odd echo to positive for every even echo. Thus, when the I method of Fig. 3 is utilized it is necessary to change the sign ¦ of every second recorded echo or to change the phase of the RF
I transmitter pulse.
I
~ ~6'7576 In accordance with the present invention, the detected spin echo resonance is measured at a plurality of successive time displaced intervals such as 64 intervals of t2, for each value of tl. Each sampled value of the detected spin echo resonance signal, for a given spin echo signal, is recorded in a corresponding location in the memory of the computer 17, as depicted in the storage data block of ~ig. 4. There are then m number of spin echo resonance signals recorded, there being one for each different value of tl. This is indicated by each different row in Fig. 4. In a simple example consider only eight resonance data values Mjk for each spin echo resonance signal and assume there are only eight echo signals, each corresponding to a different value of tl. The resultant 64 resonance data signal values are then stored in a storage block as indicated in Fig. 4.
In the next step the resonance data, indicated by Fig. 4, - is Fourier transformed from the time domain into the frequency domain with the resultant data stored in a memory block in the manner as indicated in Fig. 5. In the process of the Fourier transform of the data of the block of Fig. 4 into the block of Fig. 5, eight zero values are added to each row of the data o Fig. 4 prior to effecting the Fourier transform into the frequency domain according to the equation depicted below Figs. 4 and 5. This permits the real and imaginary portions o the resonance spectral data to be obtained and stored in the memory block as indicated in Fig. 5.
` Next, the spectral datà of the memory block of Fig. 5 is inverted according to the equations at the bottom of Fig. 6 and stored a memory block in the manner as depicted in Fi~. 6.
g .
' l l '~......... . I ~: I
,: .
106~57~j ~ Next, an additional equal number of zero values are added ¦ to each rol~ of the memory block of Fig. 6. Then the data of ¦ the expanded memory block of Fig. 6 is again Fourier transformed and stored in a memory block as indicated in Fig. 7. The 5 ¦ second Fourier transformation of ~he data of the expanded memory ¦ block of Fig. 6 into the data of block 7 is done in accordance ¦ with the Fourier transform equations depicted below Fig. 7.
¦ Next, the spectral data in the frequency domain of Fig. 7 ¦ is inverted and storedin the inverted form as shown in the 10 ¦ memory block of Fig. 8. The inversions are in accordance ¦ with the equations depicted below Fig. 8.
¦ Next, the real and imaginary spectral data of the memory ¦ block of Fig. 8 is converted into absolute value spectral data ¦ for positive frequency values of ~ and stored in a memory 15 ¦ block as depicted in Fig. 9. The transformation of the positive ¦ frequency data from the memory block of Fig. 8 into the data block of Fig. 9 is in accordance with the equation below Fig. 9.
In this step of the process, the memory data block is reduced to a matrix 8x8. Similarly, the negative frequency ~ data is converted from the data of the memory block of Pig. 8 into the data of the memory block of Fig. 10 in accordance with the equation below Fig. 10.
Next, the data from the memory blocks of Figs. 9 and 10 -is combined into one memory data block with equal increments of frequency along both the ~ and ~ axes as shown in Fig. 11.
In Fig. 11 it is assumed that the ~ increment ~ is four times the magnitude of the ~ increment A~ and therefore interpolated values Sl/4, Sl/2, S3/4 are inserted into the rows extending In he ~ direction so that equal increments of ~ -10-., I
~ ~ 106~576 frequency are obtained along both the ~ and ~ axes. Inter-polation is merely a linear interpolation between the two adja-cent previously recorded signal data values. In addition, the negative ~ frequency data is interpolated and stored in the memory block of Fig. 11. Equalizing the frequency increments along the orthogonal ~ and ~ axes facilitates a 45 projection of the resonance data of the block of Fig. 11.
The projection, which is the next step, is shown in Fig. 12 and is readily accomplished by shifting the data of successive rows in the ~ direction to the right by one increment. The projection is obtained by summing in the ~ direction the various signal values for a given column. The projected sums are thell displayed along the ~ axis to obtain a greatly simp]ified resonance spectrum as shown by waveform c of Fig. 13.
Waveform c is a spectrum free of spin-spin multiplet structure.
Waveform a of Fig. 13 shows the resonance data spectrum as obtained with all of the spin-spin coupling effects present and corresponding to a projection of the spectral data of the memory block of Fig. 11 onto the ~ axis. The coupled spectrum ~a) shows six chemically shifted groups with each group split into ovcrlapping multiplot lines due to homonuclear spin-spin coupling. The plot o Fig. 13 is an isometric projection and display of the resonance data in the block of ~ig. 11.
In a next step, the resonance data of the block of Fig. 11 is projected onto the ~ axis by summing the values of each row onto the ~ axis to obtain a J coupling spectrum of the sample under analysis. The J coupling spectrum yields important and valuable data as previously disclosed in the aforecited U.S, Patent 3,753,081.
.,`, . , .
,~, . .
~.` ~ , . . I
. .
In the method as above described, it may be advantageous to suppress weak background nuclear magnetic resonant signals.
This can be done before the aforedescribed projections by using a rather coarse digitization of the data which suppresses weak signals. In addition, time averaging obtained during the projection smooths the digitization steps almost completely.
In addition, it is often not desirable to record complete 2D spectra of complicated molecules. This may require too much memory. Therefore, in such cases a 2D spread of only a selected spectral region is employed. This is conveniently accomplished by recording complete echos of the entire spectrum and trans-. forming the same in the first Fourier transform step of Fig. 5.After the first Fourier transorm step, the interesting spectral range is selected, stored and used in the second Fourier transformation of Fig. 7. As another alternative, an analog filter is employed to select the resonances which are of particular interest in the response. Also, the number of echos which have to be recorded to obtain sufficient resolution in the ~ direction is rather limited due to the small range of coupling constants in the ~ direction.
The advantage of the present invclltion is that it permits great simplificatiGn ln the recorded spectra for purposes of analysis without losing information. The technique permits a two-dimensional spread of complicated nuclear magnetic resonance spectra, i.e., of biomolecules or synthetic polymers. At the same time, the techn.ique permits complete homonuclear decoupling of the NMR spectra. This is particularly useful for proton spectra and has a special utility in connection with blochemical application hich are ~enerally limited to proton spectroscopy -i2~
.' I .
and which spectra are particu~arly complicated by the numerous ¦ spin-spin couplings. The method of the present invention permits, ¦ for the first time, the operator to untangle extremely compli-¦ cated spectra by means of the two dimensional spread and 5 ¦ complete decoupling. The method of the presént invention is ¦ particularly useful with regard to relatively weakly coupled ¦ spin-spin systems. Most biological applications use high-field ¦ spectrometers, i.e., superconducting high-field spectrometers ¦ in the Larmor resonance range of 220-360 ~egahertz. At these 10 ¦ high-fields, most spectra are weakly coupled.
¦ The theory behind the present invention is that echo ¦ amplitudes in spin echo experiments are not effected by the ¦ chemical shift and reflect exclusively the effects of nuclear ¦ spin-spin coupling constants and of relaxation phenomena as 15 I long as the couplings are sufficiently weak. The free decay ¦ of the individual echos, on the other hand, is governed by ¦ the complete nuclear Hamaltonian. The free decay of an echo ¦ is composed of the various magnetization vectors Mjk~tl,t2).
They describe thc observable transverse magnetization of the ¦ resonance line k in the multiplet belonging to a set j oE
magnetically equivalent nuclei with thc Zeeman requency ~j, I Mjk~tl,t2) ~ Mjk~O,O) cos (~L jktl ~ jkt2) ~Xp(-tl/T2jk-t2/T2jk) I with ~jk ~j + ~ jk. The multiplet splitting is denoted by i~ l ~ jk = 2~ Jjlmlk with the coupling constants Jil and the 25 1 magnetic quantum nuMbers mlk of nucleus 1 . T2jk is the trans-verse relaxation time of resonance line jk and T2jk includes ¦ additionally the effects of magnetic field inhomogeneity. The two time parameters tl and t2 are defined in Fig. 1.
'.~ . ' .
~ 1~- .~
. . ~ . . . ~
.
,-- ~ - .. ..
I 1~67576 l To obtain a ~D J resolved spectrum~ a complete set o~
¦ echos for various tl values is recorded. A two-dimensional I Fourier transformation of ~Ittl,t2) produces the 2D spectrum ¦ S(~ ,~ ). The contribution of the magnetization component - 5 1 Mjk(tl,t2) to the absolute value ¦S¦~ "~ ) is given for ~0 by lS I jk(~ '~ )=2Mjk(0,0)[1/T2jk+(~ k) ] [1/T2~k ¦ 2 jk A 2D J-resolved spectrum of a composite sample is shown I in Fig. 13 at (c). Each pea~ of the original spectrum (Fig. 13 10 ¦ at ~a) is represented by a peak with the proper intensity in ¦ the 2D plot. The selectivity along ~ is given by the unper-¦ turbed resonance ~requencies ~jk~ whereas the spread in the ~
¦ direction is determined exclusively by the multiplet splitting ~~ jk~ A comparison with the original spectrum shows that the 15 ¦ multiplet resolution has been significantly enhanced, making I the analysis even of very complicated patterns rat'ner easy.
¦ It is crucial to notice that the peaks k of each multiplet ¦ j lie on a straight line passing through the point ~j on the I; ¦ ~ axis. By mea]ls of a projection of the 2D spec~rum along '!~ 20 1 tllis direction onto the ~ axis it is now possible to obtain ¦ a completely decoupled spectrum. This is demonstrated by Fig.
13 at c which clearly shows the six peaks corresponding to the six sets of nonequivalent protons in the sample. The obtained I resolution is severely limited by the 64 samples used to 25 ¦ represent the spectrum. There is no principle limitation to obt~in much better resolution by using a larger data array.
1.
I
~ ,' .
.
:: .~, .
_ 1 1067S'76 It is by no means necessary to record a complete 2D
spectrum to obtain information about a particular frequency region. It is easily possible to select a narrow portion of ¦ the spectrum after the first Fourier transformation of the 5 1 various echos with respect to t2.
A computer flow chart for the computer program is shown in Figs. 14-16 and the suitable computer sourcc program in Varian assembly language for use with a DAS assembler and a 620 series ¦ Varian data machine computer, commercially available from 10 ¦ Varian Data Machines of Irvine, California.
I r~
Claims (14)
1. In a method of gyromagnetic resonance spectroscopy, operative upon an assembly of first and second groups of gyromagnetic resonators immersed in a polarizing magnetic field wherein the magnetization vectors of said gyromagnetic resonators process about said polarizing magnetic field, said first and second groups being coupled through spin-spin interaction, said spin-spin interaction causing a defocusing effect whereby the phase coherents of precession among said resonators is progressively lost, the steps of:
periodically tipping the magnetization vectors of said spin-spin coupled first and second groups of gyromagnetic resonators with respect to said polarizing magnetic field;
periodically flipping the magnetization vectors of said precessing resonators relative to the direction of the polarizing magnetic field for reversing the defocusing effect of the precession effect of the precession of said resonators about the polarizing magnetic field to obtain spin echo resonance of said resonators;
periodically detecting the spin echo resonance of said resonators; and changing the period of time t between tipping of said magnetization vectors and detection of the spin echo resonance for a series of said periodically detected spin echo resonances to obtain detected gyromagnetic resonance data about said first and second spin-spin coupled groups of resonators from which to derive simplified gyromagnetic resonance spectral data.
periodically tipping the magnetization vectors of said spin-spin coupled first and second groups of gyromagnetic resonators with respect to said polarizing magnetic field;
periodically flipping the magnetization vectors of said precessing resonators relative to the direction of the polarizing magnetic field for reversing the defocusing effect of the precession effect of the precession of said resonators about the polarizing magnetic field to obtain spin echo resonance of said resonators;
periodically detecting the spin echo resonance of said resonators; and changing the period of time t between tipping of said magnetization vectors and detection of the spin echo resonance for a series of said periodically detected spin echo resonances to obtain detected gyromagnetic resonance data about said first and second spin-spin coupled groups of resonators from which to derive simplified gyromagnetic resonance spectral data.
2. The method of Claim 1 wherein the step of periodi-cally detecting the spin echo resonance of said resonators includes the step of detecting respective ones of said spin echo resonances at a plurality m of time displaced intervals t2 to derive sets of spin echo resonance data as a function of n values of the time t1.
3. The method of Claim 2 including the step of Fourier transforming said sets of spin echo resonance data with respect to t2 from the time domain into the frequency domain to obtain sets of spin echo resonance data in the frequency domain .omega.2 as a function of thq time t1 and as a function of .omega.2, where .omega.2 is related to ?.
4. The method of Claim 3 including the step of Fourier transforming the Fourier transformed sets of spin echo resonance data in the frequency domain with respect to t1 into corresponding sets of spin echo resonances in the frequency domain as a function of the frequency .omega.1 which is related to the inverse of t1 to obtain sets of spin echo resonance spectral data as a function of both .omega.1 and .omega.2.
5. The method of Claim 4 including the steps of projecting said sets of spin echo resonance data as a function of both .omega.1 and .omega.2 onto an axis in the .omega.1, .omega.2 plane to obtain simplified gyromagnetic resonance spectral data.
6. The method of Claim 5 wherein said .omega.1 and .omega.2 axes are orthogonal and wherein said projection is onto the .omega.2 axis along a line at approximately 45° to said .omega.2 axis.
7. The method of Claim 5 wherein said .omega.1 and .omega.2 axes are orthogonal to each other and said projection is taken generally orthogonal to and onto said .omega.1 axis.
8. In a gyromagnetic resonance spectrometer, operative upon an assembly of first and second groups of gyromagnetic resonators immersed in a polarizing magnetic field wherein the magnetization vectors of said gyromagnetic resonators precess about said polarizing magnetic field, said first and second groups being coupled through spin-spin inter-action, said spin-spin interaction causing a defocusing effect whereby the phase coherence of precession among said resonators is progressively lost:
means for periodically tipping the magnetization vectors of said spin-spin coupled first and second groups of gyromagnetic resonators with respect to said polarizing magnetic field;
means for periodically flipping the magnetization vectors of said precessing resonators relative to the direction of the polarizing magnetic field for reversing the defocusing effect of the precession of said resonators about said polarizing magnetic field to obtain spin echo resona-nces of said resonators;
means for periodically detecting the spin echo resonance of said resonators; and means for changing the period of time t between tipping of said magnetization vectors and detection of the spin echo resonance for a series of said periodically detected spin echo resonances to obtain detected gyromagnetic resonance data about said first and second spin-spin coupled groups of resonators from which to derive simplified gyromagnetic resonance spectral data.
means for periodically tipping the magnetization vectors of said spin-spin coupled first and second groups of gyromagnetic resonators with respect to said polarizing magnetic field;
means for periodically flipping the magnetization vectors of said precessing resonators relative to the direction of the polarizing magnetic field for reversing the defocusing effect of the precession of said resonators about said polarizing magnetic field to obtain spin echo resona-nces of said resonators;
means for periodically detecting the spin echo resonance of said resonators; and means for changing the period of time t between tipping of said magnetization vectors and detection of the spin echo resonance for a series of said periodically detected spin echo resonances to obtain detected gyromagnetic resonance data about said first and second spin-spin coupled groups of resonators from which to derive simplified gyromagnetic resonance spectral data.
9. The apparatus of Claim 8 wherein said means for periodically detecting the spin echo resonance of said resonators includes means for detecting respective ones of said echo resonances at a plurality m of time displaced intervals of t2 to derive sets of spin echo resonance data as a function of n values of the time t1.
10. The apparatus of Claim 9 including means for Fourier transforming said sets of spin echo resonance data with respect to t2 from the time domain into the frequency domain to obtain sets of spin echo resonance data in the frequency domain .omega.2 as a function of the time t1 and as a function of .omega.2 where .omega.2 is related to ?.
11. The apparatus of Claim 10 including means for Fourier transforming the Fourier transform sets of spin echo resonance data in the frequency domain with respect to t1 in the corresponding sets of spin echo resonance data in the frequency domain as a function of the frequency .omega.1 which is related to the inverse of t1 to obtain sets of spin echo resonance spectral data as a function of both .omega.1 and .omega.2.
12. The apparatus of Claim 11 including means for projecting said sets of spin echo resonance data as a function of both .omega.1 and .omega.2 onto an axis in the .omega.1, .omega.2 plane to obtain simplified gyromagnetic resonance spectral data.
13. The apparatus of Claim 12 wherein said .omega.1 and .omega.2 axes are orthogonal and wherein said projection is taken onto the .omega.2 axis along a line at approximately 45° to said .omega.2 axis.
14. The apparatus of Claim 12 wherein said .omega.1 and .omega.2 axes are orthogonal to each other and said projection is taken generally orthogonal to and onto said .omega.1 axis.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US05/685,962 US4068161A (en) | 1976-05-13 | 1976-05-13 | Gyromagnetic resonance spectroscopy employing spin echo spin-spin decoupling and two-dimensional spreading |
Publications (1)
Publication Number | Publication Date |
---|---|
CA1067576A true CA1067576A (en) | 1979-12-04 |
Family
ID=24754371
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CA278,232A Expired CA1067576A (en) | 1976-05-13 | 1977-05-12 | Gyromagnetic resonance spectroscopy employing spin echo spin-spin decoupling and two-dimensional spreading |
Country Status (8)
Country | Link |
---|---|
US (1) | US4068161A (en) |
JP (2) | JPS5310490A (en) |
CA (1) | CA1067576A (en) |
CH (1) | CH616239A5 (en) |
DE (1) | DE2721011A1 (en) |
FR (1) | FR2351413A1 (en) |
GB (1) | GB1566481A (en) |
NL (1) | NL7705359A (en) |
Families Citing this family (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4068161A (en) * | 1976-05-13 | 1978-01-10 | Varian Associates, Inc. | Gyromagnetic resonance spectroscopy employing spin echo spin-spin decoupling and two-dimensional spreading |
US4168462A (en) * | 1977-10-20 | 1979-09-18 | Varian Associates, Inc. | Indirect detection of magnetic resonance by heteronuclear two-dimensional spectroscopy |
US4134058A (en) * | 1977-11-28 | 1979-01-09 | Varian Associates, Inc. | Selective detection of multiple quantum transitions in nuclear magnetic resonance |
US4297637A (en) * | 1978-07-20 | 1981-10-27 | The Regents Of The University Of California | Method and apparatus for mapping lines of nuclear density within an object using nuclear magnetic resonance |
US4238735A (en) * | 1979-02-21 | 1980-12-09 | Varian Associates, Inc. | Indirect detection of nuclear spins of low gyromagentic ratio coupled to spins of high gyromagnetic ratio |
US4379262A (en) * | 1979-08-10 | 1983-04-05 | Picker International Limited | Nuclear magnetic resonance systems |
US4291271A (en) * | 1979-11-01 | 1981-09-22 | Phillips Petroleum Company | Method for determining pore size distribution and fluid distribution in porous media |
US4345207A (en) * | 1980-04-24 | 1982-08-17 | The United States Of America As Represented By The Secretary Of The Navy | Method and apparatus for obtaining enhanced NMR signals |
US4383219A (en) * | 1980-08-19 | 1983-05-10 | Indiana University Foundation | Nuclear magnetic resonance spatial mapping |
US4443761A (en) * | 1981-06-16 | 1984-04-17 | National Research Development Corporation | NMR Spectroscopy |
JPS58142251A (en) * | 1982-02-19 | 1983-08-24 | Jeol Ltd | Measuring method for nuclear magnetic resonance |
US4521732A (en) * | 1982-03-16 | 1985-06-04 | Pegg David T | Pulse sequence for use in performing nuclear magnetic resonance spectroscopy |
US4543529A (en) * | 1983-01-04 | 1985-09-24 | Nmr Imaging, Inc. | Method of tuning an NMR apparatus |
US4502007A (en) * | 1983-01-31 | 1985-02-26 | Mee Gary L | Method of measuring NMR spin-spin relaxation time (T2) using spin-echos |
EP0155978B1 (en) * | 1984-03-29 | 1988-10-12 | Oxford Research Systems Limited | Method of operating a nuclear magnetic resonance spectrometer |
ATE82398T1 (en) * | 1984-04-05 | 1992-11-15 | Varian Associates | MAGNETIC NUCLEAR RESONANCE PULSE SEQUENCE FOR SPATIAL SELECTIVITY. |
US4766382A (en) * | 1984-05-17 | 1988-08-23 | Jeol Ltd. | Two-dimensional nuclear magnetic resonance spectrometry |
FI81204C (en) * | 1984-12-20 | 1990-09-10 | Instrumentarium Oy | Procedure for mapping the material properties of object s if to be investigated |
US4661775A (en) * | 1985-07-15 | 1987-04-28 | Technicare Corporation | Chemical shift imaging with field inhomogeneity correction |
US4721911A (en) * | 1985-07-26 | 1988-01-26 | Siemens Aktiengesellschaft | Nuclear magnetic resonance tomography apparatus |
US4701708A (en) * | 1986-08-01 | 1987-10-20 | General Electric Company | Polarization transfer by selective homonuclear technique for suppression of uncoupled spins in NMR spectroscopy |
DE3837317A1 (en) * | 1988-11-03 | 1990-05-10 | Philips Patentverwaltung | NUCLEAR RESONANCE SPECTROSCOPY METHOD AND ARRANGEMENT FOR IMPLEMENTING THE METHOD |
CN103645453B (en) * | 2013-12-23 | 2016-03-09 | 厦门大学 | A kind of method obtaining the monomer element one dimension localization spectrum eliminating scalar coupling modulation |
US10014906B2 (en) * | 2015-09-25 | 2018-07-03 | Microsemi Semiconductor (U.S.) Inc. | Acoustic echo path change detection apparatus and method |
CN106093099B (en) * | 2016-06-06 | 2018-06-29 | 厦门大学 | A kind of method for obtaining high-resolution two dimension J and decomposing spectrum |
CN107144591B (en) * | 2017-06-23 | 2018-11-30 | 厦门大学 | A method of measurement is independently with nuclear spin to indirect coupling mode |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3529235A (en) * | 1968-07-10 | 1970-09-15 | E H Research Lab Inc | Magnetic induction spectrometer employing a pair of coupled resonant cavities |
US3651396A (en) * | 1969-09-04 | 1972-03-21 | Bell Telephone Labor Inc | Fourier transform nuclear magnetic resonance spectroscopy |
US3648156A (en) * | 1970-05-26 | 1972-03-07 | Perkin Elmer Ltd | Quiet slot scanning |
US3753081A (en) * | 1971-12-30 | 1973-08-14 | Varian Associates | Gyromagnetic resonance method and apparatus for obtaining spin-spin coupling constants |
US3968424A (en) * | 1974-08-01 | 1976-07-06 | Varian Associates | Fourier transform NMR spectroscopy employing a phase modulated rf carrier |
US4068161A (en) * | 1976-05-13 | 1978-01-10 | Varian Associates, Inc. | Gyromagnetic resonance spectroscopy employing spin echo spin-spin decoupling and two-dimensional spreading |
-
1976
- 1976-05-13 US US05/685,962 patent/US4068161A/en not_active Expired - Lifetime
-
1977
- 1977-05-10 DE DE19772721011 patent/DE2721011A1/en not_active Ceased
- 1977-05-11 GB GB19742/77A patent/GB1566481A/en not_active Expired
- 1977-05-12 FR FR7714532A patent/FR2351413A1/en not_active Withdrawn
- 1977-05-12 CA CA278,232A patent/CA1067576A/en not_active Expired
- 1977-05-13 NL NL7705359A patent/NL7705359A/en not_active Application Discontinuation
- 1977-05-13 JP JP5447377A patent/JPS5310490A/en active Granted
- 1977-05-13 CH CH603977A patent/CH616239A5/de not_active IP Right Cessation
-
1985
- 1985-12-26 JP JP60292378A patent/JPS61262643A/en active Granted
Also Published As
Publication number | Publication date |
---|---|
JPH0115812B2 (en) | 1989-03-20 |
JPS61262643A (en) | 1986-11-20 |
JPS5310490A (en) | 1978-01-30 |
GB1566481A (en) | 1980-04-30 |
US4068161A (en) | 1978-01-10 |
FR2351413A1 (en) | 1977-12-09 |
CH616239A5 (en) | 1980-03-14 |
DE2721011A1 (en) | 1977-12-01 |
NL7705359A (en) | 1977-11-15 |
JPH0236901B2 (en) | 1990-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA1067576A (en) | Gyromagnetic resonance spectroscopy employing spin echo spin-spin decoupling and two-dimensional spreading | |
US4318043A (en) | Method and apparatus for rapid NMR imaging of nuclear densities within an object | |
US4431968A (en) | Method of three-dimensional NMR imaging using selective excitation | |
CA1052861A (en) | Gyromagnetic resonance fourier transform zeugmatography | |
US4297637A (en) | Method and apparatus for mapping lines of nuclear density within an object using nuclear magnetic resonance | |
US4689567A (en) | NMR Fourier imaging from multiple echoes | |
US5343151A (en) | Method for automatically shimming a high resolution NMR magnet | |
US4649346A (en) | Complex quotient nuclear magnetic resonance imaging | |
JPS58151545A (en) | Method of forming nuclear magnetic resonance picture | |
US4689566A (en) | NMR phase encoding using phase varying rf pulses | |
JPS6193941A (en) | Two-dimensional magnetic rotary resonance spectroscopic analyzer | |
JPH0332756B2 (en) | ||
US4567440A (en) | Vivo P-31 NMR imaging of phosphorus metabolites | |
Bodenhause et al. | Separation of dipolar and quadrupolar splittings in single crystal nitrogen-14 NMR | |
Emsley | [5] Selective pulses and their applications to assignment and structure determination in nuclear magnetic resonance | |
JPH07265281A (en) | Mr imaging system | |
US4814709A (en) | Method of selecting specific region of sample | |
EP0541636B1 (en) | Improvements in, or relating to magnetic resonance spectroscopy and imaging | |
US5317262A (en) | Single shot magnetic resonance method to measure diffusion, flow and/or motion | |
US5168229A (en) | Multidimensional nmr spectroscopy using switched acquisition time gradients for multiple coherence transfer pathway detection | |
US5262723A (en) | Method and apparatus for obtaining pure-absorption two-dimensional lineshape data for multidimensional NMR spectroscopy using switched acquisition time gradients | |
US4833408A (en) | NMR chemical shift imaging method with influence of intensity distribution of static magnetic field removed | |
US3753081A (en) | Gyromagnetic resonance method and apparatus for obtaining spin-spin coupling constants | |
US4672319A (en) | Multiple pulse excitation in NMR imaging | |
JPS62139641A (en) | Nmr imaging apparatus |