Rapid determination of the RF pulse flip angle and spin-lattice relaxation time for materials imaging

For samples with T1s longer than 10s, calibration of the RF probe and a measurement of T1 can be very time-consuming. A technique is proposed for use in imaging applications where one wishes to rapidly obtain information about the RF flip angle and sample T1 prior to imaging. The flip angle measurement time is less than 1s for a single scan. Prior knowledge of the RF flip angle is not required for the measurement of T1. The resulting time savings in measuring the values of flip angle and T1 are particularly significant in the case of samples with very long T1 and short T2*. An imaging extension of the technique provides RF flip angle mapping without the need for incrementing the pulse duration, i.e., RF mapping can be performed at fixed RF amplifier output.     

An improved 3-D Look--Locker imaging method for T(1) parameter estimation

The 3-D Look-Locker (LL) imaging method has been shown to be a highly efficient and accurate method for the volumetric mapping of the spin lattice relaxation time T(1). However, conventional 3-D LL imaging schemes are typically limited to small tip angle RF pulses (5 degrees ), thereby improving the SNR and the accuracy of the method. In phantom studies, a mean T(1) measurement accuracy of less than 2% (0.2-3.1%) using a tip angle of 10 degrees was obtained for a range of T(1) from approximately 300 to 1,700 ms with a measurement time increase of only 15%. This accuracy compares favorably with the conventional 3-D LL method that provided an accuracy between 2.2% and 7.3% using a 5 degrees flip angle.     

Single-slice mapping of ultrashort T(2)

In this communication we present a method for single-slice mapping of ultrashort transverse relaxation times T(2). The RF pulse sequence consists of a spin echo preparation of the magnetization followed by slice-selective ultrashort echo time (UTE) imaging with radial k-space sampling. In order to keep the minimum echo time as small as possible, avoid out-of-slice contamination and signal contamination due to unwanted echoes, the implemented pulse sequence employs a slice-selective 180° RF refocusing pulse and a 4-step phase cycle. The slice overlap of the two slice-selective RF pulses was investigated. An acceptable Gaussian slice profile could be achieved by adjusting the strength of the two slice-selection gradients. The method was tested on a short T(2) phantom consisting of an arrangement of a roll of adhesive tape, an eraser, a piece of modeling dough made of Plasticine?, and a 10% w/w agar gel. The T(2) measurements on the phantom revealed exponential signal decays for all samples with T(2)(adhesive tape)=(0.5 ± 0.1)ms, T(2)(eraser)=(2.33 ± 0.07)ms, T(2)(Plasticine?)=(2.8 ± 0.06)ms, and T(2)(10%agar)=(9.5 ± 0.83)ms. The T(2) values obtained by the mapping method show good agreement with the T(2) values obtained by a non-selective T(2) measurement. For all samples, except the adhesive tape, the effective transverse relaxation time T(2)(?) was significantly shorter than T(2). Depending on the scanner hardware the presented method allows mapping of T(2) down to a few hundreds of microseconds. Besides investigating material samples, the presented method can be used to study the rapidly decaying MR-signal from biological tissue (e.g.: bone, cartilage, and tendon) and quadrupolar nuclei (e.g.: (23)Na, (35)Cl, and (17)O).     

Multiband excitation pulses for hyperpolarized 13C dynamic chemical-shift imaging

Hyperpolarized 13C offers high signal-to-noise ratios for imaging metabolic activity in vivo, but care must be taken when designing pulse sequences because the magnetization cannot be recovered once it has decayed. It has a short lifetime, on the order of minutes, and gets used up by each RF excitation. In this paper, we present a new dynamic chemical-shift imaging method that uses specialized RF pulses designed to maintain most of the hyperpolarized substrate while providing adequate SNR for the metabolic products. These are multiband, variable flip angle, spectral-spatial RF pulses that use spectral selectivity to minimally excite the injected prepolarized 13C-pyruvate substrate. The metabolic products of lactate and alanine are excited with a larger flip angle to increase SNR. This excitation was followed by an RF amplitude insensitive double spin-echo and an echo-planar flyback spectral-spatial readout gradient. In vivo results in rats and mice are presented showing improvements over constant flip angle RF pulses. The metabolic products are observable for a longer window because the low pyruvate flip angle preserves magnetization, allowing for improved observation of spatially varying metabolic reactions.     

Slice profile effects in 2D slice-selective MRI of hyperpolarized nuclei

