Radiofrequency pulse design in parallel transmission under strict temperature constraints

Actual patient position versus safety models: Specific Absorption Rate implications of initial head position for Ultrahigh Field Magnetic Resonance Imaging

Specific absorption rate (SAR) relates power absorption to tissue heating, and therefore is used as a safety constraint in magnetic resonance imaging (MRI). This study investigates the implications of initial head positioning on local and whole-head SAR. A virtual body model was simulated at 161 positions inside an eight-channel parallel-transmit (pTx) array. On-axis displacements and rotations of up to 20 mm/degrees and off-axis axial/coronal translations were investigated. Single-channel, radiofrequency (RF) shimming (i.e., single-spoke pTx) and multispoke pTx pulses were designed for seven axial, five coronal and five sagittal slices at each position (the slices were consistent across all positions). Whole-head and local SAR were calculated using safety models consisting of a single (centred) body position, multiple representative positions and all simulated body positions. Positional mismatches between safety models and actual positions cause SAR underestimation. For axial imaging, the actual peak local SAR was up to 4.2-fold higher for both single-channel and 5-spoke pTx, 3.5-fold higher for 3-/4-spoke pTx, and 2-fold higher for RF shimming and 2-spoke pTx, compared with that calculated using the centred body position. For sagittal and coronal imaging, the underestimation of peak local SAR was up to 5.2-fold and 3.8-fold, respectively. Using all body positions to estimate SAR prevented SAR underestimation but yielded up to 11-fold SAR overestimation for RF shimming. Local SAR of single-channel and pTx multispoke pulses showed considerable dependence on the initial patient position. RF shimming yielded much lower sensitivity to positional mismatches for axial imaging but not for sagittal and coronal imaging. This was deemed attributable to the higher degrees-of-freedom of control offered by the investigated coil array for axial imaging. Whole-head SAR is less sensitive to positional mismatches compared with local SAR. Nevertheless, whole-head SAR increased by up to 80% for sagittal imaging. Local and whole-head SAR were observed to be more sensitive to positional mismatches in the axial plane, because of larger variations in coil-tissue proximity. Using all possible body positions in the safety model may become substantially over-conservative and limit imaging performance, especially for the RF shimming mode for axial imaging.

Specific absorption rate implications of within-scan patient head motion for ultra-high field MRI

PURPOSE

This study investigates the implications of all degrees of freedom of within-scan patient head motion on patient safety.

METHODS

Electromagnetic simulations were performed by displacing and/or rotating a virtual body model inside an 8-channel transmit array to simulate 6 degrees of freedom of motion. Rotations of up to 20° and displacements of up to 20 mm including off-axis axial/coronal translations were investigated, yielding 104 head positions. Quadrature excitation, RF shimming, and multi-spoke parallel-transmit excitation pulses were designed for axial slice-selection at 7T, for seven slices across the head. Variation of whole-head specific absorption rate (SAR) and 10-g averaged local SAR of the designed pulses, as well as the change in the maximum eigenvalue (worst-case pulse) were investigated by comparing off-center positions to the central position.

RESULTS

In their respective worst-cases, patient motion increased the eigenvalue-based local SAR by 42%, whole-head SAR by 60%, and the 10-g averaged local SAR by 210%. Local SAR was observed to be more sensitive to displacements along right-left and anterior-posterior directions than displacement in the superior-inferior direction and rotation.

CONCLUSION

This is the first study to investigate the effect of all 6 degrees of freedom of motion on safety of practical pulses. Although the results agree with the literature for overlapping cases, the results demonstrate higher increases (up to 3.1-fold) in local SAR for off-axis displacement in the axial plane, which had received less attention in the literature. This increase in local SAR could potentially affect the local SAR compliance of subjects, unless realistic within-scan patient motion is taken into account during pulse design.

Parallel transmission RF pulse design with strict temperature constraints

RF safety in parallel transmission (pTx) is generally ensured by imposing specific absorption rate (SAR) limits during pTx RF pulse design. There is increasing interest in using temperature to ensure safety in MRI. In this work, we present a local temperature correlation matrix formalism and apply it to impose strict constraints on maximum absolute temperature in pTx RF pulse design for head and hip regions. Electromagnetic field simulations were performed on the head and hip of virtual body models. Temperature correlation matrices were calculated for four different exposure durations ranging between 6 and 24 min using simulated fields and body-specific constants. Parallel transmission RF pulses were designed using either SAR or temperature constraints, and compared with each other and unconstrained RF pulse design in terms of excitation fidelity and safety. The use of temperature correlation matrices resulted in better excitation fidelity compared with the use of SAR in parallel transmission RF pulse design (for the 6 min exposure period, 8.8% versus 21.0% for the head and 28.0% versus 32.2% for the hip region). As RF exposure duration increases (from 6 min to 24 min), the benefit of using temperature correlation matrices on RF pulse design diminishes. However, the safety of the subject is always guaranteed (the maximum temperature was equal to 39°C). This trend was observed in both head and hip regions, where the perfusion rates are very different.

