Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

In-vivo and numerical analysis of the eigenmodes produced by a multi-level Tic-Tac-Toe head transmit array for 7 Tesla MRI

  • Tales Santini ,

    Contributed equally to this work with: Tales Santini, Yujuan Zhao

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Bioengineering, University of Pittsburgh, Pittsburgh, PA, United States of America

  • Yujuan Zhao ,

    Contributed equally to this work with: Tales Santini, Yujuan Zhao

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Bioengineering, University of Pittsburgh, Pittsburgh, PA, United States of America

  • Sossena Wood,

    Roles Conceptualization, Data curation, Visualization, Writing – review & editing

    Affiliation Department of Bioengineering, University of Pittsburgh, Pittsburgh, PA, United States of America

  • Narayanan Krishnamurthy,

    Roles Software, Writing – review & editing

    Affiliation Department of Bioengineering, University of Pittsburgh, Pittsburgh, PA, United States of America

  • Junghwan Kim,

    Roles Conceptualization, Writing – review & editing

    Affiliation Department of Bioengineering, University of Pittsburgh, Pittsburgh, PA, United States of America

  • Nadim Farhat,

    Roles Data curation, Writing – review & editing

    Affiliation Department of Bioengineering, University of Pittsburgh, Pittsburgh, PA, United States of America

  • Salem Alkhateeb,

    Roles Data curation, Validation, Writing – review & editing

    Affiliation Department of Bioengineering, University of Pittsburgh, Pittsburgh, PA, United States of America

  • Tiago Martins,

    Roles Software, Writing – review & editing

    Affiliation Department of Bioengineering, University of Pittsburgh, Pittsburgh, PA, United States of America

  • Minseok Koo,

    Roles Data curation, Writing – review & editing

    Affiliation Department of Bioengineering, University of Pittsburgh, Pittsburgh, PA, United States of America

  • Tiejun Zhao,

    Roles Conceptualization, Software, Writing – review & editing

    Affiliation Siemens Medical Solutions, Pittsburgh, PA, United States of America

  • Howard J. Aizenstein,

    Roles Conceptualization, Supervision, Writing – review & editing

    Affiliations Department of Bioengineering, University of Pittsburgh, Pittsburgh, PA, United States of America, Department of Psychiatry, University of Pittsburgh Medical Center, Pittsburgh, PA, United States of America

  • Tamer S. Ibrahim

    Roles Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    tibrahim@pitt.edu

    Affiliations Department of Bioengineering, University of Pittsburgh, Pittsburgh, PA, United States of America, Department of Psychiatry, University of Pittsburgh Medical Center, Pittsburgh, PA, United States of America

Abstract

Radio-frequency (RF) field inhomogeneities and higher levels of specific absorption rate (SAR) still present great challenges in ultrahigh-field (UHF) MRI. In this study, an in-depth analysis of the eigenmodes of a 20-channel transmit Tic-Tac-Toe (TTT) RF array for 7T neuro MRI is presented. The eigenmodes were calculated for five different Z levels (along the static magnetic field direction) of the coil. Four eigenmodes were obtained for each Z level (composed of 4 excitation ports), and they were named based on the characteristics of their field distributions: quadrature, opposite-phase, anti-quadrature, and zero-phase. Corresponding finite-difference time-domain (FDTD) simulations were performed and experimental B1+ field maps were acquired using a homogeneous spherical phantom and human head (in-vivo). The quadrature mode is the most efficient and it excites the central brain regions; the opposite-phase mode excites the brain peripheral regions; anti-quadrature mode excites the head periphery; and the zero-phase mode excites cerebellum and temporal lobes. Using this RF array, up to five eigenmodes (from five different Z levels) can be simultaneously excited. The superposition of these modes has the potential to produce homogeneous excitation with full brain coverage and low levels of SAR at 7T MRI.

Introduction

Ultrahigh-field (UHF) magnetic resonance imaging (MRI) can be exploited for medical research and applications through its higher resolution anatomical imaging, inherent higher contrast, and improved spectroscopy. However, there are technical and physical challenges associated with UHF imaging that have not been completely addressed yet: a) the inhomogeneous distribution of the circularly polarized transmit fields (B1+), responsible for excitation [14]; b) the potentially higher power deposition in the tissues [5, 6]; c) the absence of commercial transmit body coil integrated into the system (commonly seen at lower fields); and d) the difficulty to supervise the local specific absorption rate (SAR) [7].

Several designs of radio-frequency (RF) transmit arrays have been proposed to improve the RF (B1+ and SAR) performance at UHF MRI [811]. A major advantage of these multichannel systems is that the channels of the RF arrays can be manipulated to operate at specific amplitudes and phases, having the potential to optimize a certain characteristic of the RF fields distribution (usually improving B1+ homogeneity and/or efficiency and minimizing SAR.) To determine these operational points, some techniques have been applied, among them the eigenmodes approach; for instance, the two-dimension image uniformity of a spherical phantom was 10% by linearly combining four harmonics modes [12]. Moreover, two time-interleaved acquisitions using different modes have shown improvement in the homogeneity without increasing the time of acquisition [13]. Eigenmode approaches have also been utilized to analyze the signal-to-noise ratio (SNR) behavior of phased array receive coils [14, 15] and to increase the acceleration factor in parallel imaging [16].

