Abstract:
A method for performing magnetic resonance measurements on a sample includes applying a first predetermined number of pulse trains for excitation, each pulse train having a constant amplitude and including a second predetermined number of pulses spaced by a predetermined time interval. The pulse trains are modulated by a bent function. After each pulse, data is sampled. Preferably a square number of pulses is generated being constant in power, and the Walsh transform of the sequence of their phases is constant in power, so that the power of the excitation in time and frequency domain is constant. The method can reduce power requirements and may permit undercutting specific absorption rate (SAR) limits due to the small excitation power necessary to create time signals with reasonable signal to noise ratio.

Description:
BACKGROUND 
       [0001]    The invention relates to performing magnetic resonance measurements such as spectroscopic NMR (nuclear magnetic resonance) measurements as well as in MRI (magnetic resonance imaging). 
         [0002]    In NMR experiments spins in a sample are aligned by a static magnetic field. The aligned spins are then excited by applying electromagnetic pulses to the sample. This can be done in immobile spectrometers which often comprise a superconducting magnet or with portable NMR devices which are e.g. known by the trademark NMR-Mouse in the market. In particular when imaging human beings it is important to limit the excitation power which is introduced by the electromagnetic pulses to the patient for to undercut specific absorption rate (SAR) limits. Furthermore, when performing NMR measurements with mobile NMR devices it is necessary to limit excitation power to increase battery life. Therefore, low excitation power is desirable. Nevertheless, in NMR techniques known from prior art the drawback of low excitation power is an insufficient response peak power which results in a low signal to noise ratio and necessitates to repeat the measurements frequently while summing up the results to increase signal to noise. Similar problems arise in other magnetic resonance methods as electron spin resonance (ESR) and nuclear quadrupole resonance (NQR). 
       SUMMARY 
       [0003]    It is desired to solve the problems known from prior art at least partly and in particular to provide a method for performing magnetic resonance measurements and a respective apparatus to allow measurements with a good signal to noise ratio at low excitation peak power and high response peak power. 
         [0004]    According to a disclosed method, magnetic resonance measurements are performed on a sample, in which a first predetermined number of pulse trains each comprising a second predetermined number of pulses spaced by a predetermined time interval are applied for excitation, each pulse having a constant amplitude. The pulse trains are modulated by a bent function (known in this context as Frank sequences as well), wherein after each pulse data is sampled. 
         [0005]    It is understood that the term magnetic resonance measurements covers at least Nuclear Magnetic Resonance (NMR) measurements, Nuclear Quadrupole Resonance (NQR) measurements, and Electron Spin Resonance (ESR) measurements. The term “the pulse trains are modulated by a bent function” is preferably understood such that a predetermined square number m 2  of pulses are generated which form the first predetermined number of pulse trains. These pulse trains are constant in power, and the Walsh transform of the sequence of their phases is constant in power, too. This means that the spectral power density of the sequence is constant. The term amplitude is understood as the sum of the squares of the real and imaginary components of the pulse trains in the rotating coordinate frame representation. 
         [0006]    It is preferred to chose a positive natural number m, and a sequence of numbers z(n), n=0, . . . , m 2 −1 is generated such that all z(n) have values ranging between 0, 1, . . . , m−1. 
         [0007]    The sequence of these numbers has the additional property that its Walsh transform, i.e. the sequence Z(N) defined as 
         [0000]    
       
         
           
             
               Z 
                
               
                 ( 
                 N 
                 ) 
               
             
             := 
             
               
                 ∑ 
                 
                   n 
                   = 
                   0 
                 
                 
                   
                     m 
                     2 
                   
                   - 
                   1 
                 
               
                
               
                 
                   exp 
                    
                   
                     ( 
                     
                       
                         
                           2 
                            
                           π 
                         
                         m 
                       
                        
                       
                         z 
                          
                         
                           ( 
                           n 
                           ) 
                         
                       
                     
                     ) 
                   
                 
                  
                 
                   exp 
                    
                   
                     ( 
                     
                       
                         - 
                         
                           
                             2 
                              
                             π 
                              
                             
                                 
                             
                              
                              
                           
                           
                             m 
                             2 
                           
                         
                       
                        
                       nN 
                     
                     ) 
                   
                 
               
             
           
         
       
     
         [0000]    for N=0, . . . , m 2 −1 has constant absolute value. Sequences that exhibit this property are called bent functions in this context. The term i denotes the imaginary unit, which means that the square of i equals −1. 
         [0008]    One example for the pulse trains is given by 
         [0000]        z ( n ):= n   1   n   0 , where  n=mn   1   +n   0    
         [0000]    which are the unique natural numbers in the range 0, . . . , m−1 that allow the m-ary representation of n=0, . . . , m 2 −1. 
         [0009]    z(n) as defined above can be varied e.g. in at least one of the following ways:
   (a) z 1 =z(n)+a, where a is an arbitrary number in 0, . . . , m−1;   (b) z 2 =z(n+a), where a an arbitrary number in 0, . . . , m 2 −1;   (c) z 3 =z(n·b), where b is a number between 0, . . . , m 2 −1 such that the greatest common divisor of b and m is 1; and   (d) z 4 =z(n)·b, where b is a number between 0, . . . , m−1 such that the greatest common divisor of b and m is one.   
 
         [0014]    It is to be understood that it is possible to combine the possibilities (a) to (d) in such way that the function z in any of the possibilities (a) to (d) can be replaced by any the function z 1 , z 2 , z 3  and z 4 . 
         [0015]    For the NMR measurement a sequence of m 2  pulses is generated from z in the following way: all pulses have the same and preferably very small flip angle θ and the pulse direction, i.e. the angle between the pulse axis and the x-axis is given by 
         [0000]    
       
         
           
             
               
                 φ 
                  
                 
                   ( 
                   n 
                   ) 
                 
               
               = 
               
                 
                   
                     2 
                      
                     π 
                   
                   m 
                 
                  
                 
                   z 
                    
                   
                     ( 
                     n 
                     ) 
                   
                 
               
             
             , 
             
               
                 for 
                  
                 
                     
                 
                  
                 n 
               
               = 
               0 
             
             , 
             … 
              
             
                 
             
             , 
             
               
                 m 
                 2 
               
               - 
               1. 
             
           
         
       
     
         [0016]    Using this pulse sequence for an NMR measurement allows for a homogeneous excitation of the frequencies within a frequency window determined by the pulse spacing in the time domain. This leads to the possibility to use low excitation power while having maximum response power and thus a good signal to noise ratio. Data acquisition is preferably done after each pulse applied to the probe. 
         [0017]    Other bent functions of length m 2  could be used for NMR measurements with the same benefit. An example for such a bent function comes from the mathematical theory of Gaussian sums and is given by the sequence z(n)=n 2 , n=0, . . . , m 2 −1. 
         [0018]    Likewise, a sequence of pulse trains with constant small flip angle θ is generated where 
         [0000]    
       
         
           
             
               
                 φ 
                  
                 
                   ( 
                   n 
                   ) 
                 
               
               = 
               
                 
                   
                     2 
                      
                     π 
                   
                   
                     m 
                     2 
                   
                 
                  
                 
                   z 
                    
                   
                     ( 
                     n 
                     ) 
                   
                 
               
             
             , 
             
               
                 for 
                  
                 
                     
                 
                  
                 n 
               
               = 
               0 
             
             , 
             … 
              
             
                 
             
             , 
             
               
                 m 
                 2 
               
               - 
               1. 
             
           
         
       
     
         [0019]    Therefore, the partition of pulse phases that this function generates is finer than that of the example given above. 
         [0020]    The fact that the Walsh transform of all these functions is constant in absolute value implies that all frequencies in the frequency domain are covered with equal intensity, particularly also the 0 frequency. This allows scanning a time grid wherein adjacent points are separated by the time interval. 
         [0021]    In one embodiment the bent function may be chosen such that the pulse trains comprise a series of the first predetermined number n 1  of wavelet functions, where each pulse train is defined by the second predetermined number n 0  of discrete points in time. 
         [0022]    Preferably the ranges m of the first number n 1  and the second number n 0  are identical. This results in one important example of bent functions that can be understood as wavelets in terms of discrete representations of rotating waves 
         [0000]        x ( t )=exp{ i 2π ft}   (1)
 
         [0000]    over finite durations of time t and having a frequency f. The term i denotes the imaginary unit, whereas the square of i is −1. In a discrete representation the frequencies f can be written as 
         [0000]        f=n   1 /( mΔt )  (2)
 
         [0000]    where m is an integer, Δt is the predetermined time interval, and n 1  is a counter value. The time t is an integer multiple of the predetermined time interval: 
         [0000]        t=n   0   Δt   (3)
 
         [0023]    If x(t) is understood as a discrete representation of a rotating wave then the integer n 1  counts the number of turns of the wave within a fixed time (m Δt) whereas the integer n 0  counts the time increments within a rotating wave of given frequency f. 
         [0024]    This means the pulse trains can be described as 
         [0000]        x ( nΔt )=exp{ i 2π n   1   n   0   /m}.   (4)
 
         [0000]      with 
         [0000]        n=mn   1   +n   0 . 
         [0000]      For 
         [0000]      0≦ n   1   ≦m− 1 and 0≦ n   0   ≦m− 1
 
         [0000]    the largest values of the first and the second number are both m so that the pulse trains are m wavelets each defined at m successive values of n 0  As the absolute value of x is always one, and the envelope of the wavelets is rectangular. A Fourier transformation of the pulse trains reveals that the resulting power spectrum is constant. 
         [0025]    Therefore, a homogeneous excitation of the frequencies within the frequency window may be possible. Therefore, by exciting the spins with a multitude of constant-amplitude, small flip angle pulses the phases of which are discrete approximations of helices constituting complex wavelets defined on a discrete time grid, it is possible to use minimum excitation power while having maximum response power and thus a good signal to noise ratio. Due to the intermittent excitation and sampling of data without longer delays for to await relaxation phenomena, fast data sampling is possible. In this context it is preferred to sample one data point after each pulse. 
         [0026]    The use of wavelet functions, in particular as described above, has been found to be advantageous as it is possible to measure results with sufficient resolution with reasonable limitations on the hardware used to generate the pulse trains and sample the respective data. In particular, constant amplitude radio frequency (rf) pulse trains can be generated with simple nonlinear electronic components including rf amplifiers, while the pulse phase can easily be set with high accuracy using modern rf electronics. 
         [0027]    The spectrum of excitation frequencies may cover zero and integer multiples of a frequency interval. 
         [0028]    This allows a well-defined coverage of the frequency domain. The frequency interval of one wavelet is the reciprocal value of the product of the predetermined number of pulse trains and the predetermined time interval. 
         [0029]    The time interval may be determined such that by each pulse the net magnetization of the sample is rotated by a small flip angle. 
         [0030]    It is understood that the term small flip angle covers flip angles of less than 0.5°, in particular less than 0.2°, in particular less than 0.1° or less than 0.01° or even less than 0.005°. The rotation of the net magnetization by a small flip angle allows the rapid sampling of data points. 
         [0031]    The method may allow for rapid sampling of spectra and/or of imaging data if a respective magnetic field gradient is present with low excitation power and good signal to noise ratio. 
         [0032]    A An apparatus for executing the disclosed method may include: 
         [0000]    a) a magnet arrangement for creating a static magnetic field in a sample;
 
