Implant‐friendly MRI of deep brain stimulation electrodes at 7 T

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.


INTRODUCTION
2][3][4][5][6] MRI examinations are essential for these patients either as a standard diagnostic tool or for monitoring their existing condition.However, MRI conditions apply to DBS patients who can be very restrictive in terms of specific absorption rate 7 due to the significant safety risks.These risks are primarily related to induced RF currents on elongated metallic leads, which may critically increase the temperature at the electrode's tip and cause damage to the surrounding tissues 8,9 during the RF excitation.Therefore, it is crucial to develop strategies that mitigate the induced RF currents on DBS electrodes and RF heating of these devices.1][12][13] Despite the promising impact of the trajectory management techniques on the reduction of RF heating at the tip, these techniques may not be a practical solution to all cases, such as patients with existing DBS implants.8][19][20][21][22][23][24] Eryaman et al. 17 theoretically demonstrated that implant-friendly (IF) RF excitation is possible by generating a plane with minimal electric field (E-field) without significant alteration in the transmit (Tx) magnetic field.They showed that the linear excitation mode of a birdcage coil could be an example of such IF excitation.Kazemivalipour et al. 20 used a rotating linearly polarized birdcage coil to steer the low-E-field slab-like region within the head.They numerically and experimentally showed that the induced current on the electrode and RF heating at the tip can be significantly decreased when the electrode is maximally contained within this low-E-field region.
2][23][24][25] Guerin et al. 22 presented a pTx pulse design technique that minimizes the RF excitation inhomogeneity while constraining the specific absorption rate in the tissue surrounding the tip of the DBS electrode.For five different DBS patient models, they numerically showed that their proposed pTx strategy used with an eight-channel coil at 3 T can reduce the tip heating by a factor of 18 on average compared with the commonly used birdcage coil.Etezadi-Amoli et al. 24 introduced null-mode excitation, which minimizes the induced RF current on the electrode.They performed a feasibility study at 1.5 T on a single wire using a four-channel transmit array coil and demonstrated that the heating at the tip of the wire can be reduced below the detectable level.
Despite the extensive investigation at 1.5 T and 3 T, the potential for pTx to minimize implant heating has not been investigated in detail at ultrahigh fields (defined as ≥ 7 T).As for the importance of ultrahigh-field MRI for DBS patients, various significant benefits of 7T MRI for DBS procedures have recently been shown. 26,27Duchin et al. 26 used 7T MR images to create patient-specific anatomical 3D models of subthalamic nuclei to enhance presurgical DBS targeting as well as postsurgical visualization of the DBS lead position and orientation.Schrock et al. 27 took advantage of 7T MRI techniques to detect the mispositioning of an implant in a patient with severe side effects, and subsequently reverse those side effects by repositioning the implant.Moreover, contrary to popular belief, Kazemivalipour et al. 28 numerically demonstrated that higher-field MRI is not necessarily more dangerous than lower fields for patients with conductive leads.Additionally, with the recent Food and Drug Administration approval of multiple 7T MRI scanners for clinical use, 29,30 investigating IF strategies at 7 T carries a high implication.
In this study, motivated by the benefits of 7T MRI for DBS patients 26,27 and the success of previously proposed IF strategies, 17,[20][21][22][23][24][25] implant-friendly excitation solutions are investigated at 7 T.For that purpose, the IF excitation modes of a 16-channel pTx coil 31 are calculated and used to mitigate the RF heating at the tip of a commercial DBS electrode (directional lead, Infinity DBS system; Abbott Laboratories, Chicago, IL, USA).The proposed strategy is experimentally tested and validated by measuring the temperature increase (ΔT) at the tip of both unilateral and bilateral electrodes in a head-shaped uniform phantom undergoing an MRI scan at 7 T.

