Fast-kz three-dimensional tailored radiofrequency pulse for reduced B1 inhomogeneity

Multiphoton parallel transmission (MP-pTx): Pulse design methods and numerical validation

PURPOSE

We propose and evaluate multiphoton parallel transmission (MP-pTx) to mitigate flip angle inhomogeneities in high-field MRI. MP-pTx is an excitation method that utilizes a single, conventional birdcage coil supplemented with low-frequency (kHz) irradiation from a multichannel shim array and/or gradient channels. SAR analysis is simplified to that of a conventional birdcage coil, because only the radiofrequency (RF) field from the birdcage coil produces significant SAR.

METHODS

MP-pTx employs an off-resonance RF pulse from a conventional birdcage coil supplemented with oscillatingz $$ z $$ -directed fields from a multichannel shim array and/or the gradient coils. We simulate the ability of MP-pTx to create uniform nonselective brain excitations at 7 T using realisticB 1 + $$ {\mathrm{B}}_1^{+} $$ andΔ B 0 $$ \Delta {\mathrm{B}}_0 $$ field maps. The RF, shim array, and gradient waveform's amplitudes and phases are optimized using a genetic algorithm followed by sequential quadratic programming.

RESULTS

A 1 ms MP-pTx excitation using a 32-channel shim array with current constrained to less than 50 Amp-turns reduced the transverse magnetization's normalized root-mean-squared error from 29% for a conventional birdcage excitation to 6.6% and was nearly 40% better than a 1 ms birdcage coil 5 kT-point excitation with optimized kT-point locations and comparable pulse power.

CONCLUSION

The MP-pTx method resembles conventional pTx in its goals and approach but replaces the parallel RF channels with cheaper, low-frequency shim channels. The method mitigates high-field flip angle inhomogeneities to a level better than 3 T CP-mode and comparable to 7 T pTx while retaining the straightforward SAR characteristics of conventional birdcage excitations, as low-frequency shim array fields produce negligible SAR.

Hybrid algorithms for SAR matrix compression and the impact of post-processing on SAR calculation complexity

PURPOSE

This study proposes faster virtual observation point (VOP) compression as well as post-processing algorithms for specific absorption rate (SAR) matrix compression. Furthermore, it shows the relation between the number of channels and the computational burden for VOP-based SAR calculation.

METHODS

The proposed new algorithms combine the respective benefits of two different criteria for determining upper boundedness of SAR matrices by the VOPs. Comparisons of the old and new algorithms are performed for head coil arrays with various channel counts. The new post-processing algorithm is used to post-process the VOP sets of nine arrays, and the number of VOPs for a fixed median relative overestimation is compared.

RESULTS

The new algorithms are faster than the old algorithms by a factor of two to more than 10. The compression efficiency (number of VOPs relative to initial number of SAR matrices) is identical. For a fixed median relative overestimation, the number of VOPs increases logarithmically with the number of RF coil channels when post-processing is applied.

CONCLUSION

The new algorithms are much faster than previous algorithms. Post-processing is very beneficial for online SAR supervision of MRI systems with high channel counts, since for a given number of VOPs the relative SAR overestimation can be lowered.

Universal modes: Calibration-free time-interleaved acquisition of modes

PURPOSE

To introduce universal modes by applying the universal pulse concept to time-interleaved acquisition of modes (TIAMO), thereby achieving calibration-freeB 1 + $$ {B}_1^{+} $$ inhomogeneity mitigation for body imaging at ultra-high fields.

METHODS

Two databases of different RF arrays were used to demonstrate the feasibility of universal modes. The first comprised 31 cardiac in vivo data sets acquired at 7T while the second consisted of 6 simulated 10.5T pelvic data sets. Subject-specific solutions and universal modes were computed and subsequently evaluated alongside predefined default modes. For the cardiac database, subdivision into subpopulations was investigated. The optimization was performed using least-squares (LS) TIAMO and acquisition modes optimized for refocused echoes (AMORE). Finally, universal modes based on simulated pelvis data were applied in vivo at 10.5T.

RESULTS

In all studied cases, the universal modes yield improvements over the predefined default modes of up to 51% (cardiac) and 30% (pelvic) in terms of median excitation error when using two modes. The subpopulation-specific cardiac solutions revealed a further improvement of universal modes at the expense of increased errors when applied outside the appropriate subpopulation. Direct application of simulation-based universal modes in vivo resulted in up to a 14% reduction in excitation error compared to default modes and up to a 34% reduction in peak 10 g local specific absorption rate (SAR) compared to subject-specific solutions.

CONCLUSIONS

Universal modes are feasible for calibration-freeB 1 + $$ {B}_1^{+} $$ inhomogeneity mitigation at ultra-high fields. In addition, simulation-based solutions can be applied directly in vivo, eliminating the need for large in vivo databases.

Bayesian optimization of gradient trajectory for parallel-transmit pulse design

PURPOSE

Spoke pulses improve excitation homogeneity in parallel-transmit MRI. We propose an efficient global optimization algorithm, Bayesian optimization of gradient trajectory (BOGAT), for single-slice and simultaneous multislice imaging.

THEORY AND METHODS

BOGAT adds an outer loop to optimize kT-space positions. For each position, the RF coefficients are optimized (e.g., with magnitude least squares) and the cost function evaluated. Bayesian optimization progressively estimates the cost function. It automatically chooses the kT-space positions to sample, to achieve fast convergence, often coming close to the globally optimal spoke positions. We investigated the typical features of spokes cost functions by a grid search with field maps comprising 85 slabs from 14 volunteers. We tested BOGAT in this database, and prospectively in a phantom and in vivo. We compared the vendor-provided Fourier transform approach with the same magnitude least squares RF optimizer.

RESULTS

The cost function is nonconvex and seen empirically to be piecewise smooth with discontinuities where the underlying RF optimum changes sharply. BOGAT converged to within 10% of the global minimum cost within 30 iterations in 93% of slices in our database. BOGAT achieved up to 56% lower flip angle RMS error (RMSE) or 55% lower pulse energy in phantoms versus the Fourier transform approach, and up to 30% lower RMSE and 29% lower energy in vivo with 7.8 s extra computation.

CONCLUSION

BOGAT efficiently estimated near-global optimum spoke positions for the two-spoke tests, reducing flip-angle RMSE and/or pulse energy in a computation time (˜10 s), which is suitable for online optimization.

Investigation of alternative RF power limit control methods for 0.5T, 1.5T, and 3T parallel transmission cardiac imaging: A simulation study

PURPOSE

To investigate safety and performance aspects of parallel-transmit (pTx) RF control-modes for a body coil atB 0 ≤ 3 T $$ {B}_0\le 3\mathrm{T} $$ .

METHODS

Electromagnetic simulations of 11 human voxel models in cardiac imaging position were conducted forB 0 = 0.5 T $$ {B}_0=0.5\mathrm{T} $$ ,1.5 T $$ 1.5\mathrm{T} $$ and3 T $$ 3\mathrm{T} $$ and a body coil with a configurable number of transmit channels (1, 2, 4, 8, 16). Three safety modes were considered: the 'SAR-controlled mode' (SCM), where specific absorption rate (SAR) is limited directly, a 'phase agnostic SAR-controlled mode' (PASCM), where phase information is neglected, and a 'power-controlled mode' (PCM), where the voltage amplitude for each channel is limited. For either mode, safety limits were established based on a set of 'anchor' simulations and then evaluated in 'target' simulations on previously unseen models. The comparison allowed to derive safety factors accounting for varying patient anatomies. All control modes were compared in terms of theB 1 + $$ {B}_1^{+} $$ amplitude and homogeneity they permit under their respective safety requirements.

RESULTS

Large safety factors (approximately five) are needed if only one or two anchor models are investigated but they shrink with increasing number of anchors. The achievableB 1 + $$ {B}_1^{+} $$ is highest for SCM but this advantage is reduced when the safety factor is included. PCM appears to be more robust against variations of subjects. PASCM performance is mostly in between SCM and PCM. Compared to standard circularly polarized (CP) excitation, pTx offers minorB 1 + $$ {B}_1^{+} $$ improvements if local SAR limits are always enforced.

CONCLUSION