b) a coil for exciting, refocusing and sampling a spin signal in the sample, and
 
c) a controller that causes the magnet arrangement and coil to apply the pulse trains as described above to the sample.
 
         [0033]    According to an improvement the magnet arrangement comprises at least one magnet of at least one of the following kinds of magnets: 
         [0000]    a) a superconducting magnet;
 
b) a permanent magnet; and
 
c) an electromagnet.
 
         [0034]    The apparatus and the method can be used in low-field NMR spectroscopy with NMR frequencies of 5 MHz or more, e.g. 10 MHz and more and in high-field NMR spectroscopy and imaging with NMR frequencies of hydrogen of 300 MHz (Megahertz) and more, in particular of 600 MHz and more. In particular in imaging measurements in static magnetic fields of 3 Tesla and more, in particular of 7 T and more, and in particular 11 T and more, it may be possible to undercut specific absorption rate (SAR) limits with reasonable signal to noise ratios. Therefore, high resolution imaging in particular of human beings is possible. 
         [0035]    According to a further improvement the magnet arrangement is a single-sided magnet arrangement. 
         [0036]    It is preferred that these single-sided magnet arrangements have a sufficiently homogeneous stray field, so that an impulse response can be observed in the form of a free induction decay. 
         [0037]    According to a further improvement the magnet arrangement is a small permanent magnet arrangement with a sufficiently homogeneous field in the bore. It is understood that a small permanent magnet arrangement has a diameter of at most 30 cm (centimeters), preferably at most 20 cm, in particular at most 10 cm or even at most 5 cm or 3 cm. 
         [0038]    The method and apparatus may allow the measurement of spectra with reasonable signal to noise ratio with an excitation power significantly lower than that of the pulse sequences known in the art of Fourier NMR. Therefore, in particular movable NMR tools or devices can be used with a reasonably low consumption of electric energy. Compared to prior art pulse sequences the method and apparatus may enable the reduction of the excitation power by a factor of at least the order of four and thus a reduction of the energy consumption by a factor of at least the order of two with comparable results. This allows a significant reduction of battery capacities needed for movable NMR tools. Furthermore, the excitation power is reduced. This is particularly important for NMR tools used in bio-sensing applications, in sensing and control of chemical processes, in bore-hole applications, in airport security applications such as NMR fluid screening and NQR explosives and drug detection, in conducting NMR at contaminated sites like e.g. in nuclear contaminated environments or the like, in single-sided NMR of large or precious objects, and in NMR imaging in an emergency vehicle. One preferred application is in bore-hole NMR e.g. analyzing fluids in porous rocks. In these cases the measurements are performed with mobile instrumentation underground and the energy supply is critical. 
         [0039]    With the method and apparatus it may be possible to provide NMR tools having a button cell as an energy source instead of the large and heavy battery packs currently required. 
         [0040]    According to a further improvement the apparatus comprises means for generating magnetic field gradients for spatially resolved sampling of NMR signals. 
         [0041]    The means for generating magnetic field gradients can preferably comprise one or more gradient coils which are designed to create defined magnetic gradient fields at least at the sensitive volume of the probe where the coil and the sample are situated. The term sample is understood to comprise the material which is to be measured and can e.g. be a liquid sample in a container, i.e. e.g. a glass tube or the like, a human body or a part of it and so forth. 
         [0042]    The method and apparatus can be used e.g. in NMR, NQR, MRI, and ESR measurements. One use for the NQR detection is to detect the presence of explosives and/or narcotics. The disclosed method and apparatus may also be used for borehole NMR measurements, well-logging measurements, and/or measurements below the earth surface in porous environments comprising fluids. The term “wells” is understood to include water and/or oil wells. The term “fluid” is understood to include gases and/or liquids. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWING 
         [0043]    Embodiments will now be described in detail by way of example only with reference to the accompanying drawings, in which the following is depicted schematically: 
           [0044]      FIG. 1  an embodiment of an NMR pulse sequence; 
           [0045]      FIG. 2  an example of pulses of one pulse train; 
           [0046]      FIG. 3  an example of an NMR signal and the respective NMR spectrum; 
           [0047]      FIG. 4  an embodiment of an apparatus; 
           [0048]      FIG. 5  a first example of a magnet arrangement which can be used in an apparatus; and 
           [0049]      FIG. 6  a second example of a magnet arrangement which can be used in an apparatus. 
       
    
    
     DETAILED DESCRIPTION 
       [0050]      FIG. 1  depicts schematically an NMR pulse train with 1024 data points. For the sake of clarity the pulse sequence is illustrated in  FIG. 1  in a two-dimensional fashion. According to the m-ary representation of n=mn 1 +n 0  the horizontal numerals refer to n 0  whereas the vertical numerals refer to n 1 . Both range from 0, . . . , 31, where for typesetting reasons the n 0  numerals are only partly furnished. 
         [0051]    The pulse trains  1  are shown with their respective real component  2  and the imaginary component  3  corresponding to the projections of the pulse trains  1  to the x and y directions in the rotating coordinate frame the frequency of which is defined by the so called Larmor frequency which is proportional to the static magnetic field at the sample. For sake of clarity the reference numerals are only depicted for one pulse train  1  (wavelet) with frequency n 1 =0 but not for the respective pulse trains  1  with n 1 =1 to 31 (according to equation (4) above). Therefore, both the first number of pulse trains  1  (wavelets) and the second number of pulses per sequence are 32. 32 time intervals  4  are depicted as being separated by vertical lines. 
         [0052]    In each time interval  4  one pulse is applied and one data point is sampled. In one excitation scan 1024 pulses are applied successively for excitation, and hence 1024 data points are sampled. For each of the 32 values of n 1 , pulse trains  1  in the form of wavelets according to equation (4) above are applied successively for excitation. Each pulse train  1  corresponding to one specific wavelet excites one frequency window determined by the time interval  4  in the time domain and the number m of pulse trains. Multiple pulse trains  1  are applied without a recycle delay. The pulses of all pulse trains  1  have constant amplitude defined as the sum of the squares of the respective real component  2  and the respective imaginary component  3 . One data point is sampled at the end of a time interval  4 . The pulse trains  1  are applied to the sample successively. First, the pulse train  1  being a pulse sequence of the respective pulses enumerated as 0 in  FIG. 1  corresponding to n 1 =0 in equation (4) is applied to the sample in time intervals with n 0 =0, . . . , 31. Next, the pulse train with n 1 =1, . . . , 31 are applied in time intervals each with n 0 =0, . . . , 31. In any combination of n 1  and n 0  one data point is sampled. Therefore, a total of 1024 (32 rows and 32 columns) data points is sampled. 
         [0053]      FIG. 2  discloses a part of the pulse sequence in  FIG. 1  in more detail. The real component  2  of one pulse train  1  is depicted in part. This pulse train  1  has the shape of one of the wavelets as discussed above. The pulse train  1  is split into a train of pulses  13 . In each time interval  4  one pulse  13  is applied and one data point  14  is sampled. The height of the respective pulse  13  is modulated according to the value of the pulse train  1  at the respective time. Each pulse  13  has a pulse length  26  being smaller than the time interval  4 . 
         [0054]      FIG. 3  depicts on the left hand side (a) a respective time signal  5  having real and imaginary components which was measured with the pulse sequence depicted in  FIG. 1 .  FIG. 3  depicts on the right hand side (b) the respective spectrum  6  gathered by a Fourier transformation of the time signal  5 . The spectrum  6  has a real and an imaginary component as well. Time signal  5  was measured from an ethanol sample with a 300 MHz ( 1 H (proton) frequency) spectrometer. The excitation power was 35 μW (microwatts). 
         [0055]      FIG. 4  depicts schematically an NMR apparatus  7 . The apparatus  7  comprises a magnet arrangement  8  which comprises e.g. superconducting magnets, electromagnets or permanent magnets. The magnet arrangement  8  can be such that the sensitive volume in which the sample  9  is inserted is generated within the magnet arrangement  8  or on one side of the magnet arrangement  8 . The magnet arrangement  8  creates a static magnetic field at the sensitive volume which has to be located inside the sample  9 . The apparatus  7  further comprises a coil  10  for exciting and sampling a spin signal in the sample  9 . The coil  10  can be a surface coil, a solenoid coil or any other coil that is able to generate and sample the respective magnetic fields at the position of the sample  9 . 
         [0056]    The apparatus  7  further comprises control means  11  which are e.g. included in a computer or electronic device. The control means  11  are used to generate a signal with the coil  10  at the position of the sample  9 . The control means  11  are able to control the frequency, the shape and amplitude of the respective signal. Furthermore, the control means  11  are able to sample the signal received by the coil  10  after the pulses of the pulse trains  1 . The coil  10  is shaped such that the sample  9  is within the so called sensitive volume of the coil  10  in the magnet arrangement  8 . 
         [0057]    The control means  11  are provided with electrical energy by a power source  12 . This power source  12  comprises e.g. for mobile applications batteries. It may be possible to reduce the capacity of the batteries significantly as by the excitation with pulse trains  1  modulated by a bent function and the intermittent data sampling in between pulses  13  of the pulse trains  1  only low excitation power and thus only a small battery capacity is necessary. In high field MRI it may be possible to undercut SAR due to the small excitation power necessary to generate time signals  5  with a reasonable signal to noise ratio. 
         [0058]      FIG. 5  depicts schematically a first example of a magnet arrangement  8  which can be used in an apparatus  7 . The magnet arrangement  8  is e.g. shaped from first magnets  15  of trapezoidal cross section and second magnets  16  of rectangular cross section. The arrows in the magnets  15 ,  16  depict schematically the magnetization of these magnets  15 ,  16 . 
         [0059]    The first magnets  15  and second magnets  16  form three ring magnets  17  aligned such that a bore  18  for the sample  9  is created. These first magnets  15  and second magnets  16  are designed and arranged such that the magnet arrangement  8  forms a Halbach configuration. This means the polarization of the first magnets  15  and second magnets  16  are arranged to generate a strong magnetic field in the bore  18  within the ring magnets  17  and a very small field outside the ring magnets  17 . This magnet arrangement  8  allows the generation of highly homogeneous fields within the ring magnets  17 . The second magnets  16  can be moved relatively to the second magnets  15  towards the bore  18  and away from the bore  18  to increase the homogeneity of the magnetic field within the bore  18 . 
         [0060]      FIG. 6  depicts a further example of a magnet arrangement  8  which can be used in an apparatus  7 . The magnet arrangement  8  generates a sensitive volume  19  for single-sided NMR. In comparison to the first example of the magnet arrangement  8  as shown in  FIG. 5  the sensitive volume  19  is in the second example on one side of the magnet arrangement  8  whereas the sensitive volume of the magnet arrangement  8  of  FIG. 5  is within the bore  18 . 
         [0061]    The magnet arrangement  8 , as depicted in  FIG. 6 , comprises two first permanent magnets  20  for providing a first magnetic field at the sensitive volume  19 . The two first permanent magnets  20  are arranged such that their polarity is antiparallel. The two first permanent magnets  20  are arranged on an iron yoke  21  connecting the two first permanent magnets  20 . Above the yoke  21  and in between the two first permanent magnets  20  the magnet arrangement  8  further comprises a first pair  22  of shim magnets, a second pair  23  of shim magnets, a third pair  24  of shim magnets and a fourth pair  25  of shim magnets. The details of how to shim the magnetic field within the sensitive volume  19  is disclosed in EP 1 944 617 A1 which is incorporated by reference. 
         [0062]    Furthermore, an RF coil  10  for excitation and detection of signals in and from the sensitive volume  19  is further situated on the side opposite to the yoke  4 . The coil  11  is in this case a surface coil. 
         [0063]    Preferably a square number m 2  of pulses is generated. These pulses are constant in power, and the Walsh transform of the sequence of their phases is constant in power, too. This means that the power of the excitation in time and frequency domain is constant. 
         [0064]    It may be possible to reduce the capacity of the power supply significantly as by the excitation with pulse trains shaped as bent functions and the intermittent data sampling in between the pulses of the pulse trains only necessitates low excitation power and thus e.g. only a small battery capacity as a power supply is necessary. In high field MRI it may be possible to undercut SAR due to the small excitation power necessary to generate time signals  5  with a reasonable signal to noise ratio. 
       REFERENCE NUMERALS 
       [0000]    
       
           1  pulse train 
           2  real component 
           3  imaginary component 
           4  time interval 
           5  time signal 
           6  spectrum 
           7  NMR apparatus 
           8  magnet arrangement 
           9  sample 
           10  coil 
           11  control means 
           12  power source 
           13  pulse 
           14  data point 
           15  first magnet 
           16  second magnet 
           17  ring magnet 
           18  bore 
           19  sensitive volume 
           20  first permanent magnet 
           21  second permanent magnet 
           22  first pair of shim magnets 
           23  second pair of shim magnets 
           24  third pair of shim magnets 
           25  fourth pair of shim magnets 
           26  pulse length