This work explores slice profile effects in 2D slice-selective gradient-echo MRI of hyperpolarized nuclei. Two different sequences were investigated: a Spoiled Gradient Echo sequence with variable flip angle (SPGR-VFA) and a balanced Steady-State Free Precession (SSFP) sequence. It is shown that in SPGR-VFA the distribution of flip angles across the slice present in any realistically shaped radiofrequency (RF) pulse leads to large excess signal from the slice edges in later RF views, which results in an undesired non-constant total transverse magnetization, potentially exceeding the initial value by almost 300% for the last RF pulse. A method to reduce this unwanted effect is demonstrated, based on dynamic scaling of the slice selection gradient. SSFP sequences with small to moderate flip angles (<40°) are also shown to preserve the slice profile better than the most commonly used SPGR sequence with constant flip angle (SPGR-CFA). For higher flip angles, the slice profile in SSFP evolves in a manner similar to SPGR-CFA, with depletion of polarization in the center of the slice.     

Magnetization recovery for signal enhancement: a fast imaging DEFT-based technique

This paper describes the development and application of a new fast MRI technique based on the DEFT principle. The sequence named MAgnetization RecoverY for Signal Enhancement (MARYSE) is composed of two completely symmetric gradient echoes separated by a 180 degrees refocusing pulse. The RF pulse scheme, 90 degrees x-180 degrees y-90 degrees -x enables restoration of the transverse magnetization along the longitudinal axis, and consequently artificially increases R1 relaxation rate. In this sequence, the period between the excitation pulse and the restoring pulse (Tem: transverse magnetization evolution time) is very short (< 10 ms). This makes possible a significant increase in signal-to-noise ratio, even with a relatively short repetition time (20 ms). Simulations were performed for different values of Tem and TR at definite T1 and T2 and for different values of T1 and T2 at constant Tem and TR. Relevant signal enhancement for species with long relaxation time constants as compared to classical gradient echo and fast spin-echo imaging was expected. In vitro studies on a fat/water phantom confirmed this simulation. Application of MARYSE to mouse brain imaging permitted to visualize almost completely cerebrospinal fluid of the ventricles, a signal usually partially saturated in fast gradient echo imaging.     

SAR reduced pulse sequences

Three techniques were considered for reducing the RF (radiofrequency) power deposition in the body while maintaining scan time efficiency: reducing the RF peak amplitude while increasing the pulse width, substituting gradient echoes for spin echoes, and reducing the flip angle of the phase reversal pulse. The use of gradient echoes was found to be the most efficient means to reduce the power delivered to the patient and to obtain rapid data acquisition. The effect upon SAR (specific absorption rate) and SNR (signal-to-noise ratio) was demonstrated on a phantom when the phase reversal pulse was reduced from the standard 180 degrees to 90 degrees. Data in the body indicated a fairly constant SNR down to a refocusing flip angle between 110 degrees and 135 degrees. An initial clinical evaluation was performed at three institutions using the method of reducing the flip angle of the phase reversal pulse. The scan with theta = 120 degrees was rated by readers in a blinded study as having acceptable diagnostic image quality while the 135 degrees scan had comparable image quality to a conventional 90 degrees - 180 degrees pulse sequence. The use of reduced phase reversal pulses was seen as an efficient protocol to obtain T1-weighted images at rapid data rates while reducing the power delivered to the body by about 40%.     

J pulses for multiplet-selective NMR

Exact product operator solutions have been obtained for the evolution of weakly coupled spin-(1/2) I(m)S(n) systems during arbitrary RF irradiation of one spin. These solutions, which completely characterize the nature of J-coupling modulation during RF pulses, show that significant exchange occurs between single-spin magnetization and two-spin product operator states when the RF field strength is comparable to the coupling. In particular, a long (t(p) = [2J](-1) s), low-power (B(1) = J/2 Hz), constant amplitude pulse applied on resonance to one spin in an IS system completely interconverts the spinstates S(z) <--> 2S(x)I(z) and S(x) <--> 2S(z)I(z) when the RF is applied to the S spins, and interconverts S(x) <--> 2S(y)I(y) in 100% yield when the RF is applied to the I spins. Thus, these "J pulses," which select a bandwidth approximately equal to J Hz, may replace any combination of a (2J)(-1) delay period and a consecutive hard 90 degrees pulse in any polarization transfer or multiple quantum sequence. Although these rectangular pulses are highly frequency selective, in general they increase the replaced (2J)(-1) period by only a modest 40%, a time saving of a factor of 5 compared to existing pulses exhibiting the same selectivity. In favorable cases, there is no increase in duration of a pulse sequence using a particular type of J pulse, the 90(J) variety, which accomplishes the third spin state transformation listed above. J pulses will be advantageous for systems subject to rapid signal loss from relaxation and more generally for the enhanced operation of pulse sequences via the use of J modulation during RF irradiation.     