THEORY
Previously, the analytical expression that relates the current, I ind , induced by an RF coil on an oblique DBS electrode (Figure 1A) to the left-handed circularly polarized component (defined as B1 + ) of the scattered magnetic field on an arbitrary transversal plane has been calculated as follows 32 : where B1 +sca is the scattered B1 + field caused by the induced current on the electrode; r 0 is the radial distance of an arbitrary observation point from the electrode; and  0 denotes the azimuthal angle of the observation point in relation to the conventional MRI scanner coordinate system (i.e., laboratory frame).As shown in Figure 1A, the rotation and obliqueness of the electrode with respect to the primary coordinate system are presented by a polar angle,  r , and an azimuthal angle,  r .In the case of an N-channel pTx RF coil, each channel can be treated as an independent transmitting RF coil.Therefore, the scattered B1 + at the observation point, B1 +sca  (r 0 ,  0 ), can be related to the induced current on the electrode, I ind  , due to the jth channel excitation through Eq. (1).This equation suggests that at any arbitrary observation point, the jth scattered B1 + field only varies with the electrode's current induced by the corresponding channel, as follows: On the one hand, the total induced current on the electrode caused by an arbitrary excitation of the pTx RF coil, I ind exc , can be expanded in the form of a superposition of induced currents due to single-channel excitations, as follows: where   denotes the complex weight of the jth transmit channel.Equation (3) can be written in the form of matrix-vector multiplication as follows: where I ind is a 1 × N row vector consisting of induced currents, I ind  ; and  is an N × 1 column vector consisting of complex weights,   .
On the other hand, previous studies 32,33 show that the induced current on the shaft of an electrode (near the tip) reliably controls the temperature increase (ΔT) at the tip of the electrode.Therefore, any complex excitation vector,  null , that is orthogonal to the row vector I ind sets Eq. ( 4) to zero and causes no heating at the tip of the electrode.Such an excitation vector can be considered an IF excitation mode for the electrode under examination.In other words, the null space 33 of the row vector I ind (i.e.,  null = null ( I ind ) )-which is an N×(N−1) matrixprovides all N−1 linearly independent IF modes of an N-channel pTx coil for a particular electrode.Note that because the row vector I ind is considered a rank-deficient matrix with the rank of 1, according to the rank-nullity theorem 34 (i.e., rank[I ind ] + nullity[I ind ] = N), there exist N−1 linearly independent vectors,  null , that set Eq. (4) to zero.These vectors,  null , are the columns of the matrix  null such that (5)

METHODS
Here, we have proposed a workflow (Figure 2) for determining the IF excitation modes of an N-channel pTx coil that can be used to safely image bilateral DBS electrodes at 7 T. Nevertheless, this workflow can simply be extended to any number of DBS electrodes if the number of transmit channels is greater than the number of electrodes (i.e., due to the rank-nullity theorem 34 ).The proposed workflow was experimentally validated using commercial DBS electrodes (directional lead, Infinity DBS system, Abbott Laboratories, Chicago, IL, USA) placed unilaterally and bilaterally inside a uniform head-shaped gel phantom (14 g/L hydroxyethyl cellulose, 2.25 g/L NaCl, and 0.25 g/L CuSO 4 with electrical properties of ε r = 78 and σ = 0.46S/m).All experiments were carried out on a 7T whole-body MRI scanner (Magnetom; Siemens Healthineers, Erlangen, Germany) using a 16-channel transmission line transceiver (Tx/Rx) array RF coil. 31Each channel of the RF coil was powered by a 1-kW amplifier (Communication Power Corp., Hauppauge, NY, USA) controlled by a remotely operated 16-channel RF phase/amplitude gain unit.