In this work, a description and an excitation paradigm are presented for a 20-channel, five-sided Tic-Tac-Toe (TTT) RF transmit array design for 7 Tesla (T) MRI [9]. The RF coil performance (B1+ and SAR) was studied using the eigenmodes approach. The modes were numerically calculated from finite-difference time-domain (FDTD) simulations and experimentally verified in-vivo and on a spherical phantom with a 7T human MRI scanner. Using the designed RF array, up to five eigenmodes can be excited simultaneously. The combination of these eigenmodes has the potential to achieve an efficient and homogeneous B1+ field distribution with low levels of SAR at UHF MRI.

Material and methods

The RF array design and construction

The TTT coil design has been applied to several UHF human MRI applications, including head [9, 1719], breast [20, 21], torso [22], and foot [23, 24]. Fig 1(A) shows the schematic diagram of a four-element 2x2 TTT transmit array design. The coil is composed of eight square-shaped transmission lines electrically connected to each other in a tic-tac-toe fashion. The outer strut was built from 8μm-thick single-sided copper sheets (Polyflon, Germany). The inner rods are composed of solid square-shaped copper (McMaster-Carr, USA) partially inserted into the outer strut, creating a squared shape coaxial transmission line. The dimensions of the outer strut are 228.6 × 228.6 × 19.0 mm3.

thumbnail
Fig 1. Coil schematic diagram, load position and regions of interest.

In (a), the schematic diagram of a four-element, 2x2 Tic-Tac-Toe array design. The copper rods (1 and 2) are partially inside the copper struts (3) providing matching and tuning to the RF coil. In (b), the assembled RF coil system, composed of 5 sets/sides of the 2x2 Tic-Tac-Toe transmit arrays (total of 20 transmit elements.) In (c), FDTD spherical (~17cm in diameter) water phantom model (108 by 108 by 108 Yee cells with isotropic resolution of ~1.6mm). The red dots indicate the excitation points of the three visible sets of the 2x2 Tic-Tac-Toe arrays. The 5 levels of the coil in Z direction are shown. In (d), the Duke Virtual Family Adult Head Model (114 by 117 by 144 Yee cells with isotropic resolution of ~1.6mm). In (e), the head model was divided into 8 different regions of interest (ROI) as indicated by the color code and the numbers.

https://doi.org/10.1371/journal.pone.0206127.g001

The excitation ports of one side (four channels) are also shown in Fig 1(A). Tuning and matching of the coil is performed by changing the length of the inner rods inside the outer struts, presenting similar performance in terms of s-parameters as the demonstrated in [24]. The RF copper shielding is located at the back of the coil struts (with a gap of 15.8mm) and it functions as the ground of a cavity resonator, being responsible for both increasing the RF efficiency and preventing RF leaking. The RF copper shielding is composed of double layer 4μm thick copper sheets (Polyflon, Germany) and it was slotted with specific patterns to reduce eddy currents while the RF performance is maintained, as demonstrated in [25]. The non-metal parts of the array were 3D printed using polycarbonate (Stratasys, USA).

Fig 1(B) shows the assembled RF coil system which is composed of five sides of the four-element 2x2 Tic-Tac-Toe transmit array (earlier described), resulting in a total of twenty transmit channels/excitation ports. The channels of the RF array were tuned and matched on the bench using the Agilent Network Analyzer Model E5062A (Santa Clara, US). While the five sides of the four-element 2x2 TTT transmit array are inherently decoupled from each other (less than -16dB), on any 2x2 side, the coupling among the adjacent transmit channels (S12 and S14) is about -9 to -11 dB, and the coupling between opposite elements (S13) is about -3 to -4 dB.

FDTD simulations

An in-house FDTD software package with an accurate transmission-line feed mechanism [26] was implemented to model the RF performance of the 20-channel TTT transmit array. This simulation package has been previously utilized and verified [20, 21, 2632]. The RF fields were calculated with the coil loaded with a homogeneous spherical phantom model (diameter = 17.1cm, conductivity = 0.46 S/m and relative permittivity = 79.0) and a human head model (18.2cm × 18.7cm × 23.0cm), which was extracted from the Virtual Family Duke Model [33]. The resolution of the models is ~1.6mm isotropic per voxel and the simulation was run until a steady state was achieved (100,000 time steps, time resolution of 3 ps). Fig 1C and 1D shows respectively the relative position of the homogenous spherical phantom and the human head model inside the RF array. The excitation points are identified by the red dots.

The B1+ field distribution was analyzed in eight regions of interest (ROI) described in Fig 1(E). The ROIs are based on human head anatomical characteristics as well as the electromagnetic characteristics of the coil.

Calculations of the eigenmodes

The current distributions induced on the RF coil elements can be controlled by manipulating the amplitude and phase of the voltages feeding the excitation ports. A specific current distribution induced on the elements of the RF coil also determines an eigenmode [34]. Consequently, the B1+ field distribution and SAR can be manipulated as a result of the superposition of fields produced by the individual elements. In this work, the set of B1+ field distributions was arranged by: (1) where C is the B1+ field matrix generated by an array with L transmit channels, n is the number of Yee cells inside the ROI. C * C gives the correlation among the channels of the array; therefore, the eigenmodes can be calculated by: (2) where v is a unitary matrix of eigenvectors; is a diagonal matrix of eigenvalues. With solutions for Eq 2, Cv is the spatially pseudo-independent fields or eigenmodes of the transmit coil; v gives the phase and amplitude of each coil channel; λi represents the field energy for eigenmode i.