PTx body coils can safely be used atB 0 ≤ 3 T $$ {B}_0\le 3\mathrm{T} $$ . Uncertainties in patient anatomy must be accounted for, however, by simulating many models.

Fast online spectral-spatial pulse design for subject-specific fat saturation in cervical spine and foot imaging at 1.5 T

OBJECTIVE

To compensate subject-specific field inhomogeneities and enhance fat pre-saturation with a fast online individual spectral-spatial (SPSP) single-channel pulse design.

METHODS

The RF shape is calculated online using subject-specific field maps and a predefined excitation k-space trajectory. Calculation acceleration options are explored to increase clinical viability. Four optimization configurations are compared to a standard Gaussian spectral selective pre-saturation pulse and to a Dixon acquisition using phantom and volunteer (N = 5) data at 1.5 T with a turbo spin echo (TSE) sequence. Measurements and simulations are conducted across various body parts and image orientations.

RESULTS

Phantom measurements demonstrate up to a 3.5-fold reduction in residual fat signal compared to Gaussian fat saturation. In vivo evaluations show improvements up to sixfold for dorsal subcutaneous fat in sagittal cervical spine acquisitions. The versatility of the tailored trajectory is confirmed through sagittal foot/ankle, coronal, and transversal cervical spine experiments. Additional measurements indicate that excitation field (B1) information can be disregarded at 1.5 T. Acceleration methods reduce computation time to a few seconds.

DISCUSSION

An individual pulse design that primarily compensates for main field (B0) inhomogeneities in fat pre-saturation is successfully implemented within an online "push-button" workflow. Both fat saturation homogeneity and the level of suppression are improved.

Velocity selective spin labeling using parallel transmission

PURPOSE

Ultra-high field (UHF) provides improved SNR which greatly benefits SNR starved imaging techniques such as perfusion imaging. However, transmit field (B1 + ) inhomogeneities commonly observed at UHF hinders the excitation uniformity. Here we show how replacing standard excitation pulses with parallel transmit pulses can improve efficiency of velocity selective labeling.

METHODS

The standard tip-down and tip-up excitation pulses found in a velocity selective preparation module were replaced with tailored non-selective kT -points pulse solutions. Bloch simulations and experimental validation on a custom-built flow phantom and in vivo was performed to evaluate different pulse configurations in circularly polarized mode (CP-mode) and parallel transmit (pTx) mode.

RESULTS

Tailored pTx pulses significantly improved velocity selective labeling fidelity and signal uniformity. The transverse magnetization normalized RMS error was reduced from 0.489 to 0.047 when compared to standard rectangular pulses played in CP-mode. Simulations showed that manipulation of time symmetry in the tailored pTx pulses is vital in minimizing residual magnetization. In addition, in vivo experiments achieved a 44% lower RF power output and a shorter pulse duration when compared to using adiabatic pulses in CP-mode.

CONCLUSION

Using tailored pTx pulses for excitation within a velocity selective labeling preparation mitigated transmit field artifacts and improved SNR and contrast fidelity. The improvement in labeling efficiency highlights the potential of using pTx to improve robustness and accessibility of flow-based sequences such as velocity selective spin labeling at ultra-high field.

Comparison of tight-fitting 7T parallel-transmit head array designs using excitation uniformity and local specific absorption rate metrics

PURPOSE

We model the performance of parallel transmission (pTx) arrays with 8, 16, 24, and 32 channels and varying loop sizes built on a close-fitting helmet for brain imaging at 7 T and compare their local specific absorption rate (SAR) and flip-angle performances to that of birdcage coil (used as a baseline) and cylindrical 8-channel and 16-channel pTx coils (single-row and dual-row).

METHODS

We use the co-simulation approach along with MATLAB scripting for batch-mode simulation of the coils. For each coil, we extracted B1 + maps and SAR matrices, which we compressed using the virtual observation points algorithm, and designed slice-selective RF shimming pTx pulses with multiple local SAR and peak power constraints to generate L-curves in the transverse, coronal, and sagittal orientations.

RESULTS

Helmet designs outperformed cylindrical pTx arrays at a constant number of channels in the flip-angle uniformity at a constant local SAR metric: up to 29% for 8-channel arrays, and up to 34% for 16-channel arrays, depending on the slice orientation. For all helmet arrays, increasing the loop diameter led to better local SAR versus flip-angle uniformity tradeoffs, although this effect was more pronounced for the 8-channel and 16-channel systems than the 24-channel and 32-channel systems, as the former have more limited degrees of freedom and therefore benefit more from loop-size optimization.

CONCLUSION

Helmet pTx arrays significantly outperformed cylindrical arrays with the same number of channels in local SAR and flip-angle uniformity metrics. This improvement was especially pronounced for non-transverse slice excitations. Loop diameter optimization for helmets appears to favor large loops, compatible with nearest-neighbor decoupling by overlap.

What can we gain from subpopulation universal pulses? A simulation-based study

PURPOSE

The aim of the study was to explore a novel methodology for designing universal pulses (UPs) that balances the benefits of a calibration-free approach with subject-specific online pulse design.

METHODS

The proposed method involves segmenting the population into subpopulations with variability in anatomical shapes and positions reduced to 75%, 50%, and 25% of their original values while keeping the mean values unchanged. An additional 25% extreme case with a large volume of interest and shifted position was included. For each group, a 5kT-points universal inversion pulse was designed and assessed by the normalized root mean square error (NRMSE) on the target longitudinal magnetization profile. The performance was compared to the conventional one-size-fits-all approach. A total of 132 electromagnetic simulations were executed to generate representative anatomies and specific absorption rate (SAR) distributions in a three-dimensional parameter space comprised of head breadth, head length, and Y-shift. The 99.9th percentile on the peak local SAR distribution was utilized to establish an intersubject variability safety margin.

RESULTS

UPs designed for subpopulations with decreased head shape and position variability reduced the anatomical safety margin by up to 20%. Furthermore, when a head was significantly different to the average case, the proposed approach improved the inversion homogeneity by up to 24%, compared to the conventional one-size-fits-all approach.

CONCLUSION

Subpopulation UPs present an opportunity to improve theB 1 + $$ {\mathrm{B}}_1^{+} $$ homogeneity and reduce anatomical SAR safety margins at 7T without additional acquisition time for calibration.

Whole-liver flip-angle shimming at 7 T using parallel-transmit kT -point pulses and Fourier phase-encoded DREAM B 1 + mapping

PURPOSE

To obtain homogeneous signal throughout the human liver at 7 T. Flip angle (FA) shimming in 7T whole-liver imaging was performed through parallel-transmit kT -point pulses based on subject-specific multichannel absoluteB 1 + $$ {\mathrm{B}}_1^{+} $$ maps from Fourier phase-encoded dual refocusing echo acquisition mode (PE-DREAM).

METHODS

The optimal number of Fourier phase-encoding steps for PE-DREAMB 1 + $$ {\mathrm{B}}_1^{+} $$ mapping was determined for a 7T eight-channel parallel-transmission system. FA shimming experiments were performed in the liver of 7 healthy subjects with varying body mass index. In these subjects, firstB 0 $$ {\mathrm{B}}_0 $$ shimming and Fourier PE-DREAMB 1 + $$ {\mathrm{B}}_1^{+} $$ mapping were performed. Subsequently, three small-flip-angle 3D gradient-echo scans were acquired, comparing a circularly polarized (CP) mode, a phase shim, and a kT -point pulse. Resulting homogeneity was assessed and compared with estimated FA maps and distributions.

RESULTS

Fourier PE-DREAM with 13 phase-encoding steps resulted in a good tradeoff betweenB 1 + $$ {\mathrm{B}}_1^{+} $$ accuracy and scan time. Lower coefficient of variation values (average [min-max] across subjects) of the estimated FA in the volume of interest were observed using kT -points (7.4 [6.6%-8.0%]), compared with phase shimming (18.8 [12.9%-23.4%], p < 0.001) and CP (43.2 [39.4%-47.1%], p < 0.001). kT -points delivered whole-liver images with the nominal FA and the highest degree of homogeneity. CP and phase shimming resulted in either inaccurate or imprecise FA distributions. Here, locations having suboptimal FA in the estimated FA maps corresponded to liver areas suffering from inconsistent signal intensity and T1 -weighting in the gradient-echo scans.

CONCLUSION

