Sequences and Simulations
Martin John Graves1

1Radiology, Cambridge University Hospitals, Cambridge, United Kingdom

Synopsis

This presentation will provide an overview of the main gradient echo based (gradient spoiled, RF spoiled and balanced steady state free precession) and conventional/fast spin echo based pulse sequences and will illustrate some methods by which their behaviour can be simulated

Introduction

The behaviour of the magnetisation in a pulse sequence is described by the Bloch equations (1) and a number of methods have been utilised to simulate the behaviour of spins in an imaging sequence such as the matrix formalism (2), the use of extended phase graphs (EPG) (3, 4), pictorial representations (5) and analytical descriptions (6, 7). The examples given in this lecture are based upon code from Brian Hargreaves' Bloch equation simulations (http://wwwmrsrl.stanford.edu/~brian/bloch/), Matthias Weigel's EPG simulations (http://epg.matthias-weigel.net/) and a Matlab implementation of the C program given in Klaus Scheffler's article (5).

Gradient Echo (GE)

A gradient echo sequence comprises a train of equally spaced radiofrequency (RF) excitation pulses and the necessary gradients to perform slice selection, phase encoding and frequency encoding. The echo signal is formed by the application of a negative dephasing gradient followed by the frequency encoding gradient. The minimum echo time (TE) is dependent upon the RF pulse duration, receiver bandwidth and available gradient performance. Usually the time between excitation pulses (TR) is short and the transverse magnetisation has not had time to fully dephase. Since the phase encoding gradient is the only gradient that changes between TRs the transverse magnetisation will have a position-dependent distribution of resonant offset phase angles across the slice, which will affect the signal created by the subsequent RF pulse. As a result, bands of varying signal intensity, known as "FLASH-bands" may appear degrading the image quality (Figure 1a). The three main gradient echo sequences differ in how they deal with this residual transverse magnetisation: a) the gradient spoiled sequence averages the resonant offset phase distribution across a voxel (Figure 1b), b) the radiofrequency (RF) spoiled sequence eliminates the transverse magnetisation (Figure 1c) whilst c) the balanced steady state free precession (bSSFP) sequence rephases the phase distribution. Example pulse sequence diagrams and images are shown in Figure 2.

Gradient Spoiling

The gradient spoiled sequence has additional gradients applied after data acquisition on the slice select and frequency encoding axes (Figure 2a) that averages the resonant offset phase distribution across a voxel. The phase encoding gradient is "rewound" so that the dephasing gradients are constant across TRs. Figure 2b shows an axial head image obtained using a gradient spoiled sequence.

RF Spoiling

The gradients are the same for both the gradient and RF spoiled sequences (Figure 2a), but the phase of the RF excitation pulse is incremented, usually quadratically, on each excitation.

$$ ϕ_j = \frac{1}{2} ϕ_0 (j^2+j+1),j=0,1,2,… $$

There are certain “magic numbers”, e.g. $$$ϕ_0=117^0$$$, for the for the phase increment that can produce excellent spoiling across a range of T1, T2 and excitation flip angles (8) and is equal to the steady-state transverse magnetisation for a perfectly spoiled sequence given by the Ernst equation

$$ M_{ss}=M_0 (1-exp(-TR/T_1 ))/(1-cos(α) exp(-TR/T_1 ) ) \cdot sin(α) $$

Balanced Steady State Free Precession(bSSFP)

The bSSFP sequence (Figure 2d) has all gradients balanced, i.e. there is zero net gradient area on each axis over a TR period. Any resonant offset phase errors will therefore not be averaged. This means that the signal response is dependent upon the phase offset. Figure 3 shows a simulation of the signal behaviour as a function of resonant offset phase angle for two excitation schemes. In Fig 3a and 3b there is no phase offset applied between subsequent RF pulses. In Figure 3c and 3d the phase of the RF excitation is alternated between 0 and 180°. In this example a large left-right shim offset of 14.7Hz/cm was applied. The resonant offset is given by $$$θ=2π\cdot Δf\cdot TR$$$ and the signal nulls or bands appear at $$$θ=±π$$$ , or when $$$Δf=±\frac{1}{2\cdot TR} $$$ . In this case the TR was 5ms resulting in bands appearing at $$$±$$$100Hz . Figure 2e shows the image obtained using a bSSFP sequence. Note an off resonance band can be seen anteriorly. The steady state signal for a bSSFP sequence, assuming that $$$TR≪T_1,T_2$$$, is given by