The transmit array was grouped into five levels of four elements along the static magnetic field (Z) direction: Top_Level, Level_1, Level_2, Level_3, and Level_4 (see Fig 1C and 1D). The eigenmodes were calculated in each Z level of the transmit array by applying Eq 2 on the simulated B1+ field distributions; thus totaling four different excitation field patterns per level and 20 in total. Since the magnetic field distribution and SAR are two major concerns for 7T imaging, the attributes of the modes and coil Z levels were evaluated using three criteria:

  1. average B1+ intensities inside each ROI for each Z level and mode, scaled for 1W input power per channel (totaling 4W for one Z level);
  2. B1+ homogeneity calculated by the coefficient of variation (CV) inside each ROI for each Z level and mode;
  3. average and peak SAR over the whole head volume (from the top of the neck) for each Z level and mode, scaled for 1W input power per channel (totaling 4 W for one Z level).

Please note that IEC/FDA limits the SAR in 3.2 W/kg for 10g of tissue inside the human head [35]. SAR levels were therefore evaluated in terms of average SAR over the whole head volume, peak SAR over any 10g of tissues, and safety excitation efficiency (SEE) [36], defined as average B1+ intensity over the combined volume of all eight ROIs divided by the average SAR over the whole head volume .

The eigenmodes were combined using an optimization of the 20-channel B1+ fields. The optimization aims at minimizing the CV of the B1+ field distribution within the ROI that encapsulates the whole head above and including the cerebellum and excluding the nasal cavities. The resultant field distribution was then scaled by 1W of total input power and the SEE was calculated based on the average B1+ field in the ROI divided by the square root of the average SAR for the whole head.

MRI experiments

The FDTD calculated eigenmodes were experimentally verified using the constructed 20-channel transmit array. The MR experiments were conducted using a 7 Tesla MRI scanner (Siemens MAGNETOM, Germany). This study was approved by the University of Pittsburgh’s Institutional Review Board (IRB PRO17030036). One healthy volunteer was scanned after signing a written informed consent. The phantom MRI imaging experiment and the in-vivo study were conducted by acquiring relative B1+ maps using Turbo Flash MRI sequence; the outputs of this MRI sequence are: 1) the B1+ distribution for each transmit channel (scaled to the square root of the sum of the square of all connected transmitting channels); and 2) the relative phases. The sequence parameters used were: TE/TR = 2.34/160ms, resolution 3.2mm isotropic, flip angle 12 degrees. The scanner is equipped with 8 channels in the parallel transmit (pTx) mode with 1kW power amplifier per channel (8kW in total). These 8 transmit-channels were connected to the RF array in 2 Z Levels (each level has 4 channels) in each B1+ mapping experiment. Level_1 (most homogeneous level) was always connected in addition to another level (Fig 1D) per one B1+ mapping measurement. The 4 transmit-channels not connected to Level_1 were manually changed to another level until all the B1+ maps were acquired for all (5) Z levels. A transmit/receive (T/R) switch box was used to receive the signal from all 20 channels for any B1+ mapping acquisition. The transmit channels of the coil that were not used in a specific B1+ mapping acquisition were terminated with 50Ω loads through the T/R box.

Results

Calculation of the eigenmodes

By applying Eq 2 on the FDTD-simulated B1+ fields, the phases and amplitudes of the eigenmodes were obtained for each Z level of the transmit array; the results are presented in Table 1. Four modes were identified, and these modes presented uniformly distributed relative phase shifts and constant amplitudes among the 4 channels of each Z level: Mode_1 (named as quadrature) presents phase increments of ~90°; Mode_2 (opposite-phase) has increments of ~180°; Mode_3 (anti-quadrature) presents increments of ~270°; and Mode_4 (zero-phase) has increments of ~0° or ~360°. There were minor phase (<8°) and amplitude (<8%) deviations among the 5 different Z levels, which is impacted by the position of the load inside the RF coil. For practical purposes, the phases were kept as multiples of 90° and the amplitudes were considered to be the same for all channels.

thumbnail
Table 1. FDTD-calculated relative phases and amplitudes associated with the Eigenmodes of the array’s five Z levels.

The coil was loaded with the homogeneous spherical phantom.

https://doi.org/10.1371/journal.pone.0206127.t001

B1+ field and SAR comparisons of the eigenmodes

B1+ field intensities and homogeneity of the eigenmodes.

The FDTD-calculated B1+ field distribution of all modes for all the Z levels are presented in Fig 2(A). When comparing the eigenmodes in different Z Levels, the following observations are noted:

  1. Mode_1 (quadrature) generally provides high B1+ intensity in the central regions of the head/brain with the bright spot generally moving along the Z direction for distinctive Z levels;
  2. Mode_2 (opposite-phase) generates peripheral brain excitation;
  3. Mode_3 (anti-quadrature) generally excites the periphery of the head;
  4. Mode_4 (zero-phase) excites the lower brain (cerebellum and temporal lobes).

The B1+ field phase distribution maps are shown in Fig 2(B) (note that −2π = 2π, i.e., the intense blue color is equal to the intense red color in the colorbar).

