Validating subject-specific RF and thermal simulations in the calf muscle using MR-based temperature measurements

PURPOSE

Ongoing discussions occur to translate the safety restrictions on MR scanners from specific absorption rate (SAR) to thermal dose. Therefore, this research focuses on the accuracy of thermal simulations in human subjects during an MR exam, which is fundamental information in that debate.

METHODS

Radiofrequency (RF) heating experiments were performed on the calves of 13 healthy subjects using a dedicated transmit-receive coil while monitoring the temperature with proton resonance frequency shift (PRFS) thermometry. Subject-specific models and one generic model were used for electromagnetic and thermal simulations using Pennes' bioheat equation, with the blood equilibration constant equaling zero. The simulations were subsequently compared with the experimental results.

RESULTS

The mean B1+ equaled 15 µT in the center slice of all volunteers, and 95% of the voxels had errors smaller than 2.8 µT between the simulation and measurement. The intersubject variation in RF power to achieve the required B1+ was 11%. The resulting intersubject variation in median temperature rise was 14%. Thermal simulations underestimated the median temperature increase on average, with 34% in subject-specific models and 28% in the generic model.

CONCLUSIONS

Although thermal measures are directly coupled to tissue damage and therefore suitable for RF safety assessment, insecurities in the applied thermal modeling limit their estimation accuracy. Magn Reson Med 77:1691-1700, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

7 T renal MRI: challenges and promises

The progression to 7 Tesla (7 T) magnetic resonance imaging (MRI) yields promises of substantial increase in signal-to-noise (SNR) ratio. This increase can be traded off to increase image spatial resolution or to decrease acquisition time. However, renal 7 T MRI remains challenging due to inhomogeneity of the radiofrequency field and due to specific absorption rate (SAR) constraints. A number of studies has been published in the field of renal 7 T imaging. While the focus initially was on anatomic imaging and renal MR angiography, later studies have explored renal functional imaging. Although anatomic imaging remains somewhat limited by inhomogeneous excitation and SAR constraints, functional imaging results are promising. The increased SNR at 7 T has been particularly advantageous for blood oxygen level-dependent and arterial spin labelling MRI, as well as sodium MR imaging, thanks to changes in field-strength-dependent magnetic properties. Here, we provide an overview of the currently available literature on renal 7 T MRI. In addition, we provide a brief overview of challenges and opportunities in renal 7 T MR imaging.

Enhancement of Magnetic Resonance Imaging with Metasurfaces

It is revealed that the unique properties of ultrathin metasurface resonators can improve magnetic resonance imaging dramatically. A metasurface formed when an array of metallic wires is placed inside a scanner under the studied object and a substantial enhancement of the radio-frequency magnetic field is achieved by means of subwavelength manipulation with the metasurface, also allowing improved image resolution.

A 16-channel combined loop-dipole transceiver array for 7 Tesla body MRI

PURPOSE

To develop a 16-channel transceive body imaging array at 7.0 T with improved transmit, receive, and specific absorption rate (SAR) performance by combining both loop and dipole elements and using their respective and complementary near and far field characteristics.

METHODS

A 16-channel radiofrequency (RF) coil array consisting of eight loop-dipole blocks (16LD) was designed and constructed. Transmit and receive performance was quantitatively investigated in phantom and human model simulations, and experiments on five healthy volunteers inside the prostate. Comparisons were made with 16-channel microstrip line (16ML) and 10-channel fractionated dipole antenna (10DA) arrays. The 16LD was used to acquire anatomic and functional images of the prostate, kidneys, and heart.

RESULTS

The 16LD provided > 14% improvements in the signal-to-noise ratio (SNR), peak B1+, B1+ transmit, and SAR efficiencies over the 16ML and 10DA in simulations inside the prostate. Experimentally, the 16LD had > 20% higher SNR and B1+ transmit efficiency compared with other arrays, and achieved up to 51.8% higher peak B1+ compared with 10DA.

CONCLUSION

Combining loop and dipole elements provided a body imaging array with high channel count and density while limiting inter-element coupling. The 16LD improved both near and far-field performance compared with existing 7.0T body arrays and provided high-quality MRI of the prostate kidneys and heart. Magn Reson Med 77:884-894, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

