Performance of receive head arrays versus ultimate intrinsic SNR at 7 T and 10.5 T.

PURPOSE: We examined magnetic field dependent SNR gains and ability to capture them with multichannel receive arrays for human head imaging in going from 7 T, the most commonly used ultrahigh magnetic field (UHF) platform at the present, to 10.5 T, which represents the emerging new frontier of >10 T in UHFs. METHODS: Electromagnetic (EM) models of 31-channel and 63-channel multichannel arrays built for 10.5 T were developed for 10.5 T and 7 T simulations. A 7 T version of the 63-channel array with an identical coil layout was also built. Array performance was evaluated in the EM model using a phantom mimicking the size and electrical properties of the human head and a digital human head model. Experimental data was obtained at 7 T and 10.5 T with the 63-channel array. Ultimate intrinsic SNR (uiSNR) was calculated for the two field strengths using a voxelized cloud of dipoles enclosing the phantom or the digital human head model as a reference to assess the performance of the two arrays and field depended SNR gains. RESULTS: uiSNR calculations in both the phantom and the digital human head model demonstrated SNR gains at 10.5 T relative to 7 T of 2.6 centrally, ˜2 at the location corresponding to the edge of the brain, ˜1.4 at the periphery. The EM models demonstrated that, centrally, both arrays captured ˜90% of the uiSNR at 7 T, but only ˜65% at 10.5 T, leading only to ˜2-fold gain in array SNR in going from 7 to 10.5 T. This trend was also observed experimentally with the 63-channel array capturing a larger fraction of the uiSNR at 7 T compared to 10.5 T, although the percentage of uiSNR captured were slightly lower at both field strengths compared to EM simulation results. CONCLUSIONS: Major uiSNR gains are predicted for human head imaging in going from 7 T to 10.5 T, ranging from ˜2-fold at locations corresponding to the edge of the brain to 2.6-fold at the center, corresponding to approximately quadratic increase with the magnetic field. Realistic 31- and 63-channel receive arrays, however, approach the central uiSNR at 7 T, but fail to do so at 10.5 T, suggesting that more coils and/or different type of coils will be needed at 10.5 T and higher magnetic fields.

A workflow for predicting radiofrequency-induced heating around bilateral deep brain stimulation electrodes in MRI.

BACKGROUND: Heating around deep brain stimulation (DBS) in magnetic resonance imaging (MRI) occurs when the time-varying electromagnetic (EM) fields induce currents in the electrodes which can generate heat and potentially cause tissue damage. Predicting the heating around the electrode contacts is important to ensure the safety of patients with DBS implants undergoing an MRI scan. We previously proposed a workflow to predict heating around DBS contacts and introduced a parameter, equivalent transimpedance, that is independent of electrode trajectories, termination, and radiofrequency (RF) excitations. The workflow performance was validated in a unilateral DBS system. PURPOSE: To predict RF heating around the contacts of bilateral (DBS) electrodes during an MRI scan in an anthropomorphic head phantom. METHODS: Bilateral electrodes were fixed in a skull phantom filled with hydroxyethyl cellulose (HEC) gel. The electrode shafts were suspended extracranially, in a head and torso phantom filled with the same gel material. The current induced on the electrode shaft was experimentally measured using an MR-based technique 3 cm above the tip. A transimpedance value determined in a previous offline calibration was used to scale the shaft current and calculate the contact voltage. The voltage was assigned as a boundary condition on the electrical contacts of the electrode in a quasi-static (EM) simulation. The resulting specific absorption rate (SAR) distribution became the input for a transient thermal simulation and was used to predict the heating around the contacts. RF heating experiments were performed for eight different lead trajectories using circularly polarized (CP) excitation and two linear excitations for one trajectory. The measured temperatures for all experiments were compared with the simulated temperatures and the root-mean-squared errors (RMSE) were calculated. RESULTS: The RF heating around the contacts of both bilateral electrodes was predicted with ≤ 0.29°C of RMSE for 20 heating scenarios. CONCLUSION: The workflow successfully predicted the heating for different bilateral DBS trajectories and excitation patterns in an anthropomorphic head phantom.