thumbnail
Fig 2. Simulated B1+ field and SAR distributions of the Eigenmodes in the Duke head model for each level (shown in Fig 1).

The central slices in sagittal, coronal, and axial planes are shown. In (a), the amplitude of B1+ field distributions, in μT for 1W input power per channel (total 4W as each level contains 4 channels). For the four Eigenmodes per level, the colorbar is scaled from 0 to the maximum. In (b), the phase of the B1+ field distribution in radians. In (c), the SAR distributions in W/kg for 10g of tissues per 1W input power per channel (total 4W). The coil was loaded with the Duke Virtual Family Adult Head Model.

https://doi.org/10.1371/journal.pone.0206127.g002

The values of B1+ field intensities for all modes, levels, and ROIs are presented in Fig 3, from which we can note that:

  1. Top_level produces an efficient excitation in the upper head (ROIs 1, 2, 5, 6, and 8) when operating in Mode_1, presenting an average B1+ intensity of 0.73μT for 1W input power per channel (total 4 W) in these ROIs.;
  2. Level_1 and Level_4 are also efficient operating in Mode_1, producing an average B1+ of 0.54μT in the ROIs 1, 2, 3, 5, 6, and 8;
  3. Levels 2 and 3 produces an efficient excitation in the lower brain (ROIs 3, 4, and 7) when operating in Mode_4, presenting an average B1+ of 0.48μT in these regions.

The CV of the B1+ field intensities over the specified eight ROIs for different modes and levels are shown in Fig 4.

thumbnail
Fig 3. Average B1+ intensities calculated inside the 8 different regions of interest (ROIs) shown in Fig 1(E) for each Z level of the RF array shown in Fig 1(D).

The scale is in μT for 1W input power per channel (total 4W).

https://doi.org/10.1371/journal.pone.0206127.g003

thumbnail
Fig 4. Coefficient of variation (standard deviation over the mean of B1+ field distribution) calculated inside the 8 different regions of interest (ROIs) shown in Fig 1(E) for each Z level of the RF array shown in Fig 1(D).

https://doi.org/10.1371/journal.pone.0206127.g004

SAR comparison for the eigenmodes at different Z levels.

The numerically calculated SAR distributions for all eigenmodes from all Z levels are shown in Fig 2(C). Preferable modes present higher average B1+ intensity and lower peak and average SAR. The following observations are noted:

  1. the SAR distribution significantly varies for different eigenmodes and Z levels;
  2. the highest SAR regions usually correspond to lower intensities of B1+;
  3. Top_level operating in Mode_1 produces the highest peak SAR, but it is also B1+ efficient;
  4. Levels 1 and 4 produce homogeneous SAR distribution when operating in Mode_1;
  5. Mode_4 produces higher levels of SAR in the lower brain regions (except Top_level);
  6. Mode_1 usually produces low levels of average and peak SAR (except in Top_level) and high levels of SEE.

Fig 5 shows the average/peak SAR and SEE values for all Z levels and eigenmodes of the transmit array.

thumbnail
Fig 5. SAR evaluation of the Eigenmodes for each Z level of the RF array shown in Fig 1(D).

In a) the average SAR per 1 W input power per channel (total 4 W). In b) the peak SAR per 1 W input power per channel (total 4 W). In c), the safety excitation efficiency (SEE) (the B1+ field is averaged over a volume that encapsulates all eight regions of interest.) The results are presented for the Duke Virtual Family Adult Head Model.

https://doi.org/10.1371/journal.pone.0206127.g005

Experimental verification

Fig 6 shows the simulated and measured B1+ maps for the four eigenmodes excited by each of the five Z levels of the 20-channel transmit array. Fig 6(A) and 6(C) show, respectively, the simulated and measured data in the homogeneous spherical phantom. Fig 6(B) shows the simulated B1+ maps in the Duke head model. For a visualization resembling the in-vivo acquired data, a limited number of tissues are shown: tissues distant from the brain (e.g., tongue muscle) or tissues which produce low MR signal (e.g., bone) was removed from the Fig 6(B), although the simulations were conducted using the complete Duke head model. Fig 6(D) shows the in-vivo acquired eigenmodes. The results show excellent agreement between the simulated and measured data.

thumbnail
Fig 6. Experimental verification and simulated B1+ field distributions of the Eigenmodes for the homogenous spherical phantom and the human head, showing the central sagittal view.

In (a), the simulated B1+ field distributions in the homogeneous spherical phantom with relative permittivity of 79 and conductivity 0.41 S/m. In (b), the simulations in the Duke Virtual Family Adult Head Model. In (c), B1+ maps acquired in the homogeneous phantom with relative permittivity of 79 and conductivity 0.41 S/m. In (d), in-vivo human B1 maps. All maps are scaled to the square root of the sum of the square of all connected transmitting channels.

https://doi.org/10.1371/journal.pone.0206127.g006

Combination of the eigenmodes

Fig 7 shows the combination of the modes by minimizing the CV of the B1+ fields in the ROI. The values in the ROI (composed by the whole head above and including the cerebellum and excluding the nasal cavities) are: = 16.6%, = 3.51, SEE = 1.48 (defined as mean B1+ in the ROI divided by the square root of the SAR for the whole head), mean B1+ = 0.23μT for 1W total input power.

