TR-FOCI pulse optimisations for SS-SI VASO

CBV-fMRI with VASO is highly dependent on a good inversion contrast. It gives it its CBV sensitivity and is also responsible for most of the VASO specific pitfalls (e.g. inflow, CSF etc. ). And thus, it should be optimized as much as possible.

In this blog post, I want to describe the most important features of a reliable inversion pulse for the application of VASO at 7T with a head transmit coil.

Pulse features

One of the most SAR-efficient adiabatic inversion pulse for 7T the TR-FOCI pulse class developed in Nottingham. While the original TR-FOCI pulse is optimized to achieve a sharp slab profile, VASO is not so much limited by sharpness. Instead, inversion in VASO it rather limited by SAR-efficiency and inversion range. Thus, it can be beneficial to reduce the bandwidth of the TR-FOCI pulse for VASO applications.

Screenshot 2018-11-29 16.47.51.png
Changes in VASO that save a bit of SAR. The lower bandwidth makes the inversion-profile less sharp. However, it allows smaller amplitudes (reduces SAR) at the same adiabaticity level.

Required B1+field

In the original TR-FOCI paper, Hurley et al. state that a peak amplitude of 7μT are sufficient to achieve a complete inversion. In my own validation experiments (figure below), I found that this is only true in phantoms with unrealistically small off-resonances. In-vivo, however, I find that in actual experiments at least 7.7 μT are necessary. When additional off resonances are present, like in feeding arteries at the base of the brain, I would recommend to add an extra safety margin. This is particularly necessary because we use lower bandwidths (4.7 kHz) than the original TR-FOCI pulse (6.4 kHz)

Screenshot 2018-11-30 11.35.46.png
Slab profile of the tr-FOCI pulse as a function of the magnetic field strength B1 in oil. It is apparent that the inversion efficiency reaches its maximal level for RF field strengths above 7.7μT. The shown percent values refer to the reference value of 25.91μT.

B1-independent way of reducing the inversion efficiency to minimize inflow effects

In SS-SI VASO, inflow effects can be challenging. As such, in brain areas with very short arterial arrival times (e.g. Auditory and or insular cortex) or in combination with global respiration tasks that reduce the arterial time (e.g. breath holding and hypercapnia). In those cases, it can happen that in the time period between the inversion and the excitation, fresh (non-inverted) blood reaches the micro-vessels of the imaging slab. This fresh blood results in a FAIR-like CBF contamination that counteracts the VASO contrast and should be avoided.

Connection of reduced inversion efficiency and inflow effect. The application of a conventional full inversion can result in long blood-nulling times (TIs) and thus introduces inflow effects. Reducing the inversion efficiency, also reduces the blood-nulling time to be shorter than the arterial arrival time. The content of this figure is explained in the following talk section on Youtube

Technical implementation of the B1-independent reduction in inversion efficiency

In order to achieve an adjustable inversion efficiency in a B1-independent manner, the concept of RF phase modulation during the RF pulse can be applied. Such phase modulations are known from plane rotation pulses for flip-angle variations [Norris and Haase, 1989] and were transferred to conventional (non-planar) adiabatic pulses. The reduction in inversion efficiency can be achieved by introducing a phase skip of the RF field, B1, during the inversion exactly at the time when the frequency of the adiabatic pulse is on resonance. The corresponding shape functions of the amplitude and phase of the RF pulse is shown in the figure below.

Amplitude and phase shape functions of the partial inversion RF pulse used to avoid inflow effects in SS-SI VASO: The partial inversion pulse is based on a TR-FOCI pulse with the pulse duration D. The partial inversion efficiency is adjusted by the phase skip (black arrow). Here, an example of 60◦ is depicted.

A more intuitive illustration of the working principle of the new partial inversion pulse is illustrated blow. In conventional adiabatic inversion pulses, the magnetization precesses around the effective magnetic field at the surface of a narrow ‘cone of precession’. During the frequency sweep of adiabatic RF pulses this cone usually follows the effective magnetic field, while it changes its direction. When a phase skip of the RF pulse is applied, this ‘cone of precession’ opens up. After inversion, the ‘cone of precession’ points along anti-parallel to the z-direction. The B1-independent reduction in inversion efficiency is then given by the transverse component of the magnetization, which depends directly on the phase skip applied. The reduction of the inversion efficiency is given by the cosine of the opening angle.

In conventional adiabatic inversion pulses of a pulse duration D, the magnetization follows the effective magnetic field in a narrow ‘cone of precession’ as the RF frequency offset is inverted (upper row). For partial inversion, the first half of the pulse is played out in the conventional manner. Half way throughout the pulse, when the RF-frequency is exactly the Larmor frequency, a phase skip is introduced that opens up the ‘cone of precession’. The second half of the pulse is again played out in the conventional manner. After applying a spoiler gradient, the magnetization, pointing opposite to the external magnetic field, is reduced with respect to the conventional inversion pulse. This reduction can be qualitatively adjusted by setting the angel of the phase skip applied. 

It is important to note, however, that the reduction of the TI to minimize inflow effect, also comes along with a reduced GM signal and reduced VASO sensitivity (See PhD Thesis Fig. 4.24).

Stability across participants

Personally, I found the TR-FOCI pulse very reliable across people. In fact, I find the inter-subject variability smaller than the error of the transmit adjust calibration scan. And thus, with the 32ch Nova coil, I fix the transmit voltage to 220V independent of the persons head size. The B1+ field maps blow confirm that the threshold of 7.7μT is exceeded for all people well below the circle of Willis.

B1 maps of 20 people. The volume with field strengths above 30 μT varies strongly between the subjects. Although the region in the neck, where required threshold of 7.7 μT is exceeded, does not vary strongly between the subjects.
B1 in μT as a function of position along the z-axis for all subjects. The threshold of 7 μT is required for full inversion. The abscissa is given in 6mm voxels starting from the upper bound of the cranium. All subjects were oriented in a way that the upper bound of the cranium is placed at the second voxel.

Other comments

  • It is to note, that the first VASO versions actually had a gradient implemented to make the SS-SI-VASO real slab-selective. In practice, however, it turned out that the slab thickness should be thicker than 14cm. Thus, the slab-thickness it actually larger than the range of the transmit coil coverage. As a result, the gradient is switched off for all scans with head-transmit coil only.
  • The inversion pulse can be implemented in multiple ways in the sequence.
    • In most of the SIEMENS VASO sequences, it is implemented as an ARB_rf. This means the pulse shape is calculated in real time based on the equations given from Hurley.
    • In the SIEMENS pulse-tool with this pulse tool file.
    • In Bruker with this pulse file.
    • This pulse file has also been applied on Phillips scanner by Peiying Liu.
  • The B1/B0 and SAR constraints are amplified at even higher field strengths of 9.4T. There, additional precautions need to be considered. This is described here:

Further Reading

  • Theoretical background of TR-FOCI pulses in Chapter 2.3.3 page 22 of the Masters thesis here.
  • Regime of reliable parameters of the TR-FOCI pulse in Chapter 6.1-6.3 (page 60) of the Masters thesis here.
  • Strategies to avoid inflow with reduced inversion efficiencies in Chapter 4.3 (page 123) of the PhD thesis here.


Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s