(19)F MRSI of capecitabine in the liver at 7 T using broadband transmit-receive antennas and dual-band RF pulses

Capecitabine (Cap) is an often prescribed chemotherapeutic agent, successfully used to cure some patients from cancer or reduce tumor burden for palliative care. However, the efficacy of the drug is limited, it is not known in advance who will respond to the drug and it can come with severe toxicity. (19)F Magnetic Resonance Spectroscopy (MRS) and Magnetic Resonance Spectroscopic Imaging (MRSI) have been used to non-invasively study Cap metabolism in vivo to find a marker for personalized treatment. In vivo detection, however, is hampered by low concentrations and the use of radiofrequency (RF) surface coils limiting spatial coverage. In this work, the use of a 7T MR system with radiative multi-channel transmit-receive antennas was investigated with the aim of maximizing the sensitivity and spatial coverage of (19)F detection protocols. The antennas were broadband optimized to facilitate both the (1)H (298 MHz) and (19)F (280 MHz) frequencies for accurate shimming, imaging and signal combination. B1(+) simulations, phantom and noise measurements showed that more than 90% of the theoretical maximum sensitivity could be obtained when using B1(+) and B1(-) information provided at the (1)H frequency for the optimization of B1(+) and B1(-) at the (19)F frequency. Furthermore, to overcome the limits in maximum available RF power, whilst ensuring simultaneous excitation of all detectable conversion products of Cap, a dual-band RF pulse was designed and evaluated. Finally, (19)F MRS(I) measurements were performed to detect (19)F metabolites in vitro and in vivo. In two patients, at 10 h (patient 1) and 1 h (patient 2) after Cap intake, (19)F metabolites were detected in the liver and the surrounding organs, illustrating the potential of the set-up for in vivo detection of metabolic rates and drug distribution in the body.

Local specific absorption rate in brain tumors at 7 tesla

PURPOSE

MR safety at 7 Tesla relies on accurate numerical simulations of transmit electromagnetic fields to fully assess local specific absorption rate (SAR) safety. Numerical simulations for SAR safety are currently performed using models of healthy patients. These simulations might not be useful for estimating SAR in patients who have large lesions with potentially abnormal dielectric properties, e.g., brain tumors.

THEORY AND METHODS

In this study, brain tumor patient models are constructed based on scans of four patients with high grade brain tumors. Dielectric properties for the modeled tumors are assigned based on electrical properties tomography data for the same patients. Simulations were performed to determine SAR.

RESULTS

Local SAR increases in the tumors by as much as 30%. However, the location of the maximum 10-gram averaged SAR typically occurs outside of the tumor, and thus does not increase. In the worst case, if the tumor model is moved to the location of maximum electric field intensity, then we do observe an increase in the estimated peak 10-gram SAR directly related to the tumor.

CONCLUSION

Peak local SAR estimation made on the results of a healthy patient model simulation may underestimate the true peak local SAR in a brain tumor patient.

Radiofrequency configuration to facilitate bilateral breast (31) P MR spectroscopic imaging and high-resolution MRI at 7 Tesla

PURPOSE

High-resolution MRI combined with phospholipid detection may improve breast cancer grading. Currently, configurations are optimized for either high-resolution imaging or (31) P spectroscopy. To be able to perform both imaging as well as spectroscopy in a single session, we integrated a (1) H receiver array into a (1) H-(31) P transceiver at 7T. To ensure negligible signal loss due to coupling between elements, we investigated the use of a floating decoupling loop to enable bilateral MRI and (31) P MRS.

METHODS

Two quadrature double-tuned radiofrequency coils were designed for bilateral breast MR with active detuning at the (1) H frequency. The two coils were placed adjacent to each other and decoupled for both frequencies with a single resonant floating loop. Sensitivity of the bilateral configuration, facilitating space for a 26-element (1) H receive array, was compared with a transceiver configuration.

RESULTS

The floating loop was able to decouple the elements over 20 dB for both frequencies. Enlargement of the elements, to provide space for the receivers, and the addition of detuning electronics altered the (31) P sensitivity by 0.4 dB.

CONCLUSION