Homogeneous whole-liver 3D gradient-echo acquisitions at 7 T can be obtained with eight-channel kT -point pulses calculated based on subject-specific multichannel absolute Fourier PE-DREAMB 1 + $$ {\mathrm{B}}_1^{+} $$ maps.

Towards an integrated radiofrequency safety concept for implant carriers in MRI based on sensor-equipped implants and parallel transmission

To protect implant carriers in MRI from excessive radiofrequency (RF) heating it has previously been suggested to assess that hazard via sensors on the implant. Other work recommended parallel transmission (pTx) to actively mitigate implant-related heating. Here, both ideas are integrated into one comprehensive safety concept where native pTx safety (without implant) is ensured by state-of-the-art field simulations and the implant-specific hazard is quantified in situ using physical sensors. The concept is demonstrated by electromagnetic simulations performed on a human voxel model with a simplified spinal-cord implant in an eight-channel pTx body coil at 3 T . To integrate implant and native safety, the sensor signal must be calibrated in terms of an established safety metric (e.g., specific absorption rate [SAR]). Virtual experiments show that E -field and implant-current sensors are well suited for this purpose, while temperature sensors require some caution, and B 1 probes are inadequate. Based on an implant sensor matrix Q s , constructed in situ from sensor readings, and precomputed native SAR limits, a vector space of safe RF excitations is determined where both global (native) and local (implant-related) safety requirements are satisfied. Within this safe-excitation subspace, the solution with the best image quality in terms of B 1 + magnitude and homogeneity is then found by a straightforward optimization algorithm. In the investigated example, the optimized pTx shim provides a 3-fold higher mean B 1 + magnitude compared with circularly polarized excitation for a maximum implant-related temperature increase ∆ T imp ≤ 1 K . To date, sensor-equipped implants interfaced to a pTx scanner exist as demonstrator items in research labs, but commercial devices are not yet within sight. This paper aims to demonstrate the significant benefits of such an approach and how this could impact implant-related RF safety in MRI. Today, the responsibility for safe implant scanning lies with the implant manufacturer and the MRI operator; within the sensor concept, the MRI manufacturer would assume much of the operator's current responsibility.

Improvement in Intraoperative Image Quality in Transcranial Magnetic Resonance-Guided Focused Ultrasound Surgery Using Transmitter Gain Adjustment

INTRODUCTION

Transcranial magnetic resonance-guided focused ultrasound surgery (TcMRgFUS) has the advantage of allowing immediate evaluation of therapeutic effects after each sonication and intraoperative magnetic resonance imaging (MRI) to visualize the lesion. When the image shows that the lesion has missed the planned target and the therapeutic effects are insufficient, the target of the subsequent ablation can be finely adjusted based on the image. The precision of this adjustment is determined by the image quality. However, the current intraoperative image quality with a 3.0T MRI system is insufficient for precisely detecting the lesion. Thus, we developed and validated a method for improving intraoperative image quality.

METHODS

Because intraoperative image quality is affected by transmitter gain (TG), we acquired T2-weighted images (T2WIs) with two types of TG: the automatically adjusted TG (auto TG) and the manually adjusted TG (manual TG). To evaluate the character of images with 2 TGs, the actual flip angle (FA), the image uniformity, and the signal-to-noise ratio (SNR) were measured using a phantom. Then, to assess the quality of intraoperative images, T2WIs with both TGs were acquired during TcMRgFUS for 5 patients. The contrast-to-noise ratio (CNR) of the lesion was retrospectively estimated.

RESULTS

The images of the phantom with the auto TG showed substantial variations between the preset and actual FAs (p < 0.01), whereas on the images with the manual TG, there were no variations between the two FAs (p > 0.05). The total image uniformity was considerably lower with the manual TG than with the auto TG (p < 0.01), indicating that the image's signal values with the manual TG were more uniform. The manual TG produced significantly higher SNRs than the auto TG (p < 0.01). In the clinical study, the lesions were clearly detected in intraoperative images with the manual TG, but they were difficult to identify in images with the auto TG. The CNR of lesions in images with manual TG was considerably higher than in images with auto TG (p < 0.01).

CONCLUSION

Regarding intraoperative T2WIs using a 3.0T MRI system during TcMRgFUS, the manual TG method improved image quality and delineated the ablative lesion more clearly than the current method with auto TG.

Ultra-high field MRI: parallel-transmit arrays and RF pulse design

This paper reviews the field of multiple or parallel radiofrequency (RF) transmission for magnetic resonance imaging (MRI). Currently the use of ultra-high field (UHF) MRI at 7 tesla and above is gaining popularity, yet faces challenges with non-uniformity of the RF field and higher RF power deposition. Since its introduction in the early 2000s, parallel transmission (pTx) has been recognized as a powerful tool for accelerating spatially selective RF pulses and combating the challenges associated with RF inhomogeneity at UHF. We provide a survey of the types of dedicated RF coils used commonly for pTx and the important modeling of the coil behavior by electromagnetic (EM) field simulations. We also discuss the additional safety considerations involved with pTx such as the specific absorption rate (SAR) and how to manage them. We then describe the application of pTx with RF pulse design, including a practical guide to popular methods. Finally, we conclude with a description of the current and future prospects for pTx, particularly its potential for routine clinical use.

Parallel transmit pulse design for saturation homogeneity (PUSH) for magnetization transfer imaging at 7T

PURPOSE

This work proposes a novel RF pulse design for parallel transmit (pTx) systems to obtain uniform saturation of semisolid magnetization for magnetization transfer (MT) contrast in the presence of transmit field B1+ inhomogeneities. The semisolid magnetization is usually modeled as being purely longitudinal, with the applied B1+ field saturating but not rotating its magnetization; thus, standard pTx pulse design methods do not apply.

THEORY AND METHODS

Pulse design for saturation homogeneity (PUSH) optimizes pTx RF pulses by considering uniformity of root-mean squared B1+ , B1rms , which relates to the rate of semisolid saturation. Here we considered designs consisting of a small number of spatially non-selective sub-pulses optimized over either a single 2D plane or 3D. Simulations and in vivo experiments on a 7T Terra system with an 8-TX Nova head coil in five subjects were carried out to study the homogenization of B1rms and of the MT contrast by acquiring MT ratio maps.

RESULTS

Simulations and in vivo experiments showed up to six and two times more uniform B1rms compared to circular polarized (CP) mode for 2D and 3D optimizations, respectively. This translated into 4 and 1.25 times more uniform MT contrast, consistently for all subjects, where two sub-pulses were enough for the implementation and coil used.

CONCLUSION

The proposed PUSH method obtains more uniform and higher MT contrast than CP mode within the same specific absorption rate (SAR) budget.

Improved TSE imaging at ultrahigh field using nonlocalized efficiency RF shimming and acquisition modes optimized for refocused echoes (AMORE)

Purpose

To develop and evaluate a novel RF shimming optimization strategy tailored to improve the transmit efficiency in turbo spin echo imaging when performing time‐interleaved acquisition of modes (TIAMO) at ultrahigh fields.

Theory and Methods

A nonlocalized efficiency shimming cost function is proposed and extended to perform TIAMO using acquisition modes optimized for refocused echoes (AMORE). The nonlocalized efficiency shimming was demonstrated in brain and knee imaging at 7 Tesla. Phantom and in vivo torso imaging studies were performed to compare the performance between AMORE and previously proposed TIAMO mode optimizations with and without localized constraints in turbo spin echo and gradient echo acquisitions.

Results

The proposed nonlocalized efficiency RF shimming produced a circularly polarized‐like field with fewer signal dropouts in the brain and knee. For larger targets, AMORE was used and required a significantly lower transmitter voltage to produce a similar contrast to existing TIAMO mode design approaches for turbo spin echo as well as gradient echo acquisitions. In vivo, AMORE effectively reduced signal dropout in the interior torso while providing more uniform contrast with reduced transmit power. A local constraint further improved performance for a target region while maintaining performance in the larger FOV.

Conclusion

AMORE based on the presented nonlocalized efficiency shimming cost function demonstrated improved contrast and SNR uniformity as well as increased transmit efficiency for both gradient echo and turbo spin echo acquisitions.

Click here for author‐reader discussions

Performance and safety assessment of an integrated transmit array for body imaging at 7 T under consideration of specific absorption rate, tissue temperature, and thermal dose