IF mode calculation workflow
Figure 2 shows our previously proposed IF mode calculation workflow 35 extended from unilateral to bilateral electrode configuration.In the extended version, the IF modes for one of the electrodes were first calculated.Then, the latter modes were used in a similar process to calculate IF modes for the second electrode.Throughout this workflow, it was assumed that the incident B1 + (around the shaft on plane P 2 , where the current is measured) is approximately equal to the total B1 + underneath the tip of the electrode on plane P 1 (Figure 1C,D). 32This approximation is based on the rationale that the electrode current Workflow for calculating the implant-friendly excitation modes of an N-channel parallel-transmit (pTx) coil for safe MRI of bilateral DBS electrodes.
has almost zero impact on the B1 + field at a plane located underneath the tip.The distance between the tip and shaft is short enough (∼30 mm) relative to the wavelength (∼110-130 mm) inside the human body (with relative permittivity in the range of 60-80) at 7 T to assume a uniform incident B1 + field from the tip to the shaft.The proposed workflow (Figure 2) can be laid out in five steps, detailed as follows: 1 B 1 + maps: First, complex relative B1 + maps corresponding to each Tx channel were acquired on two axial planes-one plane (P 2 ) intersecting one of the electrodes at the shaft and the other plane (P 1 ) underneath the tip of the same electrode (Figure 1C,D)-using consecutive fast, low-flip-angle gradient-echo sequences (see Section 3.2.1). ) caused by each of the N transmit channels were calculated, the null space of the row vector I ind was calculated to find all solutions that set Eq. ( 4) to zero.The calculated null space,  null , contained N-1 column vectors, which were the IF excitation modes of the pTx coil under investigation for one of the DBS electrodes.5 IF modes: Considering that this system used a bilateral DBS configuration, a similar routine was repeated for the second electrode with a modification in the definition of per-channel excitations.At the conclusion of the first routine (Step 4), N−1 IF excitation modes were computed for the first electrode.In the subsequent procedure, these modes were used as per-channel excitations for the second electrode.Eventually, (N−2) IF excitation modes were calculated, which set Eq. ( 4) to zero at the shaft of both electrodes.

Experimental validation
For both unilateral and bilateral electrode configurations, three types of experiments were conducted: (1) acquisition of B1 + maps to calculate IF modes, (2) validation of IF modes using temperature measurement, and (3) evaluation of IF modes with RF shimming.

B 1 + maps: Per-channel versus IF modes
As laid out in Step 1 of this workflow, complex B1 + maps on two axial planes (P 1 and P 2 ; Figure 1C,D) were acquired and postprocessed to obtain the IF modes.For this purpose, a fast, low-power, relative B1 + mapping technique 36 was exploited.This technique consisted of 16 (number of pTx coil channels) consecutive small-tip-angle gradient-echo scans with TR/TE = 4500/1.52ms, FA = 5 • , in-plane resolution = 1 mm, slice thickness = 4 mm, acquisition matrix = 200 × 200, no. of averages = 1, pixel bandwidth = 735 Hz/pixel, GRAPPA acceleration factor = 2, and acquisition time = 36 s.During each scan, only a single Tx channel was active, whereas the signal was received from all 16 receive channels. 36e B1 + field generated by the induced RF current on the electrode interferes with the incident B1 + field generated by the Tx coil and causes an excessive B1 + inhomogeneity in the proximity of the electrode. 33In theory, an IF mode excitation should be free of this effect because induced currents on the electrode are eliminated.To explore the efficacy of the IF mode calculations using this theory, B1 + maps corresponding to IF mode excitations on the axial plane P 2 (Figure 1C,D) were acquired and compared with B1 + maps corresponding to single-channel excitations on the same plane.

Temperature measurement: IF modes
To validate the calculated IF modes, which intended to eliminate the RF heating around the tip of the electrodes, the progression of temperature change at the tip of the DBS electrodes was measured while the experimental setup was exposed to RF energy using the pTx coil excited with an IF mode.The temperature measurement was performed using fiber optic temperature probes (Lumasense Technologies, Santa Clara, CA, USA) and repeated for all calculated IF modes.
For the unilateral electrode experiment, a total average power of 6 W was applied using a 1-ms-square RF pulse with a duty cycle of 5%.However, in the case of the bilateral electrodes, a 3D turbo spin-echo (TSE) sequence with FA = 100 • , TR/TE = 300/105 ms, in-plane resolution = 0.35 mm, slice thickness = 1.25 mm, acquisition matrix = 576 × 576, TSE factor = 15, no. of averages = 1, and pixel bandwidth = 225 Hz/pixel was used, with a total average power of 12 W.
To demonstrate the potential risks of DBS electrodes undergoing MRI scans, the unilateral electrode was exposed to RF energy using a 6-W square pulse with four different arbitrary excitations, and the ΔT was measured.For the same purpose, the worst-case heating scenario was tested using the complex conjugate of the relative per-channel induced currents as elements of the excitation vector.The ΔT at the tip of the electrode corresponding to this excitation was also measured.Similarly, in the bilateral electrode configuration, the worst-case heating scenarios for each electrode were calculated and used in the temperature measurement experiments (i.e., the 12-W TSE sequence).