Dynamic contrast-enhanced scans of 0.7 mm isotropic, diffusion-weighted imaging, and (31) P MR spectroscopic imaging can be acquired at 7T in a single session as demonstrated in a patient with invasive ductal carcinoma.

Transmit and receive RF fields determination from a single low-tip-angle gradient-echo scan by scaling of SVD data

PURPOSE

A new method, called Transmit and Receive Patterns from Low-Tip-angle gradient-Echo Images (TRIPLET), is described which simultaneously maps the B1+ and B1- fields of a transmit/receive radiofrequency coil array. The input data are low-tip-angle gradient-echo images, which can be acquired in a relatively short scanning time.

THEORY AND METHODS

For each voxel in the field of view, a matrix can be assembled with the low-tip-angle gradient-echo image values of the radiofrequency coil array. Applying the singular value decomposition to those matrices, datasets are obtained which show a high resemblance with the true B1+ and B1- fields. These datasets are a voxel-wise scaled version of the true radiofrequency maps. The channel independent scaling parameters can be found by implicitly forcing the reconstructed fields to be solutions of the Maxwell equations. This is achieved by introducing a multipole expansion consisting of Bessel/Fourier functions.

RESULTS

Two FDTD simulated radiofrequency fields for two coil array combinations at 7 T and a measured, in vivo dataset at 7 T are investigated to illustrate the singular value decomposition analysis of the low-tip-angle gradient-echo images and to show how the B1+ and B1- fields can be reconstructed by Transmit and Receive Patterns from Low-Tip-angle gradient-Echo Images.

CONCLUSION

The Transmit and Receive Patterns from Low-Tip-angle gradient-Echo Images algorithm can convert the datasets from singular value decomposition analysis of low-tip-angle gradient-echo images to true B1+ and B1- fields.

Improved RF performance of travelling wave MR with a high permittivity dielectric lining of the bore

Application of travelling wave MR to human body imaging is restricted by the limited peak power of the available RF amplifiers. Nevertheless, travelling wave MR advantages like a large field of view excitation and distant location of transmit elements would be desirable for whole body MRI. In this work, improvement of the B1+ efficiency of travelling wave MR is demonstrated. High permittivity dielectric lining placed next to the scanner bore wall effectively reduces attenuation of the travelling wave in the longitudinal direction and at the same time directs the radial power flow toward the load. First, this is shown with an analytical model of a metallic cylindrical waveguide with the dielectric lining next to the wall and loaded with a cylindrical phantom. Simulations and experiments also reveal an increase of B1+ efficiency in the center of the bore for travelling wave MR with a dielectric lining. Phantom experiments show up to a 2-fold gain in B1+ with the dielectric lining. This corresponds to a 4-fold increase in power efficiency of travelling wave MR. In vivo experiments demonstrate an 8-fold signal-to-noise ratio gain with the dielectric lining. Overall, it is shown that dielectric lining is a constructive method to improve efficacy of travelling wave MR.

Specific absorption rate intersubject variability in 7T parallel transmit MRI of the head

Patient-specific radiofrequency shimming in high-field MRI strengthens the need for online, patient-specific specific absorption rate (SAR) monitoring. Numerical simulation is currently most effective for this purpose but may require a patient-specific dielectric model. To investigate whether a generic model may be combined with a safety factor to account for variation within the population, generic SAR behavior is studied for 7T MRI of the head. For six detailed head models, radiofrequency fields were simulated for an eight-channel parallel transmit array. SAR behavior is studied through comparison of the eigenvalues/eigenvectors of the local Q-matrices. Furthermore, numerical radiofrequency shimming experiments without and with SAR constraints were performed where SAR during optimization was evaluated on a generic model. In both cases, the ability of different generic models to predict actual SAR levels was evaluated. The largest eigenvalue distribution is comparable between models. Radiofrequency shimming without constraints improves the |B +1| homogeneity while the SAR increases substantially. Imposing constraints on SAR during optimization, estimating SAR on a generic model, was effective. A safety factor of 1.4 was found to be sufficient. Generic SAR behavior makes a generic head model a practical alternative to patient-specific models and allows effective |B +1| shimming with SAR constraints.