SAR Simulations & Safety

At ultra-high fields, the assessment of radiofrequency (RF) safety presents several new challenges compared to low-field systems. Multi-channel RF transmit coils in combination with parallel transmit techniques produce time-dependent and spatially varying power loss densities in the tissue. Further, in ultra-high-field systems, localized field effects can be more pronounced due to a transition from the quasi stationary to the electromagnetic field regime. Consequently, local information on the RF field is required for reliable RF safety assessment as well as for monitoring of RF exposure during MR examinations. Numerical RF and thermal simulations for realistic exposure scenarios with anatomical body models are currently the only practical way to obtain the requisite local information on magnetic and electric field distributions as well as tissue temperature. In this article, safety regulations and the fundamental characteristics of RF field distributions in ultra-high-field systems are reviewed. Numerical methods for computation of RF fields as well as typical requirements for the analysis of realistic multi-channel RF exposure scenarios including anatomical body models are highlighted. In recent years, computation of the local tissue temperature has become of increasing interest, since a more accurate safety assessment is expected because temperature is directly related to tissue damage. Regarding thermal simulation, bio-heat transfer models and approaches for taking into account the physiological response of the human body to RF exposure are discussed. In addition, suitable methods are presented to validate calculated RF and thermal results with measurements. Finally, the concept of generalized simulation-based specific absorption rate (SAR) matrix models is discussed. These models can be incorporated into local SAR monitoring in multi-channel MR systems and allow the design of RF pulses under constraints for local SAR.

Toward imaging the body at 10.5 tesla

PURPOSE

To explore the potential of performing body imaging at 10.5 Tesla (T) compared with 7.0T through evaluating the transmit/receive performance of similarly configured dipole antenna arrays.

METHODS

Fractionated dipole antenna elements for 10.5T body imaging were designed and evaluated using numerical simulations. Transmit performance of antenna arrays inside the prostate, kidneys and heart were investigated and compared with those at 7.0T using both phase-only radiofrequency (RF) shimming and multi-spoke pulses. Signal-to-noise ratio (SNR) comparisons were also performed. A 10-channel antenna array was constructed to image the abdomen of a swine at 10.5T. Numerical methods were validated with phantom studies at both field strengths.

RESULTS

Similar power efficiencies were observed inside target organs with phase-only shimming, but RF nonuniformity was significantly higher at 10.5T. Spokes RF pulses allowed similar transmit performance with accompanying local specific absorption rate increases of 25-90% compared with 7.0T. Relative SNR gains inside the target anatomies were calculated to be >two-fold higher at 10.5T, and 2.2-fold SNR gain was measured in a phantom. Gradient echo and fast spin echo imaging demonstrated the feasibility of body imaging at 10.5T with the designed array.

CONCLUSION

While comparable power efficiencies can be achieved using dipole antenna arrays with static shimming at 10.5T; increasing RF nonuniformities underscore the need for efficient, robust, and safe parallel transmission methods. Magn Reson Med 77:434-443, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Parallel transmission for ultrahigh-field imaging

The development of MRI systems operating at or above 7 T has provided researchers with a new window into the human body, yielding improved imaging speed, resolution and signal-to-noise ratio. In order to fully realise the potential of ultrahigh-field MRI, a range of technical hurdles must be overcome. The non-uniformity of the transmit field is one of such issues, as it leads to non-uniform images with spatially varying contrast. Parallel transmission (i.e. the use of multiple independent transmission channels) provides previously unavailable degrees of freedom that allow full spatial and temporal control of the radiofrequency (RF) fields. This review discusses the many ways in which these degrees of freedom can be used, ranging from making more uniform transmit fields to the design of subject-tailored RF pulses for both uniform excitation and spatial selection, and also the control of the specific absorption rate. © 2015 The Authors. NMR in Biomedicine published by John Wiley & Sons Ltd.

Subject- and resource-specific monitoring and proactive management of parallel radiofrequency transmission

PURPOSE

Develop a practical comprehensive package for proactive management of parallel radiofrequency (RF) transmission.

METHODS