纵向弛豫对核磁共振线型的影响

根据经典 Bloch方程的解析解以及考虑辐射阻尼效应的Bloch方程的数值解,通过解析分析和数值模拟,从理论上研究了在射频场扰动下以及在辐射阻尼效应的作用下纵向弛豫对核磁共振线型的影响.结果表明:①射频场的扰动和辐射阻尼效应将导致纵向磁化与横向磁化的耦合,从而使纵向弛豫对线型产生了一定的影响.②在射频场的扰动下,峰强和线宽分别为2M0sin(θ)T1T2/(T1+T2)和(T1+T2)/(2πT1T2),即纵向弛豫将使谱线的峰强增大、线宽变窄,且影响程度随着比值T2/T1的减小而增大,峰强最大可增加1倍而线宽最多可减小1/2.③在强辐射阻尼效应的作用下,纵向弛豫会使谱线的峰强降低,降低的幅度与扳转角θ以及比值T2/T1密切相关.当θ从0到3π/4时,降低的幅度均较小,只有当θ＞3π/4时,降低的幅度才开始逐渐变大,且当 θ接近π时,降低的幅度急剧增大.谱线峰强降低的幅度与T2/T1呈较严格的正比关系,即T1越接近T2,峰强下降得越显著.     

Spin-lattice relaxation and a fast T1-map acquisition method in MRI with transient-state magnetization

The magnetization under the spin-lattice relaxation and the nuclear magnetic resonance radiofrequency (RF) pulses is calculated for a signal RF pulse train and for a sequence of multiple RF pulse-trains. It is assumed that the transverse magnetization is zero when each RF pulse is applied. The result expressions can be grouped into two terms: a decay term, which is proportional to the initial magnetization M0, and a recovery term, which has no M0 dependence but strongly depends on the spin-lattice relaxation and the equilibrium magnetization Meq. In magnetic resonance pulse sequences using magnetization in transient state, the recovery term produces artifacts and can seriously degrade the function of the preparation sequence for slice selection, contrast weighting, phase encoding, etc. This work shows that the detrimental effect can be removed by signal averaging in an eliminative fashion. A novel fast data acquisition method for constructing the spin-lattice relaxation (T1) map is introduced. The method has two features: (i) By using eliminative averaging, the curve to fit the T1 value is a decay exponential function rather than a recovery one as in conventional techniques; therefore, the measurement of Meq is not required and the result is less susceptible to the accuracy of the inversion RF pulse. (ii) The decay exponential curve is sampled by using a sequence of multiple pulse-trains. An image is reconstructed from each train and represents a sample point of the curve. Hence a single imaging sequence can yield multiple sample points needed for fitting the T1 value in contrast to conventional techniques that require repeating the imaging sequence for various delay values but obtain only one sample point from each repetition.     

Generalized RF Pulse Train with Insensitivity to B 1 Inhomogeneity

Adiabatic inversion recovery radiofrequency (RF) pulse techniques are used to address B ₁ inhomogeneity; however, the specific absorption rates of these techniques are significantly higher than that of non-adiabatic RF pulse techniques. In addition, time efficiency is poorer because of the required longer inversion recovery time. Therefore, an RF pulse train with three subpulses was previously developed and reported. The purpose of this article was to generalize the RF pulse train for tissues with different T ₁ relaxation times and in a different application. The RF pulse train B ₁ insensitivities and frequency responses were calculated with different T ₁ relaxation times and different subpulse durations using the Bloch equation. The previously reported optimal flip angle (FA) combination was used. When using the optimal FA combination, the RF pulse train B ₁ insensitivity did not change even if the T ₁ relaxation times and the subpulse durations did change. In other words, the optimal FA combination does not require adjustments according to the T ₁ and subpulse duration. The RF pulse train frequency responses with these subpulses can be dramatically improved even if the inherent subpulse frequency response is poor. This finding will facilitate RF pulse train technique implementation on magnetic resonance imaging scanners.     

Diffusion in porous building materials with high internal magnetic field gradients

An algorithm to calculate NMR signals of a multi-echo pulse sequence with arbitrary position dependent B0 and B1 fields taking into account relaxation and spin-diffusion is presented. The multi-echo pulse sequence consists of an initial RF pulse ("90 degrees " RF pulse) and a series of L refocusing RF pulses with arbitrary phases and flip-angles. The calculation is exact and takes into account all the magnetization pathways that contribute to the signal on a predefined spatial grid. The theoretical prediction is verified experimentally using a high field NMR microscopy system. The algorithm was implemented in a simulation program in order to optimize the design of an inside-out MR intra-vascular catheter that is used for characterization of vessel wall tissue. Measured data obtained with the catheter are in good agreement with the theoretical prediction of the simulation.     