thumbnail
Fig 7. An example of the combination of the modes of the Tic-Tac-Toe coil (20 Tx channels).

The ROI represents the entire head above and including the cerebellum and excluding the nasal cavities.

https://doi.org/10.1371/journal.pone.0206127.g007

Discussion

In UHF MRI, as the wavelength of the electromagnetic waves inside the tissues gets closer, in size, to the body parts being scanned, inhomogeneities become a major issue, as it can affect the image quality, creating voids and low contrast regions (especially in high flip-angle sequences). In the case of brain imaging, this situation is usually accentuated in the lower brain regions such as cerebellum and temporal lobes [37]. There are several works suggesting the use of two modes to increase the homogeneity of the B1+ field distribution [13, 38, 39]. Another work suggested that the coefficient of variation of a 2D image can reach 10% by using four eigenmodes in a homogeneous spherical phantom using a birdcage RF coil [12]. At many instances, the application of these methods can come at a significant elevation of time of acquisition, elevated SAR, and difficulties in simultaneously exciting several distinct modes of a coil [34, 40, 41].

The freedom to manipulate the current distribution of different coil elements potentially contributes to the generation of a homogeneous B1+ fields distribution [4244]. However, coil arrays typically show the capability to control current distributions only at XY plane, while current distribution are not very commonly controlled in the Z direction. It is worth noting that there are some coil designs that can potentially generate current control along the Z direction. Some examples are: the multi-rows/rings transmit arrays that allow parallel transmission approaches [4547]; the rotating RF coil approach [48, 49]; and the spiral volume coil [50].

In this work, the eigenmodes of a 20-channel Tic-Tac-Toe RF array were studied. The RF array is composed of five excitation levels located at different positions along the static magnetic field. For each level (composed of four ports located in the same XY plane), there are four distinctive modes (with 90° phase-shift multiples) that can be generated, calculated using Eq 2. Using power splitters and phase shifters (1-to-N network), up to 5 different modes can be excited simultaneously in a single image acquisition (since each Z level can present a different excitation mode), potentially improving B1+ homogeneity and reducing SAR levels. Thus, 1024 possible combinations can be implemented using the five Z levels of the RF array, if the amplitudes of the channels are fixed.

It is important to analyze the field distribution of the eigenmodes provided by an RF head array so that a target homogeneous/low SAR excitation can be achieved. In terms of B1+ distribution inside the human head, our results show that: Mode_1 (quadrature) is the most efficient, producing center brightness at different Z levels. However, voids are observed in some regions in the lower brain (such as the cerebellum and temporal lobe regions). Mode_2 (Opposite-phase) produces low signal in the center and excites mostly the periphery regions of the brain. Mode_3 excites regions in the head periphery (mostly skin and skull) and can have localized functions such as fat suppression (extracerebral lipids from skins and skull can be suppressed to reduce the influence from this region and leave the central brain regions unaffected). Mode_4 excites mostly the lower brain regions (cerebellum and temporal lobes). The analysis also shows that Mode_1 of Level_1 and Level_4 can excite relatively uniform B1+ distributions, with CV = 22% and 22% inside ROI8 (upper head) for Level_1 and Level_4 respectively.

While there can be many solutions for the RF excitation that achieve a satisfactory signal fidelity to the targeted excitation pattern (e.g., homogeneous B1+ field), minimizing the local SAR is also an important target for the coil design and operation. In this work, the average and peak SAR was compared for different Z levels and eigenmodes. It is important to note that the SAR distribution presented in this work is an outcome of the phases and amplitudes determined by the eigenmodes, which were calculated using only the B1+ fields. Therefore, lower levels of SAR can be achieved if SAR constraints are included. Mode_1 produces lower average (< 0.11 W/kg) and peak SAR (except Top_Level), combined with efficient B1+ in the upper head, leading to a high SEE (>1.5 ) as seen in the Levels 1, 2 and 4. This is higher than the double row loop array (0.76 [36, 51]) and the birdcage coil (0.89 [52]) for instance. While Top_level produces a high peak SAR, it is B1+ efficient, resulting in SEE of ~1.5. Mode_2 (opposite-phase) produces a relatively high SEE for Levels 2 and 3, with relatively high brain peripheral excitation. Although Mode_4 (zero-phase) presents low levels of SEE, it has efficient B1+ in the low brain regions which are challenging at UHF MRI.

The simulations were experimentally verified by acquiring the individual channels B1+ maps in the homogeneous spherical phantom and in-vivo in human subjects. The field distributions of the eigenmodes were then calculated and compared with the simulated fields (Fig 6). The modes are highly consistent between simulations and experiments. Small differences may be due to differences in the head/phantom position in simulations and experiments, differences in the tuning of the RF coil elements, differences in the hardware of the transmitting channels, and differences in the human head model and the subject scanned. Discrepancies in the B1+ maps between the phantom and the head can be mostly attributed to dielectric effects–that commonly occur in homogeneous water phantoms [53]–and to the anatomical differences between the two models.