In this study, the performance of an integrated body-imaging array for 7 T with 32 radiofrequency (RF) channels under consideration of local specific absorption rate (SAR), tissue temperature, and thermal dose limits was evaluated and the imaging performance was compared with a clinical 3 T body coil. Thirty-two transmit elements were placed in three rings between the bore liner and RF shield of the gradient coil. Slice-selective RF pulse optimizations for B₁ shimming and spokes were performed for differently oriented slices in the body under consideration of realistic constraints for power and local SAR. To improve the B₁⁺ homogeneity, safety assessments based on temperature and thermal dose were performed to possibly allow for higher input power for the pulse optimization than permissible with SAR limits. The results showed that using two spokes, the 7 T array outperformed the 3 T birdcage in all the considered regions of interest. However, a significantly higher SAR or lower duty cycle at 7 T is necessary in some cases to achieve similar B₁⁺ homogeneity as at 3 T. The homogeneity in up to 50 cm-long coronal slices can particularly benefit from the high RF shim performance provided by the 32 RF channels. The thermal dose approach increases the allowable input power and the corresponding local SAR, in one example up to 100 W/kg, without limiting the exposure time necessary for an MR examination. In conclusion, the integrated antenna array at 7 T enables a clinical workflow for body imaging and comparable imaging performance to a conventional 3 T clinical body coil.

B1 field map synthesis with generative deep learning used in the design of parallel-transmit RF pulses for ultra-high field MRI

Spoke trajectory parallel transmit (pTX) excitation in ultra-high field MRI enables B₁⁺ inhomogeneities arising from the shortened RF wavelength in biological tissue to be mitigated. To this end, current RF excitation pulse design algorithms either employ the acquisition of field maps with subsequent non-linear optimization or a universal approach applying robust pre-computed pulses. We suggest and evaluate an intermediate method that uses a subset of acquired field maps combined with generative machine learning models to reduce the pulse calibration time while offering more tailored excitation than robust pulses (RP). The possibility of employing image-to-image translation and semantic image synthesis machine learning models based on generative adversarial networks (GANs) to deduce the missing field maps is examined. Additionally, an RF pulse design that employs a predictive machine learning model to find solutions for the non-linear (two-spokes) pulse design problem is investigated. As a proof of concept, we present simulation results obtained with the suggested machine learning approaches that were trained on a limited data-set, acquired in vivo. The achieved excitation homogeneity based on a subset of half of the B₁⁺ maps acquired in the calibration scans and half of the B₁⁺ maps synthesized with GANs is comparable with state of the art pulse design methods when using the full set of calibration data while halving the total calibration time. By employing RP dictionaries or machine-learning RF pulse predictions, the total calibration time can be reduced significantly as these methods take only seconds or milliseconds per slice, respectively.

Post-processing algorithms for specific absorption rate compression

PURPOSE

Compression of local specific absorption rate (SAR) matrices is essential for enabling SAR monitoring and efficient pulse calculation in parallel transmission. Improvements in compression result in lower error margin and/or lower number of virtual observation points (VOPs). The purpose of this work is to introduce two algorithms for post-processing of already compressed VOP sets. One calculates individual overestimation matrices for the VOPs to reduce overestimation, the other identifies redundant VOPs.

METHODS

The first algorithm was evaluated for VOP sets calculated for three different transmit arrays with either 8 or 16 channels. For each array, two different overestimation matrices were used to generate the VOP sets. Each post-processed VOP set was evaluated using one million random excitation vectors and the results compared to the VOP set before post-processing. The second algorithm was evaluated by utilizing the same random excitation vectors and comparing the results after removal of the redundant VOPs with the results before removal to verify that these were identical.

RESULTS

The first algorithm reduced the mean overestimation by up to four fifths compared to the original set, while keeping the number of VOPs constant. The second algorithm decreased the number of VOPs generated by a compression with Eichfelder and Gebhardt's algorithm by more than 40% in 40% of the investigated cases and by more than 20% in 73% of the investigated cases.

CONCLUSION

Two post-processing algorithms are presented that enhance previously compressed VOP sets by improving the accuracy per number of VOPs.

Local excitation universal parallel transmit pulses at 9.4T

PURPOSE

To demonstrate that the concept of "universal pTx pulses" is applicable to local excitation applications.

METHODS

A database of B₀ / B 1 + maps from eight different subjects was acquired at 9.4T. Based on these maps, universal pulses that aim at local excitation of the visual cortex area in the human brain (with a flip angle of 90° or 7°) were calculated. The remaining brain regions should not experience any excitation. The pulses were designed with an extension of the "spatial domain method." A 2D and a 3D target excitation pattern were tested, respectively. The pulse performance was examined on non-database subjects by Bloch simulations and in vivo at 9.4T using a GRE anatomical MRI and a presaturated TurboFLASH B 1 + mapping sequence.

RESULTS

The calculated universal pulses show excellent performance in simulations and in vivo on subjects that were not contained in the design database. The visual cortex region is excited, while the desired non-excitation areas produce the only minimal signal. In simulations, the pulses with 3D target pattern show a lack of excitation uniformity in the visual cortex region; however, in vivo, this inhomogeneity can be deemed acceptable. A reduced field of view application of the universal pulse design concept was performed successfully.

CONCLUSIONS

The proposed design approach creates universal local excitation pulses for a flip angle of 7° and 90°, respectively. Providing universal pTx pulses for local excitation applications prospectively abandons the need for time-consuming subject-specific B₀ / B 1 + mapping and pTx-pulse calculation during the scan session.

Robust RF shimming and small-tip-angle multispoke pulse design with finite-difference regularization

PURPOSE

A new regularizer is proposed for the magnitude least-squares optimization algorithm, to ensure robust parallel transmit RF shimming and small-tip-angle multispoke pulse designs for ultrahigh-field MRI.

METHODS

A finite-difference regularization term is activated as an additional regularizer in the iterative magnitude-least-squares based pulse design algorithm when an unwanted flip angle null distribution is detected. Both simulated and experimental B 1 + maps from different transmit arrays and different human subjects at 7 T were used to evaluate the proposed algorithm. The algorithm was further demonstrated in experiment with dynamic multislice RF shimming for a single-shot gradient-echo EPI for human functional MRI at 7 T.

RESULTS

The proposed finite-difference regularizer effectively prevented excitation null to be formed for RF shimming and small-tip-angle multispoke pulses, and improved the latter with a monotonic trade-off relationship between flip angle error and RF power. The proposed algorithm was demonstrated to be effective with several head-array geometries by simulation and with a commercial head array with 12 healthy human subjects by experiment. During a functional MRI scan at 7 T with dynamic RF shimming, the proposed algorithm ensured high image SNR throughout the human brain, compared with near-complete local signal loss by the conventional magnitude-least-squares algorithm.

CONCLUSION

Using finite-difference regularization to avoid unwanted solutions, the robustness of RF shimming and small-tip-angle multispoke pulse design algorithms are improved, with better flip angle homogeneity and a monotonic trade-off relationship between flip angle error and RF power.

Calibration-free regional RF shims for MRS

PURPOSE

Achieving a desired RF transmit field ( B 1 + ) in small regions of interest is critical for single-voxel MRS at ultrahigh field. Radio-frequency (RF) shimming, using parallel transmission, requires B 1 + mapping and optimization, which limits its ease of use. This work aimed to generate calibration-free RF shims for predefined target regions of interest, which can be applied to any participant, to produce a desired absolute magnitude B 1 + (| B 1 + |).

METHODS

The RF shims were found offline by joint optimization on a database comprising B 1 + maps from 11 subjects, considering regions of interest in occipital cortex, hippocampus and posterior cingulate, as well as whole brain. The | B 1 + | achieved was compared with a tailored shimming approach, and MR spectra were acquired using tailored and calibration-free shims in 4 participants. Global and local 10g specific-absorption-rate deposition were estimated using Duke and Ella dielectric models.

RESULTS

There was no difference in the mean | B 1 + | produced using calibration-free versus tailored RF shimming in the occipital cortex (p = .15), hippocampus (p = .5), or posterior cingulate (p = .98), although differences were observed in the RMS error | B 1 + |. Spectra acquired using calibration-free shims had similar SNR and low residual water signal. Under identical power settings, specific-absorption-rate deposition was lower compared with operating in quadrature mode. For example, the total head specific absorption rate was around 35% less for the occipital cortex.