$$M_{ss}=M_0 \frac {sin(α)} {(1+cos(α)+(1-cos(α) ) \cdot \frac{T_1}{T_2}}$$

The optimal flip angle is given by

$$cos(α)=\frac{\frac{T_1}{T_2} -1}{\frac{T_1}{T_2} +1}$$

which results in a steady state signal response of

$$M_{ss}=\frac{1}{2} M_0 \sqrt{\frac{T_2}{T_1}}$$

Hence tissues with a high $$$\frac{T_2}{T_1}$$$ ratio, e.g. fluids, will exhibit a high signal.

Spin Echo (SE) and Fast Spin Echo (FSE)

Spin Echo

Spin echo (SE) is a basic pulse sequence comprising typically of a 90° and a 180° RF pulse. The initial 90° excitation pulse tips the longitudinal magnetisation into the transverse plane where it naturally dephases due to T2* relaxation processes. After a time τ the 180° refocusing pulse is applied which rotates the transverse magnetisation 180° about the axis along which the pulse is applied. This rotation effectively inverts the phase of the magnetisation so that it now naturally rephases, coming fully back into phase, and forming what is known as a spin echo, after a further time period τ. The time between the initial 90° excitation pulse and the centre of the echo is known as the echo time (TE) and is equal to 2τ. The effect of the refocusing pulse is to eliminate the effect of any static magnetic field non-uniformities from the T2* signal therefore leaving any signal differences due to irreversible T2 interactions. The magnetisation will then recover along the z-direction due to T1 relaxation until the application of the next 90° pulse. The time between repeat applications of the 90° pulse is the repetition time (TR). The signal equation is:

$$M=M_0 (1-exp(-TR/T_1 )\cdot exp(-TE/T_2)$$

Fast Spin Echo

Proton density or T2-weighted spin echo sequences are particularly slow since they require a long TR in order for there to be appreciable longitudinal magnetisation recovery before the next excitation pulse. Although the acquisition of other slices can be interleaved within the TR period the sequence is still not particularly efficient. Fast, or turbo, spin echo (FSE/TSE) was developed as a method for reducing the overall acquisition time (9). The basic principle is to use a Carr-Purcell-Meiboom-Gill (CPMG) multi-echo spin echo acquisition, where the transverse magnetisation is repeatedly refocused by the application of a train of refocusing pulses forming multiple spin echoes each at correspondingly increasing TEs. If each echo is individually phase encoded, then they can be used in the same image k-space resulting in a shorter acquisition time. It is necessary to arrange the amplitudes of the phase encoding gradients so that the echo nearest the desired TE (known as the effective echo time (TEeff), has the lowest amplitude phase encoding gradient, i.e. that echoes’ data will be nearest the centre of k-space. The interval between refocusing pulses is known as the Echo Spacing (ESP) and the number of refocusing pulses per TR is known as the Echo Train Length (ETL). Since multiple echoes are acquired per TR there is also a filtering effect due to the T2 decay of the echo amplitudes over the ETL. The T2 decay results in a broad point spread function which has the effect of blurring the images. A limitation of FSE is that the large power deposition associated with long echo trains. A common solution to this problem is to use reduced flip angle refocusing pulses that substantially reduce the power with only a minor impact on the signal-to-noise ratio (SNR). The use of refocusing pulses $$$\neq$$$180° do, however, give rise to a number of coherence pathways which lead to the creation of a so-called pseudo-steady-state (PSS) (10). These coherence pathways can be most easily studied using the EPG method in which the magnetisation is described as substates which are defined as groups of spins with different Mx, My and Mz (3). Figure 4 shows a hand-drawn EPG for an initial excitation pulse followed by three refocusing pulses . The essential feature of the phase graph concept is that after an RF pulse has been applied the magnetisation behaves as a composition of three parts a) dephasing transverse magnetisation, b) rephasing transverse magnetisation and c) longitudinal magnetisation (4). The dephased magnetisation after the excitation pulse is labelled F1. The first refocusing pulse will invert a component of the magnetisation, creating the F-1 substate and will also produce z-magnetisation represented by Z1 and its complex conjugate Z-1. Since the phase information does not evolve with time these longitudinal substates are represented as horizontal (red) lines. The dephasing of F1 gives rise to F2, whilst rephasing of F-1 causes the formation of a echo (F0) shown by the pink circle. Further time evolution creates higher order substates which contribute to the amplitude of each echo. The evolutions of transverse substates from larger number of pulses are usually displayed as shown in Figure 5.