RF shimming: Conventional versus IF modes
To evaluate the capability of the IF modes in RF shimming applications, a clinically relevant region of interest (ROI) (i.e., representing a ROI in the frontal lobe) was chosen.RF shimming with the goal of minimum inhomogeneity within the ROI was performed using the calculated IF modes, as well as the scanner's shimming toolbox.The IF mode shimming was performed by solving the following optimization problem using MATLAB's fminsearch function with the Nelder-Mead optimization algorithm, as follows: where B1 + ROI is an Npxl × N matrix whose columns consist of per-channel complex B1 + values on Npxl pixels within the ROI; x is an (N−L) × 1 vector (i.e., L = number of electrodes) and represents complex weights of IF modes; ||.|| denotes the  2 -norm of the vector it operates on; and CoV is the coefficient of variation defined as In ).In all cases, ΔT was measured.Also, to demonstrate the potential of the IF strategy in imaging the tissues around electrodes in DBS patients, inhomogeneity of the TSE images around electrodes (i.e., within a circle with a diameter of 30 mm) on Plane P 2 was quantitatively compared using the CoV.

IF modes
Figure 3 shows the amplitudes and phases of 15 excitation vectors-corresponding to the 16-channel pTx coil 31 -that produce IF modes for the unilateral electrode configuration.The green diamond on each plot corresponds to the predominant channel within that particular IF mode.Similarly, Figure 4 represents the 14 IF modes of this coil for the bilateral electrode configuration.

B 1 + maps: Per-channel versus IF modes
Using these IF excitations, a B1 + map corresponding to each mode was acquired on Plane P 2 and compared with that of the per-channel excitation.Figures 5 and 6 show the outcome of this comparative study for unilateral and bilateral electrode configurations, respectively.

4.2.2
Temperature measurement: IF modes Using a square RF pulse with a 6-W total average power, the unilateral electrode setup was exposed to RF energy.The ΔT at the tip of the electrode was recorded for all 15 IF modes, as shown in Figure 7A.The bilateral electrode setup had undergone an MRI scan using a TSE pulse sequence with a total average power of 12 W, and the ΔT at the tip of both electrodes corresponding to 14 IF excitation modes were measured and are presented in Figure 8A.
In addition, the ΔT was measured in the worst-case heating scenarios for both unilateral and bilateral configurations, and they are demonstrated in Figures 7B and 8B, respectively.Furthermore, Figure 7C shows the ΔT at the tip of the electrode in the unilateral configuration for four arbitrary excitation vectors.

RF shimming: Conventional versus IF modes
For both unilateral and bilateral configurations, a rectangular ROI was chosen in the right anterior of the head-shaped phantom, as respectively indicated in Figures 9 and 10.Using the calculated IF modes, as well as the scanner's shimming toolbox, RF shimming was performed within the ROI.These RF shimming scenarios resulted in different excitation vectors, which were used for 2D-TSE and 3D-TSE imaging of the unilateral and bilateral configuration, respectively.Figure 9 shows the TSE images as well as the ΔT at the tip, corresponding to the unilateral electrode configuration for both shimming scenarios.Although excitation using the scanner-calculated shimming solution caused a ΔT of about 2.6 • C, the IF excitation did not cause any measurable ΔT at contacts of the electrode.Moreover, the CoV around the shaft of the electrode was calculated to be 49% and 31% for the scanner's shimming solution and IF excitation, respectively.

F I G U R E 3
Amplitudes (blue bars) and phases (red bars) of 15 implant-friendly (IF) excitation vectors corresponding to the 16-channel parallel-transmit coil for the unilateral electrode configuration.The green diamond on each plot corresponds to the predominant channel within that particular IF mode.
Similarly, Figure 10 summarizes the results for the bilateral electrode configuration.In this case, the highest ΔT among two electrodes was measured as about 2 • C and 0.5 • C for scanner-calculated and IF shimming scenarios, respectively.The CoV around the corresponding electrode in the two shimming scenarios was calculated as 41% and 27%, respectively.