CONCLUSION

This work demonstrates that static RF shims, optimized offline for small regions, avoid the need for B 1 + mapping and optimization for each region of interest and participant. Furthermore, power settings may be increased when using calibration-free shims, to better take advantage of RF shimming.

Local SAR compression algorithm with improved compression, speed, and flexibility

PURPOSE

Local specific absorption rate (SAR) compression algorithms are essential for enabling online SAR monitoring in parallel transmission. A better compression resulting in a lower number of virtual observation points improves speed of SAR calculation for online supervision and pulse design.

METHOD

An iterative expansion of an existing algorithm presented by Lee et al is proposed in this work. The original algorithm is used within a loop, making use of the virtual observation points from the previous iteration as the starting subvolume, while decreasing the overestimation with each iteration. This algorithm is evaluated on the SAR matrices of three different simulated arrays.

RESULT

The number of virtual observation points is approximately halved with the new algorithm, while at the same time the compression time is reduced with speed-up factors of up to 2.5.

CONCLUSION

The new algorithm improves the original algorithm in terms of compression rate and speed.

Fast online-customized (FOCUS) parallel transmission pulses: A combination of universal pulses and individual optimization

PURPOSE

To mitigate spatial flip angle (FA) variations under strict specific absorption rate (SAR) constraints for ultra-high field MRI using a combination of universal parallel transmit (pTx) pulses and fast subject-specific optimization.

METHODS

Data sets consisting of B₀ , B 1 + maps, and virtual observation point (VOP) data were acquired from 72 subjects (study groups of 48/12 healthy Europeans/Asians and 12 Europeans with pathological or incidental findings) using an 8Tx/32Rx head coil on a 7T whole-body MR system. Combined optimization values (COV) were defined as combination of spiral-nonselective (SPINS) trajectory parameters and an energy regularization weight. A set of COV was optimized universally by simulating the individual RF pulse optimizations of 12 training data sets (healthy Europeans). Subsequently, corresponding universal pulses (UPs) were calculated. Using COV and UPs, individually optimized pulses (IOPs) were calculated during the sequence preparation phase (maximum 15 s). Two different UPs and IOPs were evaluated by calculating their normalized root-mean-square error (NRMSE) of the FA and SAR in simulations of all data sets. Seven additional subjects were examined using an MPRAGE sequence that uses the designed pTx excitation pulses and a conventional adiabatic inversion.

RESULTS

All pTx pulses resulted in decreased mean NRMSE compared to a circularly polarized (CP) pulse (CP = ~28%, UPs = ~17%, and IOPs = ~12%). UPs and IOPs improved homogeneity for all subjects. Differences in NRMSE between study groups were much lower than differences between different pulse types.

CONCLUSION

UPs can be used to generate fast online-customized (FOCUS) pulses gaining lower NRMSE and/or lower SAR values.

Application of Evolution Strategies to the Design of SAR Efficient Parallel Transmit Multi-Spoke Pulses for Ultra-High Field MRI

We present an evolution-strategy based approach to solve the magnitude least squares (MLS) design problem of low flip-angle slice-selective parallel transmit RF pulses for ultra-high field MRI using SAR and peak-RF-constraints. A combined transmit k-space trajectory and RF pulse weight optimization is proposed in two algorithmic steps. The first step is a coarse grid search to find an initial solution that fulfills all constraints for the subsequent multistage optimization. This avoids convergence to the next nearest local minimum. The second step attempts to refine the results using multiple evolution strategies. We compare the performance of our approach with the non-convex optimization methods described in the literature. The proposed algorithm converges for phantom and in vivo data and only requires an initial estimate of the range of suitable regularization parameters. It demonstrates improved excitation homogeneity compared to published spoke-design methods and allows optimization for homogeneity with a subsequent reduction in the SAR burden. Moreover, excitation homogeneity and the SAR burden can be balanced against each other, enabling a further reduction in SAR at the cost of minor relaxations in excitation homogeneity. This feature makes the algorithm a good candidate for SAR limited sequences in ultra-high field imaging. The algorithm is validated using phantom and in vivo measurements obtained with a 16-channel transmit array at 9.4T.

Improving FLAIR SAR efficiency at 7T by adaptive tailoring of adiabatic pulse power through deep learning B 1 + estimation

PURPOSE

The purpose of this study is to demonstrate a method for specific absorption rate (SAR) reduction for 2D T₂ -FLAIR MRI sequences at 7 T by predicting the required adiabatic radiofrequency (RF) pulse power and scaling the RF amplitude in a slice-wise fashion.

METHODS

We used a time-resampled frequency-offset corrected inversion (TR-FOCI) adiabatic pulse for spin inversion in a T₂ -FLAIR sequence to improve B 1 + homogeneity and calculated the pulse power required for adiabaticity slice-by-slice to minimize the SAR. Drawing on the implicit B 1 + inhomogeneity in a standard localizer scan, we acquired 3D AutoAlign localizers and SA2RAGE B 1 + maps in 28 volunteers. Then, we trained a convolutional neural network (CNN) to estimate the B 1 + profile from the localizers and calculated pulse scale factors for each slice. We assessed the predicted B 1 + profiles and the effect of scaled pulse amplitudes on the FLAIR inversion efficiency in oblique transverse, sagittal, and coronal orientations.

RESULTS

The predicted B 1 + amplitude maps matched the measured ones with a mean difference of 9.5% across all slices and participants. The slice-by-slice scaling of the TR-FOCI inversion pulse was most effective in oblique transverse orientation and resulted in a 1 min and 30 s reduction in SAR induced delay time while delivering identical image quality.

CONCLUSION

We propose a SAR reduction technique based on the estimation of B 1 + profiles from standard localizer scans using a CNN and show that scaling the inversion pulse power slice-by-slice for FLAIR sequences at 7T reduces SAR and scan time without compromising image quality.

Elimination of low-inversion-efficiency induced artifacts in whole-brain MP2RAGE using multiple RF-shim configurations at 7 T

The magnetization-prepared two-rapid-gradient-echo (MP2RAGE) sequence is used for structural T₁ -weighted imaging and T₁ mapping of the human brain. In this sequence, adiabatic inversion RF pulses are commonly used, which require the B₁⁺ magnitude to be above a certain threshold. Achieving this threshold in the whole brain may not be possible at ultra-high fields because of the short RF wavelength. This results in low-inversion regions especially in the inferior brain (eg cerebellum and temporal lobes), which is reflected as regions of bright signal in MP2RAGE images. This study aims at eliminating the low-inversion-efficiency induced artifacts in MP2RAGE images at 7 T. The proposed technique takes advantage of parallel RF transmission systems by splitting the brain into two overlapping slabs and calculating the complex weights of transmit channels (ie RF shims) on these slabs for excitation and inversion independently. RF shims were calculated using fast methods implemented in the standard workflow. The excitation RF pulse was designed to obtain slabs with flat plateaus and sharp edges. These slabs were joined into a single volume during the online image reconstruction. The two-slab strategy naturally results in a signal-to-noise ratio loss; however, it allowed the use of independent shims to make the B₁⁺ field exceed the adiabatic threshold in the inferior brain, eliminating regions of low inversion efficiency. Accordingly, the normalized root-mean-square errors in the inversion were reduced to below 2%. The two-slab strategy was found to outperform subject-specific k_T -point inversion RF pulses in terms of inversion error. The proposed strategy is a simple yet effective method to eliminate low-inversion-efficiency artifacts; consequently, MP2RAGE-based, artifact-free T₁ -weighted structural images were obtained in the whole brain at 7 T.

B0 and B1 inhomogeneities in the liver at 1.5 T and 3.0 T

PURPOSE

The purpose of this work is to characterize the magnitude and variability of B₀ and B₁ inhomogeneities in the liver in large cohorts of patients at both 1.5 T and 3.0 T.

METHODS

Volumetric B₀ and B₁ maps were acquired over the liver of patients presenting for routine abdominal MRI. Regions of interest were drawn in the nine Couinaud segments of the liver, and the average value was recorded. Magnitude and variation of measured averages in each segment were reported across all patients.

RESULTS