Summary

The methods outlined in this article, and presented in depth in the references, can be used to simulate a range of MRI pulse sequences. Example applications for these methods include designing magnetisation preparation schemes to rapidly stabilise the signals at the start of a spoiled gradient echo (11), bSSFP (12) or CPMG FSE (13) acquisition or developing variable flip angle schemes to manage RF power deposition and blurring (14-16).

Acknowledgements

Addenbrooke’s Charitable Trust and the NIHR comprehensive Biomedical Research Centre award to Cambridge University Hospitals NHS Foundation Trust in partnership with the University of Cambridge.

References

1. Bloch F. Nuclear induction. Phys Rev 1946; 70:460-474.

2. Jaynes E. Matrix treatment of nuclear induction. Phys Rev 1955; 98:1099-1105.

3. Hennig J. Multiecho imaging sequences with low refocusing flip angles. J Magn Reson 1988; 78:397-407.

4. Weigel M. Extended phase graphs: dephasing, RF pulses, and echoes - pure and simple. J Magn Reson Imaging 2015; 41:266-295.

5. Scheffler K. A pictorial description of steady-states in rapid magnetic resonance imaging. Concepts in Magn. Res. 1999; 11:291-304.

6. Sobol WT, Gauntt DM. On the stationary states in gradient echo imaging. J Magn Reson Imaging 1996; 6:384-398.

7. Hanicke W, Vogel HU. An analytical solution for the SSFP signal in MRI. Magn Reson Med 2003; 49:771-775.

8. Zur Y, Wood ML, Neuringer LJ. Spoiling of transverse magnetization in steady-state sequences. Magn Reson Med 1991; 21:251-263.

9. Hennig J, Nauerth A, Friedburg H. RARE imaging: a fast imaging method for clinical MR. Magn Reson Med 1986; 3:823-833.

10. Alsop DC. The sensitivity of low flip angle RARE imaging. Magn Reson Med 1997; 37:176-184. 11. Busse RF, Riederer SJ. Steady-state preparation for spoiled gradient echo imaging. Magn Reson Med 2001; 45:653-661.

12. Le Roux P. Simplified model and stabilization of SSFP sequences. J Magn Reson 2003; 163:23-37.

13. Hennig J, Scheffler K. Easy improvement of signal-to-noise in RARE-sequences with low refocusing flip angles. Rapid acquisition with relaxation enhancement. Magn Reson Med 2000; 44:983-985.

14. Hennig J, Weigel M, Scheffler K. Calculation of flip angles for echo trains with predefined amplitudes with the extended phase graph (EPG)-algorithm: principles and applications to hyperecho and TRAPS sequences. Magn Reson Med 2004; 51:68-80.

15. Worters PW, Hargreaves BA. Balanced SSFP transient imaging using variable flip angles for a predefined signal profile. Magn Reson Med 2010; 64:1404-1412.

16. Busse RF. Reduced RF power without blurring: correcting for modulation of refocusing flip angle in FSE sequences. Magn Reson Med 2004; 51:1031-1037.


Figures

Gradient echo sequence (TE/TR = 10ms/20ms). (a) Without any spoiling demonstrating FLASH-band artefact. (b) Gradient spoiling. (c) RF spoiling.

(a) Gradient echo sequences with the spoiler gradients shown in grey. The red dotted lines demonstrate the readout period. (b) Axial image using a gradient spoiled sequence (TE/TR/Flip = 11.4ms/40ms/30°). (c) Axial image using an RF spoiled sequence (TE/TR/flip = 11.4ms/23ms/30°) with a quadratic phase increment of 117°. (d) bSSFP sequence (TE/TR/Flip =1.3/4.7/35°). (e) Axial image using a bSSFP sequence.

Balanced steady state free precession. (a) Signal response with no excitation phase alternation. (b) Phantom image obtained with a large left-right shim offset. (c) Signal response using a 180° phase alternation scheme. (d) Phantom image obtained with the same shim offset using the phase alternation scheme.

Fast spin echo EPG for the first three refocusing pulses. The refocusing pulses are represented by the grey parallel lines. The red horizontal lines represent substates of z-magnetisation, whilst the green lines represent time evolution of transverse magnetisation. The blue circles represent echoes.

Extended phase graph (EPG) for a 70 pulse CPMG FSE sequence using a refocusing flip angle of 90°. Relaxation effects are ignored in this example. The colour scale is proportional to the magnitude of Fn .




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