An example of the combination (RF shimming) of the modes demonstrates a high level of the homogeneity and coverage of the B1+ field over the ROI, as demonstrated by the values of = 16.6%, and = 3.51. The low level of SAR is also demonstrated with a high level of SEE (1.48 ) even though SAR constraints were not included as a part of the RF shimming. The strong coupling between opposite channels (-3 to -4dB) can improve the load insensitivity of the array (being able to scan subjects with different head volumes/shapes and achieve similar RF characteristics), with the cost of lower transmit efficiency. Nevertheless, an example of the combination of the modes (Fig 7) shows that the transmit RF array produces enough B1+ intensity to perform inversion with a 1ms square RF pulse using 8kW power amplifier capability with ~35% loss to the coil port.

Conclusions

The eigenmode arrangement of the TTT 20-channel RF array potentially allows controlling RF excitation not only at XY plane but also along the Z direction. As five eigenmodes from different Z levels can be excited simultaneously (one per excitation level in Z), we believe that the combination of these modes can provide a full brain homogeneous B1+ excitation. Future work will include the combination/superposition [6, 5457] of these eigenmodes in order to obtain a homogeneous and efficient B1+ field distribution with low levels of SAR.

References

  1. 1. Tropp J. Image brightening in samples of high dielectric constant. J Magn Reson. 2004;167(1):12–24. pmid:14987593.
  2. 2. Ibrahim TS, Mitchell C, Abraham R, Schmalbrock P. In-depth study of the electromagnetics of ultrahigh-field MRI. Nmr Biomed. 2007;20(1):58–68. Epub 2006/09/29. pmid:17006885.
  3. 3. Vaughan JT, Garwood M, Collins CM, Liu W, DelaBarre L, Adriany G, et al. 7T vs. 4T: RF power, homogeneity, and signal-to-noise comparison in head images. Magn Reson Med. 2001;46(1):24–30. pmid:11443707.
  4. 4. Ibrahim TS, Tang L. Insight into RF power requirements and B1 field homogeneity for human MRI via rigorous FDTD approach. J Magn Reson Imaging. 2007;25(6):1235–47. Epub 2007/05/24. pmid:17520721.
  5. 5. Röschmann P. Radiofrequency penetration and absorption in the human body: Limitations to high‐field whole‐body nuclear magnetic resonance imaging. Medical physics. 1987;14(6):922–31. pmid:3696080
  6. 6. Vaughan JT, Snyder CJ, DelaBarre LJ, Bolan PJ, Tian J, Bolinger L, et al. Whole‐body imaging at 7T: preliminary results. Magnetic resonance in Medicine. 2009;61(1):244–8. pmid:19097214
  7. 7. Eichfelder G, Gebhardt M. Local specific absorption rate control for parallel transmission by virtual observation points. Magnetic resonance in medicine. 2011;66(5):1468–76. pmid:21604294.
  8. 8. Lattanzi R, Sodickson DK. Ideal current patterns yielding optimal signal‐to‐noise ratio and specific absorption rate in magnetic resonance imaging: Computational methods and physical insights. Magnetic Resonance in Medicine. 2012;68(1):286–304. pmid:22127735
  9. 9. Ibrahim T, Zhao Y, Krishnamurthy N, Raval S, Zhao T, Wood S, et al. 20-To-8 Channel Tx Array with 32-Channel Adjustable Receive-Only Insert for 7T Head Imaging. ISMRM; 2013.
  10. 10. Aussenhofer S, Webb A. An eight-channel transmit/receive array of TE 01 mode high permittivity ceramic resonators for human imaging at 7T. Journal of Magnetic Resonance. 2014;243:122–9. pmid:24818565
  11. 11. Avdievich NI, Giapitzakis IA, Pfrommer A, Henning A. Decoupling of a tight-fit transceiver phased array for human brain imaging at 9.4T: Loop overlapping rediscovered. Magnetic resonance in medicine. 2017. pmid:28603846.
  12. 12. Taracila V, Petropoulos LS, Eagan TP, Brown RW. Image uniformity improvement for birdcage‐like volume coils at 400 MHz using multichannel excitations. Concepts in Magnetic Resonance Part B: Magnetic Resonance Engineering. 2006;29(3):153–60.
  13. 13. Orzada S, Maderwald S, Poser BA, Bitz AK, Quick HH, Ladd ME. RF Excitation Using Time Interleaved Acquisition of Modes (TIAMO) to Address B(1) Inhomogeneity in High-Field MRI. Magnetic Resonance in Medicine. 2010;64(2):327–33. WOS:000280422000002. pmid:20574991
  14. 14. King SB, Varosi SM, Duensing GR. Eigenmode analysis for understanding phased array coils and their limits. Concept Magn Reson B. 2006;29B(1):42–9. ISI:000235523200005.
  15. 15. King SB, Varosi SM, Duensing GR. Optimum SNR Data Compression in Hardware Using an Eigencoil Array. Magnet Reson Med. 2010;63(5):1346–56. WOS:000277098100022. pmid:20432305
  16. 16. Wang C, Qu P, Shen GX. Potential advantage of higher-order modes of birdcage coil for parallel imaging. J Magn Reson. 2006;182(1):160–7. pmid:16829146.
  17. 17. Ibrahim TS, Hue Y-K, Boada FE, Gilbert R. Tic Tac Toe: Highly-Coupled, Load Insensitive Tx/Rx Array and a Quadrature Coil Without Lumped Capacitors. Intl Soc Mag Reson Med; 2008.
  18. 18. Ibrahim T ST, Raval S, Krishnamurthy N, Wood S, Kim J, Zhao Y, Wu X, Yacoub E, Aizenstein H, Zhao T. Towards Homogenous 7T Neuro Imaging: Findings and Comparisons between 7T TTT and NOVA RF Coil Systems. ISMRM; 2017.
  19. 19. Zhao Y, Zhao T, Krishnamurthy N, Ibrahim T. On the E-field construction/deconstruction and B1+ Efficiency/Homogeneity with Transmit Array Eigen Modes. ISMRM; 2014; Milan, Italy.
  20. 20. Kim J, Santini T, Bae KT, Krishnamurthy N, Zhao Y, Zhao T, et al. Development of a 7 T RF coil system for breast imaging. NMR in biomedicine. 2016. pmid:27859861.
  21. 21. Kim J, Krishnamurthy N, Santini T, Zhao Y, Zhao T, Bae KT, et al. Experimental and numerical analysis of B1+ field and SAR with a new transmit array design for 7T breast MRI. J Magn Reson. 2016;269:55–64. pmid:27240143.
  22. 22. Raval S, Santini T, Wood S, Krishnamurthy N, Ibrahim, TS. In-vivo (8x4) 32-ch Tx-only Body Array for UHF MR.ISMRM; 2017.
  23. 23. Santini T, Kim J, Wood S, Krishnamurthy N, Raval S, Ibrahim T. A new RF coil for foot and ankle imaging at 7T MRI. ISMRM; 2017.
  24. 24. Santini T, Kim J, Wood S, Krishnamurthy N, Farhat N, Maciel C, et al. A new RF transmit coil for foot and ankle imaging at 7T MRI. Magn Reson Imaging. 2018;45:1–6. Epub 2017/09/13. pmid:28893660.
  25. 25. Zhao Y, Zhao T, Raval SB, Krishnamurthy N, Zheng H, Harris CT, et al. Dual optimization method of radiofrequency and quasistatic field simulations for reduction of eddy currents generated on 7T radiofrequency coil shielding. Magnetic resonance in medicine. 2014. pmid:25367703.
  26. 26. Ibrahim TS, Hue YK, Tang L. Understanding and manipulating the RF fields at high field MRI. Nmr Biomed. 2009;22(9):927–36. Epub 2009/07/22. pmid:19621335.
  27. 27. Hue Y-K., Ibrahim TS, Zhao T, Boada FE, Qian Y. A Complete Modeling System with Experimental Validation for Calculating the Transmit and Receive Fields, Total Power Deposition, Input Impedance, and Coupling between Coil Elements. ISMRM 2008. p. 1193.
  28. 28. Tang L, Hue YK, Ibrahim TS. Studies of RF Shimming Techniques with Minimization of RF Power Deposition and Their Associated Temperature Changes. Concepts Magn Reson Part B Magn Reson Eng. 2011;39B(1):11–25. Epub 2011/05/25. pmid:21607117; PubMed Central PMCID: PMC3098508.
  29. 29. Zhao Y, Tang L, Rennaker R, Hutchens C, Ibrahim TS. Studies in RF Power Communication, SAR, and Temperature Elevation in Wireless Implantable Neural Interfaces. PloS one. 2013;8(11):e77759. pmid:24223123; PubMed Central PMCID: PMC3819346.
  30. 30. Zhao Y, Rennaker RL, Hutchens C, Ibrahim TS. Implanted miniaturized antenna for brain computer interface applications: analysis and design. PloS one. 2014;9(7):e103945. pmid:25079941; PubMed Central PMCID: PMC4117534.
  31. 31. Raval SB, Zhao T, Krishnamurthy N, Santini T, Britton C, Gorantla VS, et al. Ultra-high-field RF coil development for evaluating upper extremity imaging applications. NMR in biomedicine. 2016;29(12):1768–79. pmid:27809383.
  32. 32. Raval S, Santini T, Wood S, Krishnamurthy N, Ibrahim TS. In-vivo (8x4) 32-ch Tx-only Body Array for UHF MR. ISMRM; 2017.
  33. 33. Christ A, Kainz W, Hahn EG, Honegger K, Zefferer M, Neufeld E, et al. The Virtual Family—development of surface-based anatomical models of two adults and two children for dosimetric simulations. Physics in medicine and biology. 2010;55(2):N23–38. pmid:20019402.
  34. 34. Vester M, Nistler J, Oppelt R, Renz W. Using a mode concept to reduce hardware needs for multichannel transmit array. ISMRM; 2006.
  35. 35. Fiedler TM, Ladd ME, Bitz AK. SAR Simulations & Safety. Neuroimage. 2017. pmid:28336426.
  36. 36. Kozlov M, Turner R. Analysis of RF transmit performance for a 7T dual row multichannel MRI loop array. Annual International Conference of the IEEE Engineering in Medicine and Biology Society IEEE Engineering in Medicine and Biology Society Annual Conference. 2011;2011:547–53. Epub 2012/01/19. pmid:22254369.
  37. 37. Kraff O, Quick HH. 7T: Physics, safety, and potential clinical applications. J Magn Reson Imaging. 2017. pmid:28370675.
  38. 38. Wiggins OK G. C., Zakszewski E.,Alagappan V.,Wiggins C. J. and Wald L. L. A 7 Tesla Gradient Mode Birdcage Coil for Improved Temporal and Occipital Lobe SNR. ISMRM; Seattle, Washington, USA2006.
  39. 39. Orzada S, Maderwald S, Poser BA, Johst S, Kannengiesser S, Ladd ME, et al. Time-interleaved acquisition of modes: an analysis of SAR and image contrast implications. Magnetic resonance in medicine. 2012;67(4):1033–41. pmid:21858867.
  40. 40. Yazdanbakhsh P, Fester M, Oppelt R, Bitz A, Kraff O, Orzada S, et al. Variable Power Combiner for a 7T Butler Matrix Coil Array. ISMRM. 2009:397.
  41. 41. Alagappan V, Setsompop K, Nistler J, Potthast A, Schmitt F, Adalsteinsson E, et al. A Simplified 16-channel Butler Matrix for Parallel Excitation with the Birdcage Modes at 7T. Proc Intl Soc Mag Reson Med 16; 2008.
  42. 42. Pauly J, Nishimura D, Macovski A. A k-space analysis of small-tip-angle excitation. Journal of magnetic resonance. 2011;213(2):544–57. pmid:22152370
  43. 43. Setsompop K, Alagappan V, Gagoski B, Witzel T, Polimeni J, Potthast A, et al. Slice-selective RF pulses for in vivo B1+ inhomogeneity mitigation at 7 tesla using parallel RF excitation with a 16-element coil. Magnetic resonance in medicine. 2008;60(6):1422–32. pmid:19025908; PubMed Central PMCID: PMC2635025.
  44. 44. Katscher U, Bornert P, Leussler C, van den Brink JS. Transmit SENSE. Magnetic resonance in medicine. 2003;49(1):144–50. Epub 2003/01/02. pmid:12509830.
  45. 45. Gilbert KM, Curtis AT, Gati JS, Klassen LM, Menon RS. A radiofrequency coil to facilitate B-1(+) shimming and parallel imaging acceleration in three dimensions at 7 T. NMR in biomedicine. 2011;24(7):815–23. WOS:000294686900008. pmid:21834005
  46. 46. Shajan G, Kozlov M, Hoffmann J, Turner R, Scheffler K, Pohmann R. A 16-channel dual-row transmit array in combination with a 31-element receive array for human brain imaging at 9.4 T. Magnetic resonance in medicine. 2013. pmid:23483645.
  47. 47. Avdievich NI. Transceiver-Phased Arrays for Human Brain Studies at 7 T. Applied Magnetic Resonance. 2011;41(2–4):483–506. WOS:000299295600030. pmid:23516332
  48. 48. Trakic A, Li BK, Weber E, Wang H, Wilson S, Crozier S. A Rapidly Rotating RF Coil for MRI. Concept Magn Reson B. 2009;35B(2):59–66. WOS:000265720100001.
  49. 49. Trakic A, Weber E, Li BK, Wang H, Liu F, Engstrom C, et al. Electromechanical Design and Construction of a Rotating Radio-Frequency Coil System for Applications in Magnetic Resonance. Ieee Transactions on Biomedical Engineering. 2012;59(4):1068–75. WOS:000302175300019. pmid:22231668
  50. 50. Alsop DC, Connick TJ, Mizsei G. Spiral volume coil for improved RF field homogeneity at high static magnetic field strength. Magnetic Resonance in Medicine. 1998;40(1):49–54. ISI:000074302000006. pmid:9660552
  51. 51. Kozlov M., Möller HE. Safety excitation efficiency of MRI 300MHz dualrow transmit arrays. Antennas and Propagation & USNC/URSI National Radio Science Meeting; 2015; Vancouver: IEEE.
  52. 52. Collins CM, Li S, Smith MB. SAR and B1 field distributions in a heterogeneous human head model within a birdcage coil. Magnetic Resonance in Medicine. 1998;40(6):847–56. pmid:9840829
  53. 53. Hoult DI. Sensitivity and Power Deposition in a High-Field Imaging Experiment. J Magn Reson Imaging. 2000;12(1):46–67. pmid:10931564
  54. 54. Ibrahim TS, Lee R, Baertlein BA, Abduljalil AM, Zhu H, Robitaille PML. Effect of RF coil excitation on field inhomogeneity at ultra high fields: A field optimized TEM resonator. Magnetic resonance imaging. 2001;19(10):1339–47. WOS:000173435100012. pmid:11804762
  55. 55. Tang L, Hue YK, Ibrahim TS. Studies of RF Shimming Techniques with Minimization of RF Power Deposition and Their Associated Temperature Changes. Concept Magn Reson B. 2011;39B(1):11–25. WOS:000290679600002. pmid:21607117
  56. 56. van den Bergen B, van den Berg CAT, Bartels LW, Lagendijk JJW. 7 T body MRI: B-1 shimming with simultaneous SAR reduction. Physics in medicine and biology. 2007;52(17):5429–41. WOS:000249089900022. pmid:17762096
  57. 57. Van den Berg CAT, Van den Bergen B, de Kamer JBV, Raaymakers BW, Kroeze H, Bartels LW, et al. Simultaneous B-1(+) homogenization and specific absorption rate hotspot suppression using a magnetic resonance phased array transmit coil. Magnetic Resonance in Medicine. 2007;57(3):577–86. WOS:000244657200015. pmid:17326185