A total of 316 B₀ maps and 314 B₁ maps, acquired at 1.5 T and 3.0 T on a variety of GE Healthcare MRI systems in 630 unique exams, were identified, analyzed, and, in the interest of reproducible research, de-identified and made public. Measured B₀ inhomogeneities ranged (5th-95th percentiles) from -31.7 Hz to 164.0 Hz for 3.0 T (-14.5 Hz to 81.3 Hz at 1.5 T), while measured B₁ inhomogeneities (ratio of actual over prescribed flip angle) ranged from 0.59 to 1.13 for 3.0 T (0.83 to 1.11 at 1.5 T).

CONCLUSION

This study provides robust characterization of B₀ and B₁ inhomogeneities in the liver to guide the development of imaging applications and protocols. Field strength, bore diameter, and sex were determined to be statistically significant effects for both B₀ and B₁ uniformity. Typical clinical liver imaging at 3.0 T should expect B₀ inhomogeneities ranging from approximately -100 Hz to 250 Hz (-50 Hz to 150 Hz at 1.5 T) and B₁ inhomogeneities ranging from approximately 0.4 to 1.3 (0.7 to 1.2 at 1.5 T).

Local SAR compression with overestimation control to reduce maximum relative SAR overestimation and improve multi-channel RF array performance

Objective

In local SAR compression algorithms, the overestimation is generally not linearly dependent on actual local SAR. This can lead to large relative overestimation at low actual SAR values, unnecessarily constraining transmit array performance.

Method

Two strategies are proposed to reduce maximum relative overestimation for a given number of VOPs. The first strategy uses an overestimation matrix that roughly approximates actual local SAR; the second strategy uses a small set of pre-calculated VOPs as the overestimation term for the compression.

Result

Comparison with a previous method shows that for a given maximum relative overestimation the number of VOPs can be reduced by around 20% at the cost of a higher absolute overestimation at high actual local SAR values.

Conclusion

The proposed strategies outperform a previously published strategy and can improve the SAR compression where maximum relative overestimation constrains the performance of parallel transmission.

General gradient delay correction method in bipolar multispoke RF pulses using trim blips

PURPOSE

To correct with gradient trim blips for gradient delays in bipolar-spoke RF pulses in slice-selective and slab-selective excitations, compatible with tilted acquisitions and anisotropic delays.

THEORY

The effect of small gradient delays with respect to RF pulses results in a dephasing of the second RF spoke, proportional to the slab-selection gradient amplitude and the distance of the slice center from the magnet isocenter. Accordingly, adding a trim blip along the corresponding logical gradient axis between the two spokes compensates for the same dephasing, and therefore cancels the gradient delay effects, regardless of position and orientation.

METHODS

Gradient delays on different axes were first measured on a 7T system based on an imaging method. Parallel transmission universal bipolar spokes were designed offline to mitigate the RF field inhomogeneity problem in the human brain. Trim blips were used to compensate for the known delays, which was validated with flip angle and temporal SNR measurements on two different volunteers at 7 T.

RESULTS

Pulses corrected with trim blips greatly reduced gradient delay effects. Acquisitions made with corrected and noncorrected pulses showed good fidelity with simulations.

CONCLUSIONS

Unlike time or phase-shifting approaches, trim blip-based methods apply to all possible bipolar spoke scenarios such as slice excitations, slab excitations, and anisotropy in the gradient delays.

Universal nonselective excitation and refocusing pulses with improved robustness to off-resonance for Magnetic Resonance Imaging at 7 Tesla with parallel transmission

PURPOSE

In MRI at ultra-high field, the k T -point and spiral nonselective (SPINS) pulse design techniques can be advantageously combined with the parallel transmission (pTX) and universal pulse techniques to create uniform excitation in a calibration-free manner. However, in these approaches, pulse duration is typically increased as compared to standard hard pulses, and excitation quality in regions exhibiting large resonance frequency offsets often suffer. This limitation is inherent to structure of k T -point or SPINS pulse, and likely can be mitigated using parameterization-free pulse design approaches.

METHODS

The Gradient Ascent Pulse Engineering (GRAPE) algorithm was used to design parameterization-free RF and magnetic field gradient (MFG) waveforms for creating 8 ∘ excitation, up to 105 ∘ scalable refocusing and inversion, nonselectively across the brain. Simulations were performed to provide flip angle normalized root-mean-squares error (FA-NRMSE) estimations for the 8 ∘ and the 180 ∘ k T -point, SPINS, and GRAPE pulses. GRAPE pulses were tested experimentally with anatomical head scans at 7T.

RESULTS

As compared to k T -points and SPINS, GRAPE provided substantial improvement of excitation, refocusing, and inversion quality at off-resonance while at least preserving the same global FA-NRMSE performance. As compared to k T -points, GRAPE allowed for a substantial reduction of the pulse duration for the 8 ∘ excitation and the 105 ∘ refocusing.

CONCLUSIONS

Parameterization-free universal nonselective pTX-pulses were successfully computed using GRAPE. Performance gains as compared to k T -points were validated numerically and experimentally for three imaging protocols. In its current implementation, the computational burden of GRAPE limits its use to applications where pulse computations are not subject to time constraints.

Comparison of SMS-EPI and 3D-EPI at 7T in an fMRI localizer study with matched spatiotemporal resolution and homogenized excitation profiles

The simultaneous multi-slice EPI (SMS-EPI, a.k.a. MB-EPI) sequence has met immense popularity recently in functional neuroimaging. A still less common alternative is the use of 3D-EPI, which offers similar acceleration capabilities. The aim of this work was to compare the SMS-EPI and the 3D-EPI sequences in terms of sampling strategies for the detection of task-evoked activations at 7T using detection theory. To this end, the spatial and temporal resolutions of the sequences were matched (1.6 mm isotropic resolution, TR = 1200 ms) and their excitation profiles were homogenized by means of calibration-free parallel-transmission (Universal Pulses). We used a fast-event “localizer” paradigm of 5:20 min in order to probe sensorimotor functions (visual, auditory and motor tasks) as well as higher level functions (language comprehension, mental calculation), where results from a previous large-scale study at 3T (N = 81) served as ground-truth reference for the brain areas implicated in each cognitive function. In the current study, ten subjects were scanned while their activation maps were generated for each cognitive function with the GLM analysis. The SMS-EPI and 3D-EPI sequences were compared in terms of raw tSNR, t-score testing for the mean signal, activation strength and accuracy of the robust sensorimotor functions. To this end, the sensitivity and specificity of these contrasts were computed by comparing their activation maps to the reference brain areas obtained in the 3T study. Estimated flip angle distributions in the brain reported a normalized root mean square deviation from the target value below 10% for both sequences. The analysis of the t-score testing for the mean signal revealed temporal noise correlations, suggesting the use of this metric instead of the traditional tSNR for testing fMRI sequences. The SMS-EPI and 3D-EPI thereby yielded similar performance from a detection theory perspective.

Channel-combination method for phase-based |B1+| mapping techniques

PURPOSE

The aim of this study was to propose a channel combination method for |B₁⁺| mapping methods using phase difference to reconstruct |B₁⁺| map.

THEORY AND METHODS

Phase-based |B₁⁺| mapping methods commonly consider the phase difference of two scans to measure |B₁⁺|. Multiple receiver coils acquire a number of images and the phase difference at each channel is theoretically the same in the absence of noise. Affected by noise, phase difference is approximately governed by Gaussian distribution. Considering data from all channels as samples, estimation can be achieved by maximum likelihood method. With this method, all phase differences at each channel are combined into one. In this study, the proposed method is applied with Bloch-Siegert shift |B₁⁺| mapping method. Simulations are performed to illustrate the phase difference distribution and demonstrate the feasibility and facility of the proposed method. Phantom and vivo experiments are carried out at 1.5 T scanner equipped with 8-channel receiver coil. In all experiments, the proposed method is compared with weighted averaging (WA) method.

RESULTS

Simulations revealed appropriateness of approximating the distribution of phase difference to Gaussian distribution. Compared with WA method, the proposed method reduces errors of |B₁⁺| calculation. Phantom and vivo experiments provide further validation.

CONCLUSION

Considering phase noise distribution, the proposed method achieves channel combination by finding the estimation from data acquired by multiple receivers coil. The proposed method reduces |B₁⁺| reconstruction errors caused by noise.