The flipped longitudinal polarization sequence

By flipping the longitudinal magnetization with a chain of 180° pulses it is possible to effectively restore the effects of relaxation so that the same longitudinal magnetization is periodically recovered. The pulse sequence for achieving this, called Flipped LOngitudinal Polarization (FLOP), can be incorporated into any pulse sequence whenever it is desired to stop the attenuation in longitudinal magnetization caused by relaxation. We illustrate its use for fast, single-shot measurements of the longitudinal relaxation time and for three-dimensional T₁ mapping.     

Selective storage of magnetization in strongly relaxing spin systems

We have investigated the use of a "split-sinc" RF pulse to selectively store magnetization from a selected region of a sample, for later recall in imaging or spectroscopy experiments. The pulse sequence is based on an original suggestion by Post et al. (West German Patent No. P3209263.6, 13 March 1982), later implemented by Aue et al. (J. Magn. Reson. 56, 350 (1984)). We have carried out detailed numerical calculations using the Bloch equations and show that this particular sequence is robust in the face of strong transverse relaxation, and we demonstrate its application to imaging of polymer samples in shearing and extensional flow cells.     

Encoding to the longitudinal magnetization for MR imaging and flow velocity mapping

Phase-encoding to the longitudinal magnetization is implemented by adding encoding gradient pulses in the evolution period tau of the NMR pulse sequence 90 degrees+x-tau-90 degrees-x. This work focuses on the effect of the spin-lattice relaxation and its removal and on the constraint that the 90 degrees-x pulse can only transform the phase of the transverse magnetization partially to the longitudinal magnetization. Theoretical analysis shows that the encoded phase information and the spin-lattice relaxation effect are separable and the latter is identical in each repetition in collecting phase-encoding data. Thus the relaxation effect can be eliminated by subtracting a second data set whose phase information is inverted or by alternating the polarity of the relaxation contribution. From data with partial phase information, Fourier-transform image reconstruction results in mirror aliasing in which the two halves of the Fourier spectrum of positive and negative coordinates overlap. Removal of mirror aliasing requires imaging data of the orthogonal component. Nevertheless mirror aliasing is not necessarily a problem, depending on the subject of study. Phase-encoding to the longitudinal magnetization for spatial MRI and flow velocity mapping are demonstrated using the rotating ultra-fast imaging sequence (RUFIS).     

In vivo isotropic 3D diffusion tensor mapping of the rat brain using diffusion-weighted 3D MP-RAGE MRI

The purpose of this study was to examine the potential of diffusion-weighted (DW) three-dimensional (3D) MP-RAGE MRI for diffusion-tensor mapping of the rat brain in vivo. A DW-3D-MP-RAGE (3D-DWI) sequence was implemented at 2.0 T using six gradient orientations and a b value of 1000 s/mm2. In this sequence, the preparation sequence with a "90 degrees RF-motion proving gradient (MPG): MPG-180 degrees RF-MPG-90 degrees RF" pulse train (DW driven equilibrium Fourier transform) was used to sensitize the magnetization to diffusion. A centric k-space acquisition order was necessary to minimize saturation effects (T1 contamination) from tissues with short relaxation time. The image matrix was 128x128x128 (interpolated from 64x64x64 acquisitions), which resulted in small isotropic DW image data (voxel size: 0.273x0.273x0.273 mm3). Moreover, 3D-DWI-derived maps of the fractional anisotropy (FA), relative anisotropy (RA) and main-diffusion direction were completely free of susceptibility-induced signal losses and geometric distortions. Two well-known commissural fibers, the corpus callosum and anterior commissure, were indicated and shown to be in agreement with the locations of these known stereotaxic atlases. The experiment took 45 min, and shorter times should be possible in clinical application. The 3D-DWI sequence allows for in vivo 3D diffusion-tensor mapping of the rat brain without motion artifacts and susceptibility to distortion.     

General parameter relations for the Shinnar-Le Roux pulse design algorithm

The magnetization ripple amplitudes from a pulse designed by the Shinnar-Le Roux algorithm are a non-linear function of the Shinnar-Le Roux A and B polynomial ripples. In this paper, the method of Pauly et al. [J. Pauly, P. Le Roux, D. Nishimura, A. Macovski, Parameter relations for the Shinnar-Le Roux selective excitation pulse design algorithm, IEEE Transactions on Medical Imaging 10 (1991) 56-65.] has been extended to derive more general parameter relations. These relations can be used for cases outside the five classes considered by Pauly et al., in particular excitation pulses for flip angles that are not small or 90 degrees. Use of the new relations, together with an iterative procedure to obtain polynomials with the specified ripples from the Parks-McClellan algorithm, are shown to give simulated slice profiles that have the desired ripple amplitudes.     