Implant-friendly MRI of deep brain stimulation electrodes at 7 T.

PURPOSE: The purpose of this study is to present a strategy to calculate the implant-friendly (IF) excitation modes-which mitigate the RF heating at the contacts of deep brain stimulation (DBS) electrodes-of multichannel RF coils at 7 T. METHODS: An induced RF current on an implantable electrode generates a scattered magnetic field whose left-handed circularly polarizing component ( B 1 + $$ B{1}^{+} $$ ) is approximated using a B 1 + $$ B{1}^{+} $$ -mapping technique and subsequently used as a gauge for the electrode's induced current. Using this approach, the relative induced currents resulting from each channel of a multichannel RF coil on the DBS electrode were calculated. The IF modes of the corresponding multichannel coil were determined by calculating the null space of the relative induced currents. The proposed strategy was tested and validated for unilateral and bilateral commercial DBS electrodes (directional lead; Infinity DBS system, Abbott Laboratories) placed inside a uniform phantom by performing heating and imaging studies on a 7T MRI scanner using a 16-channel transceive RF coil. RESULTS: Neither individual IF modes nor shim solutions obtained from IF modes induced significant temperature increase when used for a high-power turbo spin-echo sequence. In contrast, shimming with the scanner's toolbox (i.e., based on per-channel B 1 + $$ B{1}^{+} $$ fields) resulted in a more than 2°C temperature increase for the same amount of input power. CONCLUSION: A strategy for calculating the IF modes of a multichannel RF coil is presented. This strategy was validated using a 16-channel RF coil at 7 T for unilateral and bilateral commercial DBS electrodes inside a uniform phantom.

Effect of field strength on RF power deposition near conductive leads: A simulation study of SAR in DBS lead models during MRI at 1.5 T-10.5 T.

BACKGROUND: Since the advent of magnetic resonance imaging (MRI) nearly four decades ago, there has been a quest for ever-higher magnetic field strengths. Strong incentives exist to do so, as increasing the magnetic field strength increases the signal-to-noise ratio of images. However, ensuring patient safety becomes more challenging at high and ultrahigh field MRI (i.e., ≥3 T) compared to lower fields. The problem is exacerbated for patients with conductive implants, such as those with deep brain stimulation (DBS) devices, as excessive local heating can occur around implanted lead tips. Despite extensive effort to assess radio frequency (RF) heating of implants during MRI at 1.5 T, a comparative study that systematically examines the effects of field strength and various exposure limits on RF heating is missing. PURPOSE: This study aims to perform numerical simulations that systematically compare RF power deposition near DBS lead models during MRI at common clinical and ultra-high field strengths, namely 1.5, 3, 7, and 10.5 T. Furthermore, we assess the effects of different exposure constraints on RF power deposition by imposing limits on either the B1+ or global head specific absorption rate (SAR) as these two exposure limits commonly appear in MRI guidelines. METHODS: We created 33 unique DBS lead models based on postoperative computed tomography (CT) images of patients with implanted DBS devices and performed electromagnetic simulations to evaluate the SAR of RF energy in the tissue surrounding lead tips during RF exposure at frequencies ranging from 64 MHz (1.5 T) to 447 MHz (10.5 T). The RF exposure was implemented via realistic MRI RF coil models created based on physical prototypes built in our institutions. We systematically examined the distribution of local SAR at different frequencies with the input coil power adjusted to either limit the B1+ or the global head SAR. RESULTS: The MRI RF coils at higher resonant frequencies generated lower SARs around the lead tips when the global head SAR was constrained. The trend was reversed when the constraint was imposed on B1+. CONCLUSION: At higher static fields, MRI is not necessarily more dangerous than at lower fields for patients with conductive leads. Specifically, when a conservative safety criterion, such as constraints on the global SAR, is imposed, coils at a higher resonant frequency tend to generate a lower local SAR around implanted leads due to the decreased B1+ and, by proxy, E field levels.