A 32-channel parallel transmit system add-on for 7T MRI

PURPOSE

A 32-channel parallel transmit (pTx) add-on for 7 Tesla whole-body imaging is presented. First results are shown for phantom and in-vivo imaging.

METHODS

The add-on system consists of a large number of hardware components, including modulators, amplifiers, SAR supervision, peripheral devices, a control computer, and an integrated 32-channel transmit/receive body array. B1+ maps in a phantom as well as B1+ maps and structural images in large volunteers are acquired to demonstrate the functionality of the system. EM simulations are used to ensure safe operation.

RESULTS

Good agreement between simulation and experiment is shown. Phantom and in-vivo acquisitions show a field of view of up to 50 cm in z-direction. Selective excitation with 100 kHz sampling rate is possible. The add-on system does not affect the quality of the original single-channel system.

CONCLUSION

The presented 32-channel parallel transmit system shows promising performance for ultra-high field whole-body imaging.

IMPULSE: A scalable algorithm for design of minimum specific absorption rate parallel transmit RF pulses

PURPOSE

Managing local specific absorption rate (SAR) in parallel transmission requires ensuring that the peak SAR over a large number of voxels (> 10 5 ) is below the regulatory limit. The safety risk to the patient depends on cumulative (not instantaneous) SAR thus making a joint design of all RF pulses in a sequence desirable. We propose the Iterative Minimization Procedure with Uncompressed Local SAR Estimate (IMPULSE), an efficient optimization formulation and algorithm that can handle uncompressed SAR matrices and optimize pulses for all slices jointly within a practical time frame.

THEORY AND METHODS

IMPULSE optimizes parallel transmit pulses for small-tip-angle slice selective excitation to minimize a single cost function incorporating multiple quantities (local SAR, global SAR, and per-channel power) averaged over the entire multislice scan subject to a strict constraint on excitation accuracy. Pulses for an 8-channel 7T head coil were designed with IMPULSE and compared with pulses designed using generic optimization algorithms and VOPs to assess the computation time and SAR performance benefits.

RESULTS

IMPULSE achieves lower SAR and shorter computation time compared with a VOP approach. Compared with the generic sequential quadratic programming algorithm, computation time is reduced by a factor of 5-6 by using IMPULSE. Using as many as 6 million local SAR terms, up to 120 slices can be designed jointly with IMPULSE within 45 s.

CONCLUSIONS

IMPULSE can handle significantly larger number of SAR matrices and slices than conventional optimization algorithms, enabling the use of uncompressed or partially compressed SAR matrices to design pulses for a multislice scan in a practical time frame.

Pros and cons of ultra-high-field MRI/MRS for human application

Magnetic resonance imaging and spectroscopic techniques are widely used in humans both for clinical diagnostic applications and in basic research areas such as cognitive neuroimaging. In recent years, new human MR systems have become available operating at static magnetic fields of 7 T or higher (≥300 MHz proton frequency). Imaging human-sized objects at such high frequencies presents several challenges including non-uniform radiofrequency fields, enhanced susceptibility artifacts, and higher radiofrequency energy deposition in the tissue. On the other side of the scale are gains in signal-to-noise or contrast-to-noise ratio that allow finer structures to be visualized and smaller physiological effects to be detected. This review presents an overview of some of the latest methodological developments in human ultra-high field MRI/MRS as well as associated clinical and scientific applications. Emphasis is given to techniques that particularly benefit from the changing physical characteristics at high magnetic fields, including susceptibility-weighted imaging and phase-contrast techniques, imaging with X-nuclei, MR spectroscopy, CEST imaging, as well as functional MRI. In addition, more general methodological developments such as parallel transmission and motion correction will be discussed that are required to leverage the full potential of higher magnetic fields, and an overview of relevant physiological considerations of human high magnetic field exposure is provided.

High-resolution whole-brain diffusion MRI at 7T using radiofrequency parallel transmission

PURPOSE

Investigating the utility of RF parallel transmission (pTx) for Human Connectome Project (HCP)-style whole-brain diffusion MRI (dMRI) data at 7 Tesla (7T).

METHODS

Healthy subjects were scanned in pTx and single-transmit (1Tx) modes. Multiband (MB), single-spoke pTx pulses were designed to image sagittal slices. HCP-style dMRI data (i.e., 1.05-mm resolutions, MB2, b-values = 1000/2000 s/mm² , 286 images and 40-min scan) and data with higher accelerations (MB3 and MB4) were acquired with pTx.

RESULTS

pTx significantly improved flip-angle detected signal uniformity across the brain, yielding ∼19% increase in temporal SNR (tSNR) averaged over the brain relative to 1Tx. This allowed significantly enhanced estimation of multiple fiber orientations (with ∼21% decrease in dispersion) in HCP-style 7T dMRI datasets. Additionally, pTx pulses achieved substantially lower power deposition, permitting higher accelerations, enabling collection of the same data in 2/3 and 1/2 the scan time or of more data in the same scan time.

CONCLUSION

pTx provides a solution to two major limitations for slice-accelerated high-resolution whole-brain dMRI at 7T; it improves flip-angle uniformity, and enables higher slice acceleration relative to current state-of-the-art. As such, pTx provides significant advantages for rapid acquisition of high-quality, high-resolution truly whole-brain dMRI data.

Design of spectral-spatial phase prewinding pulses and their use in small-tip fast recovery steady-state imaging

PURPOSE

Spectrally selective "prewinding" radiofrequency pulses can counteract B0 inhomogeneity in steady-state sequences, but can only prephase a limited range of off-resonance. We propose spectral-spatial small-tip angle prewinding pulses that increase the off-resonance bandwidth that can be successfully prephased by incorporating spatially tailored excitation patterns.

THEORY AND METHODS

We present a feasibility study to compare spectral and spectral-spatial prewinding pulses. These pulses add a prephasing term to the target magnetization pattern that aims to recover an assigned off-resonance bandwidth at the echo time. For spectral-spatial pulses, the design bandwidth is centered at the off-resonance frequency for each spatial location in a field map. We use these pulses in the small-tip fast recovery steady-state sequence, which is similar to balanced steady-state free precession. We investigate improvement of spectral-spatial pulses over spectral pulses using simulations and small-tip fast recovery scans of a gel phantom and human brain.

RESULTS

In simulation, spectral-spatial pulses yielded lower normalized root mean squared excitation error than spectral pulses. For both experiments, the spectral-spatial pulse images are also qualitatively better (more uniform, less signal loss) than the spectral pulse images.

CONCLUSION

Spectral-spatial prewinding pulses can prephase over a larger range of off-resonance than their purely spectral counterparts. Magn Reson Med 79:1377-1386, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Pulse sequences and parallel imaging for high spatiotemporal resolution MRI at ultra-high field

The SNR and CNR benefits of ultra-high field (UHF) have helped push the envelope of achievable spatial resolution in MRI. For applications based on susceptibility contrast where there is a large CNR gain, high quality sub-millimeter resolution imaging is now being routinely performed, particularly in fMRI and phase imaging/QSM. This has enabled the study of structure and function of very fine-scale structures in the brain. UHF has also helped push the spatial resolution of many other MRI applications as will be outlined in this review. However, this push in resolution comes at a cost of a large encoding burden leading to very lengthy scans. Developments in parallel imaging with controlled aliasing and the move away from 2D slice-by-slice imaging to much more SNR-efficient simultaneous multi-slice (SMS) and 3D acquisitions have helped address this issue. In particular, these developments have revolutionized the efficiency of UHF MRI to enable high spatiotemporal resolution imaging at an order of magnitude faster acquisition. In addition to describing the main approaches to these techniques, this review will also outline important key practical considerations in using these methods in practice. Furthermore, new RF pulse design to tackle the B1+ and SAR issues of UHF and the increased SAR and power requirement of SMS RF pulses will also be touched upon. Finally, an outlook into new developments of smart encoding in more dimensions, particularly through using better temporal/across-contrast encoding and reconstruction will be described. Just as controlled aliasing fully exploits spatial encoding in parallel imaging to provide large multiplicative gains in accelerations, the complimentary use of these new approaches in temporal and across-contrast encoding are expected to provide exciting opportunities for further large gains in efficiency to further push the spatiotemporal resolution of MRI.