DISCUSSION
In this study, a new strategy was proposed to calculate the implant-friendly modes of a pTx coil for unilateral and bilateral DBS electrode configurations.IF modes were intended to generate zero induced RF current on the electrodes at a particular distance from the tip, thereby minimizing DBS electrode heating.The proposed technique is based on well-known EM scattering, reciprocity, 37 and superposition 38 theorems.Scattering and reciprocity theorems impose a condition that an induced current on a conductor (in this case, DBS electrode) generates an EM field, which is called a scattered field and is linearly proportional to the induced current's as well as the source's excitation.In addition, the superposition theorem implies that any incident and scattered EM fields caused by an array of sources (in this case, Tx array coil) can be expanded as a linear superposition of incident and scattered fields caused by individual sources.Using these EM theorems in the DBS electrodes problem, the induced current on the electrode by individual pTx channels was related to scattered B1 + fields per individual channel.Then, all linearly independent excitation modes (i.e., IF modes) that generate zero total induced current on the electrode were determined by calculating the null space of the per-channel induced currents.
It is noteworthy that the proposed workflow takes advantage of pTx coils, which are mostly available at ultrahigh fields (defined as ≥ 7 T) and continuously

F I G U R E 4
Amplitudes (blue bars) and phases (red bars) of 14 implant-friendly (IF) excitation vectors corresponding to the 16-channel parallel-transmit coil for the bilateral electrode configuration.The green diamond on each plot corresponds to the predominant channel within that particular IF mode.
growing in number and variety, [39][40][41][42][43][44][45][46][47][48][49] to cancel out the net-induced RF currents on multiple DBS electrodes simultaneously.Theoretically, this approach is not limited to 7 T and can be applied to any pTx coils at any field strength, including lower [39][40][41] and higher [42][43][44][45][46][47][48][49] field strengths.Yet, the implementation of this approach would not be possible using commonly used single-channel or even two-port-driven coils that are commercially available at lower field strengths (i.e., 1.5 T and 3 T).In other words, based on the rank-nullity theorem, 34 to find a nonzero solution for Eq. ( 4) with L number of electrodes and N number of transmit channels, the necessary condition is N > L (i.e., number of IF modes = N -L).Note that because the number of IF modes is always less than the number of independent excitation modes of the pTx coil (i.e., N -L IF modes vs. N independent modes), the proposed technique results in fewer degrees of freedom for shimming compared with non-IF shimming techniques.Although this reduction in degrees of freedom may limit the shimming capability of a pTx coil in realistic imaging scenarios, it can be deemed an acceptable tradeoff to ensure the safety of DBS patients during MRI scans.
In the second step of the workflow, to calculate the scattered B1 + on an axial Plane P 2 retrospectively from the experimentally acquired total B1 + field on the same plane, the incident B1 + fields on Planes P 1 and P 2 were assumed to be approximately equal.This approximation is based on the condition that the distance between planes P 1 and P 2 (∼30 mm) is short enough relative to the EM wavelengths inside the human body at 7 T (∼120 mm).To verify this assumption, relative B1 + maps on two axial planes 30 mm apart were acquired (i.e., Planes P 1 and P 2 ) for 16 single-channel excitations inside the uniform head-shaped phantom without the DBS electrode.As shown in Figure S1, the normalized RMS error of B1 + maps within the indicated ROI is less than 9%.The indicated ROI represents the region around the electrodes' shafts in the bilateral configuration.It should be noted B1 + maps corresponding to per-channel (A) and implant-friendly (IF) mode (B) excitations experimentally acquired on Plane P 2 (see Figure 1) in the unilateral configuration.

F I G U R E 6
B1 + maps corresponding to per-channel (A) and implant-friendly (IF) mode (B) excitations experimentally acquired on Plane P 2 (see Figure 1) in the bilateral configuration.

