Universal pulses: a new concept for calibration-free parallel transmission
Vincent Gras1, Alexandre Vignaud1, Alexis Amadon1, Denis Le Bihan1, and Nicolas Boulant1

1Neurospin, CEA/DSV/I2BM, Gif-sur-Yvette, France

Synopsis

At ultra-high field, a drawback of parallel transmission to mitigate the RF inhomogeneity problem is the necessity to measure subject-specific field maps in order to return optimized RF pulses, thereby decreasing the time available for clinically-relevant scans. In this work, we investigate numerically and experimentally at 7T the design of "universal" kT-points pulses, which does not require the aforementioned calibration step but yet considerably improves excitation homogeneity compared to the standard circularly-polarized and RF shim modes. Such approach can simplify considerably the workflow of parallel transmission and render the potential of ultra-high field scanners more accessible to anyone in routine.

Introduction

Parallel transmission (pTx) enables excellent control of the magnetization to mitigate the transmit field inhomogeneity at ultra-high field (UHF) by means of advanced RF pulse design techniques1-4. However, to date the implementation of these methods necessitates the knowledge of the transmit RF field (B1) and the static field (ΔB0) distributions whereas the pulse design combined with the measurement and processing of those maps can easily cumulate 15 minutes, which decreases by the same amount the time available for acquiring clinically relevant data. Several studies on the other hand have suggested that the B1 distribution exhibits a reproducible pattern across different subjects5-7. In this work, we investigate numerically and experimentally for brain imaging at 7T on an 8 channel pTx system, the design of kT-points pulses4 qualified as universal, i.e. whose design does not involve the subject-specific field distributions, but which yet considerably improves performance compared to the standard circularly-polarized (CP) and RF shim modes.

Methods

A so-called field database, composed of representative ΔB0 and B1 maps, was acquired experimentally on $$$N_s=6$$$ subjects. The universal pulses, characterized by the kT-points trajectory $$$k$$$ and the RF pulse coefficients $$$b$$$, were then calculated by minimizing the objective function:

$$f(b,k) = \text{max}_{1\le i\le N_s} \| \mathcal{A}_i(b,k) - \alpha_t \|_2,$$

where $$$\alpha_t$$$ denotes the target flip angle (FA) and $$$\mathcal{A}_i(b,k)$$$ denotes the FA map generated by the pulse $$$(b, k)$$$ on subject $$$i$$$. The simultaneous optimization of the RF coefficients and the blipped k-space trajectory was performed under explicit power and SAR constraints to satisfy hardware and patient safety limits8,9, for one excitation (9°) and one inversion (180°) pulse. SAR matrices10 were obtained from electric field simulations with HFSS (Ansys, Canonsburg, PA, USA) on a generic head model and for the home-made pTx coil under study, and were then compressed with the virtual observation point technique11. An overall SAR safety margin of 2.7 was applied to account for modelling errors, anatomic variability and uncertainties in the SAR monitoring hardware12-14. The universal pulses, designed with the field database, were then blindly applied on 6 additional subjects by acquiring MPRAGE images (TR/TI=1.1/2.6 s, TE=3 ms, nominal FA=9, resolution= 1×1×1 mm3, TA=9 min), without any field map measurement or further pulse optimizations. The excitation and inversion pulses were designed with 5/7 kT-points and 0.7/4 ms respectively. For comparison, two additional MPRAGE acquisitions were performed, in one case with an adiabatic inversion and rectangular excitation pulses driven in CP mode, and in the other case with similar optimized kT-points pulses as for the universal pulses but using this time the subject-specific ΔB0 and B1 maps. The RF shim mode, with an adiabatic parametrization for the inversion pulse, was also investigated numerically. Results were compared quantitatively by computing the FA normalized root mean square errors (NRMSEs) and qualitatively by visual inspection of the MPRAGE images.

Results

The measured B1 distributions of each transmit channel are shown in Fig. 1 for subjects #1 to #12. The magnitude of the complex correlation coefficient between the subject-specific B1 distributions is 0.95 ± 0.02 for each transmit channel, suggesting that B1 varies only slightly from subject to subject despite the variability of head shape and position. The FA NRMSE values are reported in Fig. 2 for the 6 subjects constituting the field database and on whom the universal pulses were designed, and for the 6 additional subjects on whom the universal pulses were experimentally tested. Both the excitation and inversion universal pulses significantly outperform the CP and RF shim modes. Not surprisingly, the subject-based tailored pulses perform better than the universal pulses, the gain however being more marginal. The MPRAGE images obtained experimentally on two subjects are provided in Fig. 3. With the coil driven in CP mode, two problematic regions can be identified which are completely recovered with the use of the universal pulses. There is moreover no obvious degradation in image quality compared to the MPRAGE image acquired with the subject-based tailored pulses, despite the slightly poorer NRMSE. The same behavior was observed on all tested subjects #7-12. In that sense, the test of the universal pulses never failed.

Conclusion

A method to mitigate the RF field inhomogeneity problem at UHF, without subject-specific B1 and ΔB0 map measurements and online pulse design, was presented. The proof of concept was validated at 7 T with 3D brain imaging. Preliminary results suggest that the idea could be successfully applied in 2D. The results suggest that the use of parallel transmission thereby could become completely transparent to the user, hence making this technology more accessible for routine use.

Acknowledgements

The authors wish to thank Dr. Lucie Hertz-Pannier, Dr. Pierre Brugières, Dr. Franck Mauconduit and Dr. Guillaume Ferrand for valuable discussions. The research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Program (FP7/2013-2018) / ERCGrant Agreement n. 309674.

References

[1] Grissom et al. MRM 2006;56:620–629. [2] Saekho et al. MRM 2006;55:719–724. [3] Setsompop et al. MRM 2008;59:908–915. [4] Cloos et al. Magn Reson Med 2012;67:72–80. [5] Nistler et al. ISMRM, Berlin, Germany, 2007. p. 1063. [6] Moore et al. ISMRM, Salt Lake City, Utah, USA, 2013. p. 2593. [7] Kelley. ISMRM, Toronto, Canada, 2015. p. 2385. [8] Hoyos-Idrobo et al. IEEE TMI 33:739–748. [9] Gras et al. Accepted in JMR. [10] Greasslin et al. MRM 2012;67:72-80. [11] Eichfelder and Gebhardt. MRM 2011;66: 1468–1476. [12] Makris et al. Medical & Biological Engineering & Computing 2008;46:1239-1251. [13] Ferrand et al. ISMRM 2011; Abstract 3633. [14] Gumbrecht. Ph.D. thesis, Friedrich-Alexander Universität Erlangen-Nurnberg, Erlangen 2013.

Figures

Figure 1. Axial B1 maps in amplitude at isocenter. Images are sorted according to subject (#1-12 from top to bottom) and channel indexes (#1-8 from left to right).

Figure 2. FA NRMSEs for (a) the excitation and (b) the inversion pulses and for 4 pulse designs. The universal pulses were designed using only the fields acquired on subjects #1-6.

Figure 3. MPRAGE images obtained (from left to right) in CP mode, with the universal and with the subject-based tailored kT-point pulses on two subjects.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
1018