Characterization of transceive surface element designs for 7 tesla magnetic resonance imaging of the prostate: radiative antenna and microstrip

Ultra-high field magnetic resonance (≥7 tesla) imaging (MRI) faces challenges with respect to efficient spin excitation and signal reception from deeply situated organs. Traditional radio frequency surface coil designs relying on near-field coupling are suboptimal at high field strengths. Better signal penetration can be obtained by designing a radiative antenna in which the energy flux is directed to the target location. In this paper, two different radiative antenna designs are investigated to be used as transceive elements, which employ different dielectric permittivities for the antenna substrate. Their transmit and receive performances in terms of B(+)(1), local SAR (specific absorption rate) and SNR (signal-to-noise ratio) were compared using extensive electromagnetic simulations and MRI measurements with traditional surface microstrip coils. Both simulations and measurements demonstrated that the radiative element shows twofold gain in B(+)(1) and SNR at 10 cm depth, and additionally a comparable SAR peak value. In terms of transmit performance, the radiative antenna with a dielectric permittivity of 37 showed a 24% more favorable local SAR(10g avg)/(B(+)(1))(2) ratio than the radiative antenna with a dielectric permittivity of 90. In receive, the radiative element with a dielectric permittivity of 90 resulted in a 20% higher SNR for shallow depths, but for larger depths this difference diminished compared to the radiative element with a dielectric permittivity of 37. Therefore, to image deep anatomical regions effectively, the radiative antenna with a dielectric permittivity of 37 is favorable.

Fast dose calculation in magnetic fields with GPUMCD

A new hybrid imaging-treatment modality, the MRI-Linac, involves the irradiation of the patient in the presence of a strong magnetic field. This field acts on the charged particles, responsible for depositing dose, through the Lorentz force. These conditions require a dose calculation engine capable of taking into consideration the effect of the magnetic field on the dose distribution during the planning stage. Also in the case of a change in anatomy at the time of treatment, a fast online replanning tool is desirable. It is improbable that analytical solutions such as pencil beam calculations can be efficiently adapted for dose calculations within a magnetic field. Monte Carlo simulations have therefore been used for the computations but the calculation speed is generally too slow to allow online replanning. In this work, GPUMCD, a fast graphics processing unit (GPU)-based Monte Carlo dose calculation platform, was benchmarked with a new feature that allows dose calculations within a magnetic field. As a proof of concept, this new feature is validated against experimental measurements. GPUMCD was found to accurately reproduce experimental dose distributions according to a 2%-2 mm gamma analysis in two cases with large magnetic field-induced dose effects: a depth-dose phantom with an air cavity and a lateral-dose phantom surrounded by air. Furthermore, execution times of less than 15 s were achieved for one beam in a prostate case phantom for a 2% statistical uncertainty while less than 20 s were required for a seven-beam plan. These results indicate that GPUMCD is an interesting candidate, being fast and accurate, for dose calculations for the hybrid MRI-Linac modality.

Uniform prostate imaging and spectroscopy at 7 T: comparison between a microstrip array and an endorectal coil

An endorectal coil and an eight-element microstrip array were compared for prostate imaging at 7 T. An extensive radiofrequency safety assessment was performed with the use of finite difference time domain simulations to determine safe scan parameters. These simulations showed that the endorectal coil can deliver substantially more B(1)(+) to the prostate than can the microstrip array within the specific absorption rate safety guidelines. However, the B(1)(+) field of the endorectal coil is very inhomogeneous, which makes the use of adiabatic pulses compulsory for T(1) - or T(2) -weighted imaging. As a consequence, a full prostate examination is only possible in a feasible amount of time when the microstrip array is used for T(1) - and T(2) -weighted imaging, whereas the endorectal coil is required for spectroscopic imaging. The pulse parameters were optimised within the specific absorption rate guidelines and thereafter used to provide a good illustration of the possibilities of prostate imaging at 7 T.

Integrating a 1.5 T MRI scanner with a 6 MV accelerator: proof of concept

