0935

Cardiac bSSFP accelerated using minimum-TR multiband RF pulse design and GIRF-correction
Samy Abo Seada1, Anthony Price1, Joseph Hajnal1,2, and Shaihan Malik1

1Biomedical Engineering Department, School of Biomedical Engineering & Imaging Sciences, King's College London, London, United Kingdom, 2Centre for the Developing Brain, School of Biomedical Engineering & Imaging Sciences, King's College London, London, United Kingdom

Synopsis

Cardiac balanced SSFP using simultaneous multi-slice acquisition is limited by RF pulse duration and SAR. The interplay of these constraints means that the minimum pulse duration does not usually yield the minimum TR. We show that multiband RF pulses designed using VERSE to directly minimize TR can yield a reduced TR. The result is demonstrated in a multiband-2 sequence on a 3T clinical system where the TR is only 15% longer than for single-slice excitation and the time per slice is reduced by 44%.

Purpose

Cardiac bSSFP sequences are widely used for high SNR and pronounced blood-myocardial contrast. This imaging technique is often hampered by long uncomfortable breath-holds and banding artefacts due to $$$B_0$$$ inhomogeneity, both of which are reduced with shorter TR.

Minimum-duration RF pulses can in fact result in a longer TR due to SAR restrictions. Minimum-TR RF pulse design reformulates the pulse design problem to consider minimum-TR directly1. Furthermore, this technique makes use of VERSE2 to design time-variable gradients, which can lead to slice-profile distortions due to the finite temporal bandwidth of gradient systems3. This effect includes eddy currents, and can be accurately modelled by a Gradient Impulse Response Function (GIRF)4.

In this work we demonstrate how minimum-TR RF pulse design can be used in a multiband cardiac bSSFP acquisition5–7, with the additional use of GIRF-corrected VERSE RF pulses8 to overcome slice-profile distortions due imperfect gradient systems3.

Methods

The minimum-TR framework1 designs RF pulses that minimise a sequence TR using two constraints. Firstly, an SAR constraint is given by

$$ TR_{min} = E/SAR_{lim} \times \alpha$$

where $$$ E_{RF} $$$ is RF pulse energy, $$$SAR_{lim}$$$ is a SAR limit and $$$ \alpha $$$ is a conversion factor in $$$ W/kg/\mu T^2 $$$ determined from a SAR model. Secondly, time-constraints either on the total time to perform all sequence elements, or a duty-cycle limit, lead to

$$ TR_{min} = max(T_{enc} + T_{pulse}, T_{pulse}/\delta ) $$

where $$$T_{encoding}$$$ is the encoding time, $$$T_{pulse}$$$ is the pulse duration and $$$ \delta $$$ is an RF duty cycle limit. Figure 2b shows these two constraints for conventional and VERSE pulses for $$$SAR_{lim}$$$=20W/kg, encoding time $$$ T_{encoding} $$$= 1.786ms and $$$ \delta =50\%$$$. The minimum TR is found at the intersection of the curves defined by these equations. Use of VERSE to reshape the RF pulse changes the shape and alters the SAR constraint, resulting in a different optimum.

Time-variable selection gradients also lead to fidelity issues due to limited temporal bandwidth of gradient systems, which can be modelled and potentially corrected for by using a Gradient Impulse Response Function (GIRF)3,4,8. Our measured system GIRF3 roughly resembles a Lorentzian and so we approximated the GIRF using a simple Lorentzian model$$ H(t) = \frac{1}{\tau}e^{-t/\tau} $$where $$$\tau=42\mu s$$$ was determined empirically.

Time-variable selection gradients also require time-variable modulation for carrier for off-center slice excitation. For a slice-offset $$$ \Delta x $$$ the phase modulation is found by$$e^{i \theta(t)} = e^{i \gamma \Delta x \int_t^TG(s)ds}$$

For improved results, the gradient term was also found using the GIRF.

Experiments