Clinical Neuroimaging Using 7 T MRI: Challenges and Prospects

The aim of this article is to illustrate the principal challenges, from the medical and technical point of view, associated with the use of ultrahigh field (UHF) scanners in the clinical setting and to present available solutions to circumvent these limitations. We would like to show the differences between UHF scanners and those used routinely in clinical practice, the principal advantages, and disadvantages, the different UHFs that are ready be applied to routine clinical practice such as susceptibility-weighted imaging, fluid-attenuated inversion recovery, 3-dimensional time of flight, magnetization-prepared rapid acquisition gradient echo, magnetization-prepared 2 rapid acquisition gradient echo, and diffusion-weighted imaging, the technical principles of these sequences, and the particularities of advanced techniques such as diffusion tensor imaging, spectroscopy, and functional imaging at 7TMR. Finally, the main clinical applications in the field of the neuroradiology are discussed and the side effects are reported.

Advancing RF pulse design using an open-competition format: Report from the 2015 ISMRM challenge

PURPOSE

To advance the best solutions to two important RF pulse design problems with an open head-to-head competition.

METHODS

Two sub-challenges were formulated in which contestants competed to design the shortest simultaneous multislice (SMS) refocusing pulses and slice-selective parallel transmission (pTx) excitation pulses, subject to realistic hardware and safety constraints. Short refocusing pulses are needed for spin echo SMS imaging at high multiband factors, and short slice-selective pTx pulses are needed for multislice imaging in ultra-high field MRI. Each sub-challenge comprised two phases, in which the first phase posed problems with a low barrier of entry, and the second phase encouraged solutions that performed well in general. The Challenge ran from October 2015 to May 2016.

RESULTS

The pTx Challenge winners developed a spokes pulse design method that combined variable-rate selective excitation with an efficient method to enforce SAR constraints, which achieved 10.6 times shorter pulse durations than conventional approaches. The SMS Challenge winners developed a time-optimal control multiband pulse design algorithm that achieved 5.1 times shorter pulse durations than conventional approaches.

CONCLUSION

The Challenge led to rapid step improvements in solutions to significant problems in RF excitation for SMS imaging and ultra-high field MRI. Magn Reson Med 78:1352-1361, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Simultaneous use of linear and nonlinear gradients for B1+ inhomogeneity correction

The simultaneous use of linear spatial encoding magnetic fields (L-SEMs) and nonlinear spatial encoding magnetic fields (N-SEMs) in B₁⁺ inhomogeneity problems is formulated and demonstrated with both simulations and experiments. Independent excitation k-space variables for N-SEMs are formulated for the simultaneous use of L-SEMs and N-SEMs by assuming a small tip angle. The formulation shows that, when N-SEMs are considered as an independent excitation k-space variable, numerous different k-space trajectories and frequency weightings differing in dimension, length, and energy can be designed for a given target transverse magnetization distribution. The advantage of simultaneous use of L-SEMs and N-SEMs is demonstrated by B₁⁺ inhomogeneity correction with spoke excitation. To fully utilize the independent k-space formulations, global optimizations are performed for 1D, 2D RF power limited, and 2D RF power unlimited simulations and experiments. Three different cases are compared: L-SEMs alone, N-SEMs alone, and both used simultaneously. In all cases, the simultaneous use of L-SEMs and N-SEMs leads to a decreased standard deviation in the ROI compared with using only L-SEMs or N-SEMs. The simultaneous use of L-SEMs and N-SEMs results in better B₁⁺ inhomogeneity correction than using only L-SEMs or N-SEMs due to the increased number of degrees of freedom.

Homogeneous non-selective and slice-selective parallel-transmit excitations at 7 Tesla with universal pulses: A validation study on two commercial RF coils

Parallel transmission (pTx) technology, despite its great potential to mitigate the transmit field inhomogeneity problem in magnetic resonance imaging at ultra-high field (UHF), suffers from a cumbersome calibration procedure, thereby making the approach problematic for routine use. The purpose of this work is to demonstrate on two different 7T systems respectively equipped with 8-transmit-channel RF coils from two different suppliers (Rapid-Biomed and Nova Medical), the benefit of so-called universal pulses (UP), optimized to produce uniform excitations in the brain in a population of adults and making unnecessary the calibration procedures mentioned above. Non-selective and slice-selective UPs were designed to return homogeneous excitation profiles throughout the brain simultaneously on a group of ten subjects, which then were subsequently tested on ten additional volunteers in magnetization prepared rapid gradient echo (MPRAGE) and multi-slice gradient echo (2D GRE) protocols. The results were additionally compared experimentally with the standard non-pTx circularly-polarized (CP) mode, and in simulation with subject-specific tailored excitations. For both pulse types and both coils, the UP mode returned a better signal and contrast homogeneity than the CP mode. Retrospective analysis of the flip angle (FA) suggests that the FA deviation from the nominal FA on average over a healthy adult population does not exceed 11% with the calibration-free parallel-transmit pulses whereas it goes beyond 25% with the CP mode. As a result the universal pulses designed in this work confirm their relevance in 3D and 2D protocols with commercially available equipment. Plug-and-play pTx implementations henceforth become accessible to exploit with more flexibility the potential of UHF for brain imaging.

Motion-robust cardiac B1+ mapping at 3T using interleaved bloch-siegert shifts

PURPOSE

To develop and evaluate a robust motion-insensitive Bloch-Siegert shift based B1+ mapping method in the heart.

METHODS

Cardiac Bloch-Siegert B1+ mapping was performed with interleaved positive and negative off-resonance shifts and diastolic spoiled gradient echo imaging in 12 heartbeats. Numerical simulations were performed to study the impact of respiratory motion. The method was compared with three-dimensional (3D) actual flip angle imaging (AFI) and two-dimensional (2D) saturated double angle method (SDAM) in phantom scans. Cardiac B1+ maps of three different views were acquired in six healthy volunteers using Bloch-Siegert and SDAM during breath-hold and free breathing. In vivo maps were evaluated for inter-view consistency using the correlation coefficients of the B1+ profiles along the lines of intersection between the views.

RESULTS

For the Bloch-Siegert sequence, numerical simulations indicated high similarity between breath-hold and free breathing scans, and phantom results indicated low deviation from the 3D AFI reference (normalized root mean square error [NRMSE] = 2.0%). Increased deviation was observed with 2D SDAM (NRMSE = 5.0%) due to underestimation caused by imperfect excitation slice profiles. Breath-hold and free breathing Bloch-Siegert in vivo B1+ maps were visually comparable with no significant difference in the inter-view consistency (P > 0.36). SDAM showed strongly impaired B1+ map quality during free breathing. Inter-view consistency was significantly lower than with the Bloch-Siegert method (breath-hold: P = 0.014, free breathing: P < 0.0001).

CONCLUSION

The proposed interleaved Bloch-Siegert sequence enables cardiac B1+ mapping with improved inter-view consistency and high resilience to respiratory motion. Magn Reson Med 78:670-677, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Mitigation of B1+ inhomogeneity using spatially selective excitation with jointly designed quadratic spatial encoding magnetic fields and RF shimming

PURPOSE

The inhomogeneity of flip angle distribution is a major challenge impeding the application of high-field MRI. We report a method combining spatially selective excitation using generalized spatial encoding magnetic fields (SAGS) with radiofrequency (RF) shimming to achieve homogeneous excitation. This method can be an alternative approach to address the challenge of B1+ inhomogeneity using nonlinear gradients.

METHODS

We proposed a two-step algorithm that jointly optimizes the combination of nonlinear spatial encoding magnetic fields and the combination of multiple RF transmitter coils and then optimizes the locations, RF amplitudes, and phases of the spokes.

RESULTS

Our results show that jointly designed SAGS and RF shimming can provide a more homogeneous flip angle distribution than using SAGS or RF shimming alone. Compared with RF shimming alone, our approach can reduce the relative standard deviation of flip angle by 56% and 52% using phantom and human head data, respectively.

CONCLUSION

The jointly designed SAGS and RF shimming method can be used to achieve homogeneous flip angle distributions when fully parallel RF transmission is not available. Magn Reson Med 78:577-587, 2017. © 2016 International Society for Magnetic Resonance in Medicine.