Tier-based formalism for safety assessment of custom-built radio-frequency transmit coils.

The purpose of this work is to propose a tier-based formalism for safety assessment of custom-built radio-frequency (RF) coils that balances validation effort with the effort put in determinating the safety factor. The formalism has three tier levels. Higher tiers require increased effort when validating electromagnetic simulation results but allow for less conservative safety factors. In addition, we propose a new method to calculate modeling uncertainty between simulations and measurements and a new method to propagate uncertainties in the simulation into a safety factor that minimizes the risk of underestimating the peak specific absorption rate (SAR). The new safety assessment procedure was completed for all tier levels for an eight-channel dipole array for prostate imaging at 7 T and an eight-channel dipole array for head imaging at 10.5 T, using data from two different research sites. For the 7 T body array, the validation procedure resulted in a modeling uncertainty of 77% between measured and simulated local SAR distributions. For a situation where RF shimming is performed on the prostate, average power limits of 2.4 and 4.5 W/channel were found for tiers 2 and 3, respectively. When the worst-case peak SAR among all phase settings was calculated, power limits of 1.4 and 2.7 W/channel were found for tiers 2 and 3, respectively. For the 10.5 T head array, a modeling uncertainty of 21% was found based on B1 + mapping. For the tier 2 validation, a power limit of 2.6 W/channel was calculated. The demonstrated tier system provides a strategy for evaluating modeling inaccuracy, allowing for the rapid translation of novel coil designs with conservative safety factors and the implementation of less conservative safety factors for frequently used coil arrays at the expense of increased validation effort.

A workflow for predicting temperature increase at the electrical contacts of deep brain stimulation electrodes undergoing MRI.

PURPOSE: The purpose of this study is to present a workflow for predicting the radiofrequency (RF) heating around the contacts of a deep brain stimulation (DBS) lead during an MRI scan. METHODS: The induced RF current on the DBS lead accumulates electric charge on the metallic contacts, which may cause a high local specific absorption rate (SAR), and therefore, heating. The accumulated charge was modeled by imposing a voltage boundary condition on the contacts in a quasi-static electromagnetic (EM) simulation allowing thermal simulations to be performed with the resulting SAR distributions. Estimating SAR and temperature increases from a lead in vivo through EM simulation is not practical given anatomic differences and variations in lead geometry. To overcome this limitation, a new parameter, transimpedance, was defined to characterize a given lead. By combining the transimpedance, which can be measured in a single calibration scan, along with MR-based current measurements of the lead in a unique orientation and anatomy, local heating can be estimated. Heating determined with this approach was compared with results from heating studies of a commercial DBS electrode in a gel phantom with different lead configurations to validate the proposed method. RESULTS: Using data from a single calibration experiment, the transimpedance of a commercial DBS electrode (directional lead, Infinity DBS system, Abbott Laboratories, Chicago, IL) was determined to be 88 Ω. Heating predictions using the DBS transimpedance and rapidly acquired MR-based current measurements in 26 different lead configurations resulted in a <23% (on average 11.3%) normalized root-mean-square error compared to experimental heating measurements during RF scans. CONCLUSION: In this study, a workflow consisting of an MR-based current measurement on the DBS lead and simple quasi-static EM/thermal simulations to predict the temperature increase around a DBS electrode undergoing an MRI scan is proposed and validated using a commercial DBS electrode.

A nine-channel transmit/receive array for spine imaging at 10.5 T: Introduction to a nonuniform dielectric substrate antenna.