Vendor RF pulses with time-bandwidth product 2.13 were scaled to flip angles 25 to 90 degrees. Time-optimal VERSE2 was applied to create VERSE singleband pulses. VERSE multiband (vMB) 2-3 pulses were created by applying a time-dilated multiband modulation function to the VERSE SB pulses3.Slice profiles were measured in a gadolinium-doped phantom using a 2D gradient-echo sequence, on a clinical Philips 3T Achieva system. On the same system, multiband 2 cardiac bSSFP images were acquired in a single healthy volunteer using blipped-SSFP5,9. MB2 reconstruction used a SENSE-based algorithm in ReconFrame (Gyrotools, GmbH, Zurich). To overcome a scanner software issue SAR limit was set to 18.6 W/kg for phantom work.

Results & Discussion

Figure 3a shows the minimised TR for different flip-angles. The TR can always be reduced by using time-variable selection gradients, for multiband factors 1-3. Figure 3b shows the same results, expressed as time spent per slice, which directly relates to scan time reduction. When accelerating a singleband acquisition with conventional multiband 2, an average of 36% of reduction in time per slice can be achieved. However this acceleration comes at an average increase of 29% in TR, making breath-holding more challenging and increasing vulnerability to artefacts. A VERSE MB2 implementation reduces the time-per-slice by 43% with only an average increase of 15% in TR, improving banding artefacts.

Figure 4 shows off-center slice-profile measurements using MB2 and vMB2 pulses, without and with GIRF correction. GIRF-corrected RF pulses result in reduced slice sidelobes and were used for in-vivo imaging.

Figure 5 shows in-vivo images of cardiac bSSFP multiband acquisition with and without the minimum TR framework. A minimum-TR RF pulse design approach results in the optimal reduction in TR, reducing banding artefacts and resulting in shorter breath-holds.

Conclusion

This work applies minimum-TR RF pulse design for multiband-accelerated cardiac bSSFP, using GIRF-corrected RF pulses. Key to this work is that minimum RF pulse duration does not imply minimum TR. Other rapid gradient-echo applications which could benefit from this approach include TOF angiography10 and quantitative imaging sequences11.

Acknowledgements

This work was supported by the Wellcome EPSRC Centre for Medical Engineering at Kings College London (WT 203148/Z/16/Z), MRC strategic grant MR/K006355/1, King’s College London & Imperial College London EPSRC Centre for Doctoral Training in Medical Imaging (EP/L015226/1) and by the National Institute for Health Research (NIHR) Biomedical Research Centre based at Guy’s and St Thomas’ NHS Foundation Trust and King’s College London. The views expressed are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health. SJM acknowledges support from the EPSRC (EP/L00531X/1). The student stipend for SAS is partially funded by Philips Healthcare.

References

1. Abo Seada S, Beqiri A, Price AN, Hajnal JV, Malik SJ. Minimum-TR pulse design for rapid gradient echo sequences. In: Proc. Intl. Soc. Mag. Reson. Med. 26. ; 2018. p. 5800. doi: 10.1002/nbm.3701.6.

2. Lee D, Lustig M, Grissom WA, Pauly JM. Time-optimal design for multidimensional and parallel transmit variable-rate selective excitation. Magn. Reson. Med. 2009;61:1471–1479. doi: 10.1002/mrm.21950.

3. Abo Seada S, Price AN, Schneider T, Hajnal JV., Malik SJ. Multiband RF pulse design for realistic gradient performance. Magn. Reson. Med. 2018:1–15. doi: 10.1002/mrm.27411.

4. Vannesjo SJ, Haeberlin M, Kasper L, Pavan M, Wilm BJ, Barmet C, Pruessmann KP. Gradient system characterization by impulse response measurements with a dynamic field camera. Magn. Reson. Med. 2013;69:583–593. doi: 10.1002/mrm.24263.

5. Price AN, Cordero-Grande L, Malik SJ, Hajnal JV. Accelerated Cine Imaging of the Heart using Blipped Multiband SSFP. Int. Soc. Magn. Reson. Med. 2017:0631. doi: 10.1002/mrm.25897.