F I G U R E 7
Experimentally measured temperature changes at the tip of the unilateral electrode due to RF exposures corresponding to the implant-friendly (IF) modes (A), worst-case heating (B), and four arbitrary excitations (C).that this assumption was only verified for the particular pTx coil 31 used in this study.Thus, a similar verification procedure must be performed if there are any changes in the field strength and/or pTx coil.In a similar manner, the validity of this assumption in a human head imaging scenario was evaluated by in vivo B1 + mapping on axial Planes P 1 and P 2 (i.e., 30 mm apart).The corresponding results with a median normalized RMS error of 11% are presented in Figure S2.
In support of the hypothesis that an IF mode should be free of scattered B1 + fields, Figures 5 and 6 show that the total B1 + fields corresponding to IF modes form relatively uniform maps around the electrode, whereas some of the single-channel excitations present excessive anomalies around the electrode (e.g., Channels 1, 2, 3, and 5 in Figure 5 and Channels 2, 3, and 4 in Figure 6).These anomalies are the consequence of complex interferences between the incident fields and scattered fields caused by induced currents on the electrode.Nevertheless, not every single-channel excitation induces a significant current on the electrode-causing a recognizable anomaly around the electrode.Such a single-channel excitation-which interacts with the electrodes negligibly-along with minimal contributions from other channels, is likely to form one of the IF modes.For example, Channel 9 in Figure 5A does not cause an anomaly around the electrode and has almost an identical B1 + profile as IF Mode 8 in Figure 5B.This is also reflected in the excitation vector of IF Mode 8 shown in Figure 3, where channel 9 (the predominant channel indicated by the green diamond) is almost the only contributing channel to the excitation vector.A similar example in the bilateral configuration is Channel 16 in Figure 6A and IF Mode 14 in Figure 6B.The excitation vector corresponding to IF Mode 14, shown in Figure 4, supports the idea that Channel 16, with minimal contributions from other channels, can serve as an IF excitation mode.
The efficacy of the IF mode calculations was further assessed by monitoring the temperature change at the tip of the electrodes while the setup was exposed to RF energy with IF excitation modes.As shown in Figures 7 and 8, no heating at a level measurable with our temperature probes was detected.Furthermore, for both unilateral and bilateral configurations, the worst-case heating scenarios (Figures 7B and 8B) were calculated and used in heating experiments as demonstrations of potential risks of DBS electrodes undergoing MRI scans.In the worst-case scenario, the temperature at the tip of the unilateral electrode Experimentally measured temperature changes at the tips of the bilateral electrodes due to RF exposures corresponding to the implant-friendly (IF) modes (A) and worst-case heating excitations (B) for Electrodes 1 and 2.
was increased up to 3.5 • C for a 6-W RF exposure.Similarly, the worst-case scenario for the bilateral configuration caused 3.5 • C and 2.5 • C temperature increases, respectively, in Electrode 1 and Electrode 2 for a 12-W RF exposure.Because the worst-case heating scenarios were deliberately calculated, the likelihood of occurrence of such excitations in a clinical scan can be argued.To elaborate on such an argument and to further emphasize the significance of the IF excitation strategies, four arbitrary excitations with a 6-W input power were applied to the unilateral configuration, and ΔT from about 1 • C to 2.5 • C were recorded (Figure 7C).To further investigate the importance of IF excitations, some risks and drawbacks of non-IF excitations were evaluated in a bilateral electrode setup (Figure S3).In the corresponding experiment, 3D-TSE imaging was performed using a circularly polarized excitation with a total average power of 12 W.Although the temperature increase at electrodes was negligible, the imaging artifact around electrode 2 (left electrode) was considerably high.Moreover, the temperature at the tip of electrodes corresponding to some of the per-channel excitations was measured.As a result of this measurement, temperature increases were observed in some cases (e.g., Channels 14 and 15), implying a potential risk of non-IF excitations.
To attain a complete picture of the heating propensity of the DBS electrode with the pTx coil under investigation, a Monte-Carlo simulation with 10 7 random excitation vectors was performed, and corresponding ΔT at the tip of Electrode 1 in the bilateral configuration (Figure 1D) was recorded and presented as a histogram in Figure S4.To perform the Monte-Carlo simulation, the induced current on the DBS electrode, I ind exc (see Eq. [4]), corresponding to each random excitation vector, was calculated and normalized to that of the worst-case scenario.The square of the normalized induced current was used to scale the worst-case ΔT (Figure 8B) and calculate the corresponding ΔT.Although the worst-case excitation heats up the electrode by 3.5 • C, the temperature increase remains below 2 • C in 99% of cases, for this particular trajectory.
To prove the concept of IF modes and to demonstrate the use of IF modes in RF shimming scenarios, a simple shimming scenario with the objective of a uniform FA within an ROI was investigated, in which a rectangular ROI in the right anterior of the phantom was chosen, representing an ROI in a brain's frontal lobe.This shimming scenario was carried out using IF modes as well as the scanner's shimming toolbox.Although the former uses the IF excitation modes as the basis of the excitation space, the latter uses individual-channel excitations for the shimming problem.As shown in Figure 9 for the unilateral configuration, the IF mode solution did not result in any measurable heating for a 3D-TSE sequence, whereas the scanner's shim solution caused an approximate 2.5 • C temperature increase for the same pulse sequence parameters.The induced current mitigation also improved image quality in the vicinity of the electrode.In the case of the scanner's shimming toolbox (Figure 9A), the induced current on the electrode generated scattered B1 + field near the electrode.Complex interference of the scattered and incident B1 + fields caused overflips and underflips around the electrode, which resulted in significant inhomogeneity in the TSE image around the electrode.The IF solution (Figure 9C) was free of the induced current, and hence resulted in more homogeneous image intensity.
Similarly, in the case of the bilateral configuration (Figure 10), both the scanner's shimming toolbox and IF modes were used to calculate the homogeneous excitation within the rectangular ROI in the right anterior of the phantom.Using the scanner's shim solution, while Electrode 1 was heated negligibly (∼0.5 • C), Electrode 2 experienced significant heating for a tissue as sensitive as the brain (∼2 • C).It is noteworthy that in this configuration, the IF shim solution did not entirely eliminate the induced current for Electrode 2, and the residual induced current caused minor heating (∼0.5 • C).This can be justified by the fact that induced currents on electrodes in the bilateral configuration interact with each other, which could be considered in Eq. ( 1) by adding a second term as follows: The second term in Eq. ( 8) represents the scattered B1 + field by Electrode 2 around the shaft of Electrode 1.
Because the scattered field decays with 1/r (distance from the electrode), the first term is considerably more dominant around Electrode 1 compared with the second term, and vice versa.Thus, the second term in Eq. ( 8) was neglected for the simplicity of calculation.
In the proposed technique, correctly determining the location of Plane P 2 (at the electrode's shaft) is essential.On the one hand, because susceptibility artifacts are extremely high around the metallic contacts of the electrode and disrupt the B1 + mapping close to the contacts, Plane P 2 should intersect the electrode at a location sufficiently far from the metallic contacts.On the other hand, the assumption of equal-incident B1 + fields on Planes P 1 (underneath the electrode's tip) and P 2 (at the electrode's shaft) only holds if the distance between them is short enough relative to the EM wavelength.Therefore, Plane P 2 should be positioned close to the electrode's tip.These opposing conditions result in a narrow range of possible locations for Plane P 2 .In the case of the commercial electrode used throughout this work, the optimal location on the electrode's shaft for Plane P 2 was determined to be 30 mm away from the electrode's tip.
In this study, experimental validations for IF mode calculation workflow were performed in a uniform phantom.However, heterogeneity of the human head may increase the nonuniformity of B1 + fields (as shown in Figure S2), as well as susceptibility artifacts around the electrode.Both factors can degrade the efficacy of the proposed IF mode calculation workflow.In this regard, further investigations in heterogeneous media, such as animal specimens, are required.

CONCLUSIONS
In this study, a systematic workflow was proposed to determine the implant-friendly RF excitation modes of a pTx coil for unilateral and bilateral electrode configurations undergoing an MRI scan at 7 T.The proposed workflow was experimentally tested and validated by means of B1 + mapping and temperature monitoring around commercial DBS electrodes.None of the calculated IF excitation modes caused a measurable temperature increase at the tip of the electrodes.The IF modes were also used to obtain an RF shim solution within a predetermined ROI.Using clinically relevant 2D-TSE and 3D-TSE pulse sequences, corresponding solutions were compared with those of the scanner's shimming toolbox.Findings suggested that while the outcome of the scanner's shimming toolbox can cause excessive heating in electrodes, the IF shim solution does not heat up the electrode.Hence, the proposed workflow can determine implant-friendly excitations to safely scan DBS patients.

F I G U R E 1
Depiction of the deep brain stimulation (DBS) electrode and its orientation in the experimental setup.(A) The coordinate parameters used in calculations of the B1 + field of the induced current on an oblique DBS electrode.(B) Depiction of the commercial segmented DBS electrode used throughout this study.(C) A sagittal depiction of the DBS electrode's trajectory in the unilateral configuration.(D) A coronal depiction of the DBS electrodes' trajectories in the bilateral configuration.

F I G U R E 9
RF shimming within the indicated region of interest (ROI) using the scanner's shimming toolbox versus the implant-friendly (IF) modes for the unilateral configuration.(A-C) Two-dimensional turbo spin-echo (TSE) images.The red circle shows the area around the electrode, whose coefficient of variation (CoV) was chosen to evaluate image uniformity.(B-D) Experimentally measured temperature change at the tip of the electrode.F I G U R E 10 RF shimming within the indicated region of interest (ROI) using the scanner's shimming toolbox versus the implant-friendly (IF) modes for the bilateral configuration.(A-C) Three-dimensional turbo spin-echo (TSE) images.The red circle shows the area around the electrode, whose coefficient of variation (CoV) was chosen to evaluate image uniformity.(B-D) Experimentally measured temperature change at the tip of the electrodes.

Figure S2 .
Figure S2.Verification of the assumption that the incident B1 + fields on Planes P1 and P2 (see Figure1) around the electrodes are equal with a negligible error.This relative B1 + -mapping experiment was conducted in vivo in a healthy human subject head (without deep brain stimulation [DBS] electrodes).(A) A depiction of Planes P1 and P2 with respect to the head.(B) Experimentally acquired B1 + maps on Planes P1 and P2 corresponding to single-channel excitations.The region of interest (ROI) shown as red rectangular represents the region where the shaft of electrodes would be located in a hypothetical bilateral DBS configuration.The alteration of the incident B1 + fields within the ROI between two planes was quantified using normalized RMS error (NRMSE).(C) A boxplot representation of the NRMSE values throughout different single-channel excitations.Although the maximum NRMSE is 20% (corresponding to Channel 8 excitation), the median of NRMSEs is 11%. Figure S3.Evaluation of temperature increase in circularly polarized (CP) mode and per-channel excitations, as well as the image quality of the CP excitation.In all cases, a 3D turbo spin-echo (TSE) sequence (flip angle [FA] = 100 • , TR/TE = 200/105 ms, in-plane resolution = 0.35 mm, slice thickness = 1.25 mm, acquisition matrix = 576 × 432, TSE factor = 15, no. of averages = 1, pixel bandwidth = 225 Hz/pixel, with a total average power of 12 W) was used for imaging.(A) A depiction of the bilateral deep brain stimulation (DBS) electrode setup inside a uniform phantom with a contralateral configuration.

2
32attered B 1 + fields: Assuming the incident B1 + field on Plane P 2 and the total B1 + field on P 1 are approximately equal,32per-channel complex B1 the induced current by Channel 1, I ind 1 .For this purpose, N−1 minimization routines using MATLAB's fminsearch function were performed on pairs of scattered B1 + fields ( i.e., I ind 1 and I ind + maps on Plane P1 were subtracted from those on Plane P 2 to approximate the scattered B1 + fields, B1 +sca  , on Plane P 2 .3 Relative induced currents: Governed by Eq. (2), the induced current on the electrode's shaft caused by each Tx channel excitation, I ind  , is proportional to the scattered B1 + field on plane P 2 , B1 +sca  , and was calculated relative to