At the UMC Utrecht, The Netherlands, we have constructed a prototype MRI accelerator. The prototype is a modified 6 MV Elekta (Crawley, UK) accelerator next to a modified 1.5 T Philips Achieva (Best, The Netherlands) MRI system. From the initial design onwards, modifications to both systems were aimed to yield simultaneous and unhampered operation of the MRI and the accelerator. Indeed, the simultaneous operation is shown by performing diagnostic quality 1.5 T MRI with the radiation beam on. No degradation of the performance of either system was found. The integrated 1.5 T MRI system and radiotherapy accelerator allow simultaneous irradiation and MR imaging. The full diagnostic imaging capacities of the MRI can be used; dedicated sequences for MRI-guided radiotherapy treatments will be developed. This proof of concept opens the door towards a clinical prototype to start testing MRI-guided radiation therapy (MRIgRT) in the clinic.

Integrating a 1.5 T MRI scanner with a 6 MV accelerator: proof of concept

At the UMC Utrecht, The Netherlands, we have constructed a prototype MRI accelerator. The prototype is a modified 6 MV Elekta (Crawley, UK) accelerator next to a modified 1.5 T Philips Achieva (Best, The Netherlands) MRI system. From the initial design onwards, modifications to both systems were aimed to yield simultaneous and unhampered operation of the MRI and the accelerator. Indeed, the simultaneous operation is shown by performing diagnostic quality 1.5 T MRI with the radiation beam on. No degradation of the performance of either system was found. The integrated 1.5 T MRI system and radiotherapy accelerator allow simultaneous irradiation and MR imaging. The full diagnostic imaging capacities of the MRI can be used; dedicated sequences for MRI-guided radiotherapy treatments will be developed. This proof of concept opens the door towards a clinical prototype to start testing MRI-guided radiation therapy (MRIgRT) in the clinic.

Feasibility of MRI guided proton therapy: magnetic field dose effects

Many methods exist to improve treatment outcome in radiotherapy. Two of these are image-guided radiotherapy (IGRT) and proton therapy. IGRT aims at a more precise delivery of the radiation, while proton therapy is able to achieve more conformal dose distributions. In order to maximally exploit the sharp dose gradients from proton therapy it has to be combined with soft-tissue based IGRT. MRI-guided photon therapy (currently under development) offers unequalled soft-tissue contrast and real-time image guidance. A hybrid MRI proton therapy system would combine these advantages with the advantageous dose steering capacity of proton therapy. This paper addresses a first technical feasibility issue of this concept, namely the impact of a 0.5 T magnetic field on the dose distribution from a 90 MeV proton beam. In contrast to photon therapy, for MR-guided proton therapy the impact of the magnetic field on the dose distribution is very small. At tissue-air interfaces no effect of the magnetic field on the dose distribution can be detected. This is due to the low-energy of the secondary electrons released by the heavy protons.

MRI/linac integration

PURPOSE/OBJECTIVES

In radiotherapy the healthy tissue involvement still poses serious dose limitations. This results in sub-optimal tumour dose and complications. Daily image guided radiotherapy (IGRT) is the key development in radiation oncology to solve this problem. MRI yields superb soft-tissue visualization and provides several imaging modalities for identification of movements, function and physiology. Integrating MRI functionality with an accelerator can make these capacities available for high precision, real time IGRT.

DESIGN AND RESULTS

The system being built at the University Medical Center Utrecht is a 1.5T MRI scanner, with diagnostic imaging functionality and quality, integrated with a 6MV radiotherapy accelerator. The realization of a prototype of this hybrid system is a joint effort between the Radiotherapy Department of the University of Utrecht, the Netherlands, Elekta, Crawley, U.K., and Philips Research, Hamburg, Germany. Basically, the design is a 1.5 T Philips Achieva MRI scanner with a Magnex closed bore magnet surrounded by a single energy (6 MV) Elekta accelerator. Monte Carlo simulations are used to investigate the radiation beam properties of the hybrid system, dosimetry equipment and for the construction of patient specific dose deposition kernels in the presence of a magnetic field. The latter are used to evaluate the IMRT capability of the integrated MRI linac.

CONCLUSIONS

A prototype hybrid MRI/linac for on-line MRI guidance of radiotherapy (MRIgRT) is under construction. The aim of the system is to deliver the radiation dose with mm precision based on diagnostic quality MR images.