6. Wang Y, Shao X, Martin T, Moeller S, Yacoub E, Wang DJJ. Phase-cycled simultaneous multislice balanced SSFP imaging with CAIPIRINHA for efficient banding reduction. Magn. Reson. Med. 2016;76:1764–1774. doi: 10.1002/mrm.26076.

7. Landes V, Jao T, Nayak K. Practical implementation of SMS bSSFP in the Heart. In: Proc. Intl. Soc. Mag. Reson. Med. 26. ; 2018. pp. 18–20.

8. Abo Seada S, Hajnal JV, Malik SJ. A simple optimisation approach to making time efficient VERSE-multiband pulses feasible on non-ideal gradients. Proc. Intl. Soc. Mag. Reson. Med. 25 2017:5049.

9. Setsompop K, Gagoski BA, Polimeni JR, Witzel T, Wedeen VJ, Wald LL. Blipped-controlled aliasing in parallel imaging for simultaneous multislice echo planar imaging with reduced g-factor penalty. Magn. Reson. Med. 2012;67:1210–1224. doi: 10.1002/mrm.23097.

10. Schmitter S, Bock M, Johst S, Auerbach EJ, Ugurbil K, Van De Moortele PF. Contrast enhancement in TOF cerebral angiography at 7 T using saturation and MT pulses under SAR constraints: Impact of verse and sparse pulses. Magn. Reson. Med. 2012;68:188–197. doi: 10.1002/mrm.23226.

11. Hardy PA, Recht MP, Piraino D, Thomasson D. Optimization of a Dual Echo in the Steady State (DESS) free-precession sequence for imaging cartilage. J. Magn. Reson. Imaging 1996;6:329–335. doi: 10.1002/jmri.1880060212.

Figures

Top: A simplified sequence diagram with a conventional MB2 pulse results in a TR of 4.06ms.

Bottom: For the same design, a shorter TR of 3.56ms can be achieved by using minimum-TR RF pulse design, designed here using time-optimal VERSE. The minimum-TR optimization space, with the SAR and encoding constraints, is shown in Figure 2b.


a) Three example MB2 and (VERSE) vMB2 scaled to different peak amplitudes. These pulses have the same design but result in different durations. b) Minimum TR RF pulse design is achieved by finding the RF pulse duration which minimizes TR based on a SAR constraint and an image encoding constraint (red).

a) Minimum TR achieved for conventional constant gradient and VERSE pulses, for different MB-factors. Using VERSE for minimum-TR RF pulse design always reduces TR, however larger relative gain is achieved at higher flip-angles. b) Performance expressed as time per slice (120 phase-encode steps, 20 cardiac phases) is more relevant for reduced acquisition time, and thus breath-hold periods.

Top: MB and vMB RF pulses as well as their respective gradients and predicted gradient distortion using the GIRF. It is possible to correct the RF pulse using a GIRF model. Bottom: Slice profile measurements using a cylindrical phantom, in a 2D gradient-echo acquisition (TR=50ms, TE = 8.5ms, 0.55x0.575mm in-plane resolution). Low temporal bandwidth of gradient systems, as modelled by the GIRF, results in increased side-lobes. These can be corrected for by using GIRF-corrected RF pulses. Slices are shifted away from isocenter to view additional errors from non-GIRF corrected slice-shifting. GIRF-correction results in a small increase in RF energy.

Cardiac balanced SSFP acquisition using minimum TR RF pulss and multiband 2 acquisition at 3T with GIRF-corrected VERSE RF pulses. (FA=40o, slice-thickness = 7mm, 1.6x2mm resolution, no in-plane undersampling, 32channel receiver coil). From left to right, the first column shows singleband reference data. Second column shows a conventional MB2 acquisition with minimum TR. Third column shows vMB2 (VERSE) data, without using the minimum-TR framework, which leads to a 0.9ms increase in TR. The minimum-TR framework on the right is shown to reduce this approach to 2.9ms. Time per slice is calculated using 120 phase-encodes and 20 cardiac phases.

Proc. Intl. Soc. Mag. Reson. Med. 27 (2019)
0935