With a constrained optimization framework and predictive models from a prescan based multichannel calibration, we presented a method supporting design and optimization of parallel RF excitation pulses that accurately obey the forward/reflected peak and average power limits of the RF power amplifiers in parallel transmit imaging experiments and Bloch simulations. Moreover, local SAR limits were incorporated into the parallel RF excitation pulses using electromagnetic field simulations. Virtual transmit coils concept for minimization of reflected power (effecting subject-specific matching) was additionally demonstrated by leveraging experimentally calibrated power models.

RESULTS

Incorporation of experimentally calibrated power prediction models resulted in accurate compliance with prescribed hardware and global specific absorption rate (SAR) limits. Incorporation of spatial average 10 g SAR models, facilitated by simplifying numerical approximations, provided assurance of patient safety. RF pulses designed with various constraints demonstrated excellent excitation fidelity-the normalized root-mean-square error of the simulated excitation profiles was 2.6% for the fully constrained pulses, comparable to that of the unconstrained pulses. An RF shimming example showed a reduction of the reflected-to-forward power ratio to 1.7% from a conventional approach's 8.1%.

CONCLUSION

Using the presented RF pulse design method, effective proactive management of the multifaceted power and SAR limits was demonstrated in experimental and simulation studies. Magn Reson Med 76:20-31, 2016. © 2015 Wiley Periodicals, Inc.

Direct control of the temperature rise in parallel transmission by means of temperature virtual observation points: Simulations at 10.5 Tesla

PURPOSE

A method using parallel transmission to mitigate B1+ inhomogeneity while explicitly constraining the temperature rise is reported and compared with a more traditional SAR-constrained pulse design.

METHODS

Finite difference time domain simulations are performed on a numerical human head model and for a 16-channel coil at 10.5 Tesla. Based on a set of presimulations, a virtual observation point compression model for the temperature rise is derived. This compact representation is then used in a nonlinear programming algorithm for pulse design under explicit temperature rise constraints.

RESULTS

In the example of a time-of-flight sequence, radiofrequency pulse performance in some cases is increased by a factor of two compared with SAR-constrained pulses, while temperature rise is directly and efficiently controlled. Pulse performance can be gained by relaxing the SAR constraints, but at the expense of a loss of direct control on temperature.

CONCLUSION

Given the importance of accurate safety control at ultrahigh field and the lack of direct correspondence between SAR and temperature, this work motivates the need for thorough thermal studies in normal in vivo conditions. The tools presented here will possibly contribute to safer and more efficient MR exams.

Predicting long-term temperature increase for time-dependent SAR levels with a single short-term temperature response

PURPOSE

Present a novel method for rapid prediction of temperature in vivo for a series of pulse sequences with differing levels and distributions of specific energy absorption rate (SAR).

THEORY AND METHODS

After the temperature response to a brief period of heating is characterized, a rapid estimate of temperature during a series of periods at different heating levels is made using a linear heat equation and impulse-response (IR) concepts. Here the initial characterization and long-term prediction for a complete spine exam are made with the Pennes' bioheat equation where, at first, core body temperature is allowed to increase and local perfusion is not. Then corrections through time allowing variation in local perfusion are introduced.

RESULTS

The fast IR-based method predicted maximum temperature increase within 1% of that with a full finite difference simulation, but required less than 3.5% of the computation time. Even higher accelerations are possible depending on the time step size chosen, with loss in temporal resolution. Correction for temperature-dependent perfusion requires negligible additional time and can be adjusted to be more or less conservative than the corresponding finite difference simulation.

CONCLUSION

With appropriate methods, it is possible to rapidly predict temperature increase throughout the body for actual MR examinations.

A generalized slab-wise framework for parallel transmit multiband RF pulse design

PURPOSE

We propose a new slab-wise framework to design parallel transmit multiband pulses for volumetric simultaneous multislice imaging with a large field of view along the slice direction (FOVs).

THEORY AND METHODS

The slab-wise framework divides FOVs into a few contiguous slabs and optimizes pulses for each slab. Effects of relevant design parameters including slab number and transmit B1 (B1+) mapping slice placement were investigated for human brain imaging by designing pulses with global or local SAR control based on electromagnetic simulations of a 7T head RF array. Pulse design using in vivo B1+ maps was demonstrated and evaluated with Bloch simulations.

RESULTS

RF performance with respect to SAR reduction or B1+ homogenization across the entire human brain improved with increasing slabs; however, this improvement was nonlinear and leveled off at ∼12 slabs when the slab thickness reduced to ∼12 mm. The impact of using different slice placements for B1+ mapping was small.

CONCLUSION

Compared with slice-wise approaches where each of the many imaging slices requires both B1+ mapping and pulse optimization, the proposed slab-wise design framework attained comparable RF performance while drastically reducing the number of required pulses; therefore, it can be used to increase time efficiency for B1+ mapping, pulse calculation, and sequence preparation.