PURPOSE: The purpose of this study is to introduce a new antenna element with improved transmit performance, named the nonuniform dielectric substrate (NODES) antenna, for building transmit arrays at ultrahigh-field. METHODS: We optimized a dipole antenna at 10.5 Tesla by maximizing the B1+ -SAR efficiency in a phantom for a human spine target. The optimization parameters included permittivity variation in the substrate, substrate thickness, antenna length, and conductor geometry. We conducted electromagnetic simulations as well as phantom experiments to compare the transmit/receive performance of the proposed NODES antenna design with existing coil elements from the literature. RESULTS: Single NODES element showed up to 18% and 30% higher B1+ -SAR efficiency than the fractionated dipole and loop elements, respectively. The new element is substantially shorter than a commonly used dipole, which enables z-stacked array formation; it is additionally capable of providing a relatively uniform current distribution along its conductors. The nine-channel transmit/receive NODES array achieved 7.5% higher B1+ homogeneity than a loop array with the same number of elements. Excitation with the NODES array resulted in 33% lower peak 10g-averaged SAR and required 34% lower input power than the loop array for the target anatomy of the spine. CONCLUSION: In this study, we introduced a new RF coil element: the NODES antenna. NODES antenna outperformed the widely used loop and dipole elements and may provide improved transmit/receive performance for future ultrahigh field MRI applications.

Analysis and mitigation of noise in simultaneous transmission and reception in MRI.

PURPOSE: In simultaneous transmission and reception (STAR) MRI, along with the coupling of the excitation pulse to the received signal, noise, and undesired distortions (spurs) coming from the transmit chain also leak into the acquired signal and degrade image quality. Here, properties of this coupled noise and its relationship with the transmit amplifier gain, transmit chain noise density, isolation performance, and imaging bandwidth are analyzed. It is demonstrated that by utilizing a recently proposed STAR technique, the transmit noise can be reduced. The importance of achieving high isolation and careful selection of the corresponding parameters are demonstrated. THEORY AND METHODS: A cancellation algorithm, together with a vector modulator, is used for transmit-receive isolation. The scanner is modeled as a pipeline of blocks to demonstrate the noise contribution from each block. With higher isolation, coupled transmit noise can be reduced to the point that the dominant noise source becomes acquisition noise, as in the case for pulsed MRI. Amplifiers with different gain and noise properties are used in the experiments to verify the derived noise-transmit parameter relation. RESULTS: With the proposed technique, more than 80 dB isolation in the analog domain is achieved. The leakage noise and the spurs coupled from the transmit chain, are reduced. It is shown that the transmit gain plays the most critical role in determining sufficient isolation, whereas the amplifier noise figure does not contribute as much. CONCLUSION: The transmit noise and the spurs in STAR imaging are analyzed and mitigated by using a vector modulator.

Eigenmode analysis of the scattering matrix for the design of MRI transmit array coils.

PURPOSE: To obtain efficient operation modes of transmit array (TxArray) coils using a general design technique based on the eigenmode analysis of the scattering matrix. METHODS: We introduce the concept of modal reflected power and excitation eigenmodes, which are calculated as the eigenvalues and eigenvectors of SH S, where the superscript H denotes the Hermitian transpose. We formulate the normalized reflected power, which is the ratio of the total reflected power to the total incident power of TxArray coils for a given excitation signal as the weighted sum of the modal reflected power. By minimizing the modal reflected power of TxArray coils, we increase the excitation space with a low total reflection. The algorithm was tested on 4 dual-row TxArray coils with 8 to 32 channels. RESULTS: By minimizing the modal reflected power, we designed an 8-element TxArray coil to have a low reflection for 7 out of 8 dimensions of the excitation space. Similarly, the minimization of the modal reflected power of a 16-element TxArray coil enabled us to enlarge the dimension of the excitation space by 50% compared with commonly employed design techniques. Moreover, we demonstrated that the low total reflected power for some critical excitation modes, such as the circularly polarized mode, can be achieved for all TxArray coils even with a high level of coupling. CONCLUSION: Eigenmode analysis is an efficient method that intuitively provides a quantitative and compact representation of the coil's power transmission capabilities. This method also provides insight into the excitation modes with low reflection.