Slice-selective excitation with B₁⁺-insensitive composite pulses

Spatially selective excitation pulses have been designed to produce uniform flip angles in the presence of the RF and static field inhomogeneities typically encountered in MRI studies of the human brain at 7 T. Pulse designs are based upon non-selective, composite pulses numerically optimized for the desired performance over prescribed ranges of field inhomogeneities. The non-selective pulses are subsequently transformed into spatially selective pulses with the same field-insensitive properties through modification of the spectral composition of the individual sub-pulses which are then executed in conjunction with an oscillating gradient waveform. An in-depth analysis of the performance of these RF pulses is presented in terms of total pulse durations, slice profiles, linearity of in-slice magnetization phase, sensitivity to RF and static field variations, and signal loss due to T(2) effects. Both simulations and measurements in phantoms and in the human brain are used to evaluate pulses with nominal flip angles of 45° and 90°. Target slice thickness in all cases is 2mm. Results indicate that the described class of field-insensitive RF pulses is capable of improving flip-angle uniformity in 7 T human brain imaging. There appears to be a subset of pulses with durations ?10 ms for which non-linearities in the magnetization phase are minimal and signal loss due to T(2) decay is not prohibitive. Such pulses represent practical solutions for achieving uniform flip angles in the presence of the large field inhomogeneities common to high-field human imaging and help to better establish the performance limits of high-field imaging systems with single-channel transmission.     

Optimization of the ultrafast look-locker echo-planar imaging T1 mapping sequence

The Look–Locker echo-planar imaging (LL-EPI) sequence has been numerically optimized in terms of the signal-to-noise ratio in the measured value of T₁, for both single-shot (repetition time (T_R) = ∞), and dynamically repeated T₁ measurements. The sequence is optimized for the normal biologic range of T₁ (0.2 s to 2.0 s) and for a range of sequence parameters found on most magnetic resonance (MR) scanners. Both linearly and geometrically spaced magnetization sample pulse intervals were considered. For single-shot measurements, the sequence with 24 linearly spaced sample pulses, an inversion time of 0.01 s, an inter-sample pulse delay of 0.10 s, and a sample radiofrequency (RF) pulse flip angle of 25^o was found to be optimum. When the number of sample pulses was limited due to hardware limitations, different pulse sequence parameters were indicated. The optimization procedures used are appropriate for any single-shot T₁ mapping sequence variant and for any rapid T₁ mapping application. The use of an optimized Look–Locker echo-planar imaging sequence is demonstrated by an example of dynamic contrast-enhanced scanning in the brain using fast T₁ mapping.     

Dynamics of paramagnetic agents by off-resonance rotating frame technique

Off-resonance rotating frame technique offers a novel tool to explore the dynamics of paramagnetic agents at high magnetic fields (B0 > 3T). Based on the effect of paramagnetic relaxation enhancement in the off-resonance rotating frame, a new method is described here for determining the dynamics of paramagnetic ion chelates from the residual z-magnetizations of water protons. In this method, the dynamics of the chelates are identified by the difference magnetization profiles, which are the subtraction of the residual z-magnetization as a function of frequency offset obtained at two sets of RF amplitude omega(1) and pulse duration tau. The choices of omega(1) and tau are guided by a 2-D magnetization map that is created numerically by plotting the residual z-magnetization as a function of effective field angle theta and off-resonance pulse duration tau. From the region of magnetization map that is the most sensitive to the alteration of the paramagnetic relaxation enhancement efficiency R(1rho)/R1, the ratio of the off-resonance rotating frame relaxation rate constant R(1rho) verse the laboratory frame relaxation rate constant R(1), three types of difference magnetization profiles can be generated. The magnetization map and the difference magnetization profiles are correlated with the rotational correlation time tauR of Gd-DTPA through numerical simulations, and further validated by the experimental data for a series of macromolecule conjugated Gd-DTPA in aqueous solutions. Effects of hydration water number q, diffusion coefficient D, magnetic field strength B0 and multiple rotational correlation times are explored with the simulations of the magnetization map. This method not only provides a simple and reliable approach to determine the dynamics of paramagnetic labeling of molecular/cellular events at high magnetic fields, but also a new strategy for spectral editing in NMR/MRI based on the dynamics of paramagnetic labeling in vivo.