Accelerating the co-simulation method for the design of transmit array coils for MRI.

OBJECTIVE: Accelerating the co-simulation method for the design of transmit array (TxArray) coils is studied using equivalent circuit models. MATERIALS AND METHODS: Although the co-simulation method dramatically reduces the complexity of the design of TxArray coils, finding the optimum solution is not trivial since there exist many local minima in the optimization problem. We propose to utilize an equivalent circuit model of the TxArray coil to obtain a proper initial guess for the optimization process of the co-simulation method. To prove the concept, six different TxArray coils (i.e., three degenerate birdcage coils (DBC), two dual-row head coils, and one elliptical body TxArray coil) with two different loading strategies (cylindrical phantom and human head/body model) at 3 T field strength are investigated theoretically; as an example study, an eight-channel head-DBC is constructed using the obtained values. RESULTS: This approach accelerates the design process more than 20-fold for the coils that are investigated in this manuscript. CONCLUSION: A fast and accurate method for tuning and decoupling of a TxArray coil can be achieved using its equivalent circuit model combined with the co-simulation method.

Improving radiofrequency power and specific absorption rate management with bumped transmit elements in ultra-high field MRI.

PURPOSE: In this study, we investigate a strategy to reduce the local specific absorption rate (SAR) while keeping B1+ constant inside the region of interest (ROI) at the ultra-high field (B0 ≥ 7T) MRI. METHODS: Locally raising the resonance structure under the discontinuity (i.e., creating a bump) increases the distance between the accumulated charges and the tissue. As a result, it reduces the electric field and local SAR generated by these charges inside the tissue. The B1+ at a point that is sufficiently far from the coil, however, is not affected by this modification. In this study, three different resonant elements (i.e., loop coil, snake antenna, and fractionated dipole [FD]) are investigated. For experimental validation, a bumped FD is further investigated at 10.5T. After the validation, the transmit performances of eight-channel arrays of each element are compared through electromagnetic (EM) simulations. RESULTS: Introducing a bump reduced the peak 10g-averaged SAR by 21, 26, 23% for the loop and snake antenna at 7T, and FD at 10.5T, respectively. In addition, eight-channel bumped FD array at 10.5T had a 27% lower peak 10g-averaged SAR in a realistic human body simulation (i.e., prostate imaging) compared to an eight-channel FD array. CONCLUSION: In this study, we investigated a simple design strategy based on adding bumps to a resonant element to reduce the local SAR while maintaining B1+ inside an ROI. As an example, we modified an FD and performed EM simulations and phantom experiments with a 10.5T scanner. Results show that the peak 10g-averaged SAR can be reduced more than 25%.

In vivo human head MRI at 10.5T: A radiofrequency safety study and preliminary imaging results.

PURPOSE: The purpose of this study is to safely acquire the first human head images at 10.5T. METHODS: To ensure safety of subjects, we validated the electromagnetic simulation model of our coil. We obtained quantitative agreement between simulated and experimental B1+ and specific absorption rate (SAR). Using the validated coil model, we calculated radiofrequency power levels to safely image human subjects. We conducted all experiments and imaging sessions in a controlled radiofrequency safety lab and the whole-body 10.5T scanner in the Center for Magnetic Resonance Research. RESULTS: Quantitative agreement between the simulated and experimental results was obtained including S-parameters, B1+ maps, and SAR. We calculated peak 10 g average SAR using 4 different realistic human body models for a quadrature excitation and demonstrated that the peak 10 g SAR variation between subjects was less than 30%. We calculated safe power limits based on this set and used those limits to acquire T2 - and T2∗ -weighted images of human subjects at 10.5T. CONCLUSIONS: In this study, we acquired the first in vivo human head images at 10.5T using an 8-channel transmit/receive coil. We implemented and expanded a previously proposed workflow to validate the electromagnetic simulation model of the 8-channel transmit/receive coil. Using the validated coil model, we calculated radiofrequency power levels to safely image human subjects.