Next-generation MRI scanner designed for ultra-high-resolution human brain imaging at 7 Tesla.

To increase granularity in human neuroimaging science, we designed and built a next-generation 7 Tesla magnetic resonance imaging scanner to reach ultra-high resolution by implementing several advances in hardware. To improve spatial encoding and increase the image signal-to-noise ratio, we developed a head-only asymmetric gradient coil (200 mT m-1, 900 T m-1s-1) with an additional third layer of windings. We integrated a 128-channel receiver system with 64- and 96-channel receiver coil arrays to boost signal in the cerebral cortex while reducing g-factor noise to enable higher accelerations. A 16-channel transmit system reduced power deposition and improved image uniformity. The scanner routinely performs functional imaging studies at 0.35-0.45 mm isotropic spatial resolution to reveal cortical layer functional activity, achieves high angular resolution in diffusion imaging and reduces acquisition time for both functional and structural imaging.

Experimental validation of a PNS-optimized whole-body gradient coil.

PURPOSE: Peripheral nerve stimulation (PNS) limits the usability of state-of-the-art whole-body and head-only MRI gradient coils. We used detailed electromagnetic and neurodynamic modeling to set an explicit PNS constraint during the design of a whole-body gradient coil and constructed it to compare the predicted and experimentally measured PNS thresholds to those of a matched design without PNS constraints. METHODS: We designed, constructed, and tested two actively shielded whole-body Y-axis gradient coil winding patterns: YG1 is a conventional symmetric design without PNS-optimization, whereas YG2's design used an additional constraint on the allowable PNS threshold in the head-imaging landmark, yielding an asymmetric winding pattern. We measured PNS thresholds in 18 healthy subjects at five landmark positions (head, cardiac, abdominal, pelvic, and knee). RESULTS: The PNS-optimized design YG2 achieved 46% higher average experimental thresholds for a head-imaging landmark than YG1 while incurring a 15% inductance penalty. For cardiac, pelvic, and knee imaging landmarks, the PNS thresholds increased between +22% and +35%. For abdominal imaging, PNS thresholds did not change significantly between YG1 and YG2 (-3.6%). The agreement between predicted and experimental PNS thresholds was within 11.4% normalized root mean square error for both coils and all landmarks. The PNS model also produced plausible predictions of the stimulation sites when compared to the sites of perception reported by the subjects. CONCLUSION: The PNS-optimization improved the PNS thresholds for the target scan landmark as well as most other studied landmarks, potentially yielding a significant improvement in image encoding performance that can be safely used in humans.

Peripheral nerve stimulation informed design of a high-performance asymmetric head gradient coil.

PURPOSE: Peripheral nerve stimulation (PNS) limits the image encoding performance of both body gradient coils and the latest generation of head gradients. We analyze a variety of head gradient design aspects using a detailed PNS model to guide the design process of a new high-performance asymmetric head gradient to raise PNS thresholds and maximize the usable image-encoding performance. METHODS: A novel three-layer coil design underwent PNS optimization involving PNS predictions of a series of candidate designs. The PNS-informed design process sought to maximize the usable parameter space of a coil with <10% nonlinearity in a 22 cm region of linearity, a relatively large inner diameter (44 cm), maximum gradient amplitude of 200 mT/m, and a high slew rate of 900 T/m/s. PNS modeling allowed identification and iterative adjustment of coil features with beneficial impact on PNS such as the number of winding layers, shoulder accommodation strategy, and level of asymmetry. PNS predictions for the final design were compared to measured thresholds in a constructed prototype. RESULTS: The final head gradient achieved up to 2-fold higher PNS thresholds than the initial design without PNS optimization and compared to existing head gradients with similar design characteristics. The inclusion of a third intermediate winding layer provided the additional degrees of freedom necessary to improve PNS thresholds without significant sacrifices to the other design metrics. CONCLUSION: Augmenting the design phase of a new high-performance head gradient coil by PNS modeling dramatically improved the usable image-encoding performance by raising PNS thresholds.

Prediction of experimental cardiac magnetostimulation thresholds using pig-specific body models.

PURPOSE: Modern high-amplitude gradient systems can be limited by the International Electrotechnical Commission 60601-2-33 cardiac stimulation (CS) limit, which was set in a conservative manner based on electrode experiments and E-field simulations in uniform ellipsoidal body models. Here, we show that coupled electromagnetic-electrophysiological modeling in detailed body and heart models can predict CS thresholds, suggesting that such modeling might lead to more detailed threshold estimates in humans. Specifically, we compare measured and predicted CS thresholds in eight pigs. METHODS: We created individualized porcine body models using MRI (Dixon for the whole body, CINE for the heart) that replicate the anatomy and posture of the animals used in our previous experimental CS study. We model the electric fields induced along cardiac Purkinje and ventricular muscle fibers and predict the electrophysiological response of these fibers, yielding CS threshold predictions in absolute units for each animal. Additionally, we assess the total modeling uncertainty through a variability analysis of the 25 main model parameters. RESULTS: Predicted and experimental CS thresholds agree within 19% on average (normalized RMS error), which is smaller than the 27% modeling uncertainty. No significant difference was found between the modeling predictions and experiments (p < 0.05, paired t-test). CONCLUSION: Predicted thresholds matched the experimental data within the modeling uncertainty, supporting the model validity. We believe that our modeling approach can be applied to study CS thresholds in humans for various gradient coils, body shapes/postures, and waveforms, which is difficult to do experimentally.

Mapping the human connectome using diffusion MRI at 300 mT/m gradient strength: Methodological advances and scientific impact.

Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength diffusion MRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for diffusion MRI and where the field is headed in the coming years.

A Huygens' surface approach to rapid characterization of peripheral nerve stimulation.

PURPOSE: Peripheral nerve stimulation (PNS) modeling has a potential role in designing and operating MRI gradient coils but requires computationally demanding simulations of electromagnetic fields and neural responses. We demonstrate compression of an electromagnetic and neurodynamic model into a single versatile PNS matrix (P-matrix) defined on an intermediary Huygens' surface to allow fast PNS characterization of arbitrary coil geometries and body positions. METHODS: The Huygens' surface approach divides PNS prediction into an extensive pre-computation phase of the electromagnetic and neurodynamic responses, which is independent of coil geometry and patient position, and a fast coil-specific linear projection step connecting this information to a specific coil geometry. We validate the Huygens' approach by performing PNS characterizations for 21 body and head gradients and comparing them with full electromagnetic-neurodynamic modeling. We demonstrate the value of Huygens' surface-based PNS modeling by characterizing PNS-optimized coil windings for a wide range of patient positions and poses in two body models. RESULTS: The PNS prediction using the Huygens' P-matrix takes less than a minute (instead of hours to days) without compromising numerical accuracy (error ≤ 0.1%) compared to the full simulation. Using this tool, we demonstrate that coils optimized for PNS at the brain landmark using a male model can also improve PNS for other imaging applications (cardiac, abdominal, pelvic, and knee imaging) in both male and female models. CONCLUSION: Representing PNS information on a Huygens' surface extended the approach's ability to assess PNS across body positions and models and test the robustness of PNS optimization in gradient design.

Measurement of magnetostimulation thresholds in the porcine heart.

PURPOSE: Powerful MRI gradient systems can surpass the International Electrotechnical Commission (IEC) 60601-2-33 limit for cardiac stimulation (CS), which was determined by simple electromagnetic simulations and electrode stimulation experiments. Only a few canine studies measured magnetically induced CS thresholds in vivo and extrapolating them to human safety limits can be challenging. METHODS: We measured cardiac magnetostimulation thresholds in 10 healthy, anesthetized pigs using capacitors discharged into a flat spiral coil to produce damped sinusoidal waveforms with effective stimulus duration ts,eff  = 0.45 ms. Electrocardiography (ECG), blood pressure, and peripheral oximetry signals were recorded to determine threshold coil currents yielding cardiac capture. Dixon and CINE MR volumes from each animal were segmented to generate porcine-specific electromagnetic models to calculate dB/dt and E-field values in the porcine heart at threshold. For comparison, we also simulated maximum dB/dt and E-field values created by three MRI gradient systems in the heart of a human body model. RESULTS: The average dB/dt threshold estimated in the porcine heart was 1.66 ± 0.23 kT/s, which is 11-fold greater than the IEC dB/dt limit at ts,eff  = 0.45 ms, and 31-fold greater than the maximum value created by the investigated MRI gradients in the human heart. The average E-field threshold estimated in the porcine heart was 92.9 ± 13.5 V/m, which is 6-fold greater than the IEC E-field limit at ts,eff  = 0.45 ms and 37-fold greater than the maximum gradient-induced E-field in the human heart. CONCLUSION: This first measurement of cardiac magnetostimulation thresholds in pigs indicates that the IEC cardiac safety limit is conservative for the investigated stimulus duration (ts,eff  = 0.45 ms).

A 31-channel integrated "AC/DC" B0 shim and radiofrequency receive array coil for improved 7T MRI.

PURPOSE: To test an integrated "AC/DC" array approach at 7T, where B0 inhomogeneity poses an obstacle for functional imaging, diffusion-weighted MRI, MR spectroscopy, and other applications. METHODS: A close-fitting 7T 31-channel (31-ch) brain array was constructed and tested using combined Rx and ΔB0 shim channels driven by a set of rapidly switchable current amplifiers. The coil was compared to a shape-matched 31-ch reference receive-only array for RF safety, signal-to-noise ratio (SNR), and inter-element noise correlation. We characterize the coil array's ability to provide global and dynamic (slice-optimized) shimming using ΔB0 field maps and echo planar imaging (EPI) acquisitions. RESULTS: The SNR and average noise correlation were similar to the 31-ch reference array. Global and slice-optimized shimming provide 11% and 40% improvements respectively compared to baseline second-order spherical harmonic shimming. Birdcage transmit coil efficiency was similar for the reference and AC/DC array setups. CONCLUSION: Adding ΔB0 shim capability to a 31-ch 7T receive array can significantly boost 7T brain B0 homogeneity without sacrificing the array's rdiofrequency performance, potentially improving ultra-high field neuroimaging applications that are vulnerable to off-resonance effects.

Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome.

The first phase of the Human Connectome Project pioneered advances in MRI technology for mapping the macroscopic structural connections of the living human brain through the engineering of a whole-body human MRI scanner equipped with maximum gradient strength of 300 mT/m, the highest ever achieved for human imaging. While this instrument has made important contributions to the understanding of macroscale connectional topology, it has also demonstrated the potential of dedicated high-gradient performance scanners to provide unparalleled in vivo assessment of neural tissue microstructure. Building on the initial groundwork laid by the original Connectome scanner, we have now embarked on an international, multi-site effort to build the next-generation human 3T Connectome scanner (Connectome 2.0) optimized for the study of neural tissue microstructure and connectional anatomy across multiple length scales. In order to maximize the resolution of this in vivo microscope for studies of the living human brain, we will push the diffusion resolution limit to unprecedented levels by (1) nearly doubling the current maximum gradient strength from 300 mT/m to 500 mT/m and tripling the maximum slew rate from 200 T/m/s to 600 T/m/s through the design of a one-of-a-kind head gradient coil optimized to minimize peripheral nerve stimulation; (2) developing high-sensitivity multi-channel radiofrequency receive coils for in vivo and ex vivo human brain imaging; (3) incorporating dynamic field monitoring to minimize image distortions and artifacts; (4) developing new pulse sequences to integrate the strongest diffusion encoding and highest spatial resolution ever achieved in the living human brain; and (5) calibrating the measurements obtained from this next-generation instrument through systematic validation of diffusion microstructural metrics in high-fidelity phantoms and ex vivo brain tissue at progressively finer scales with accompanying diffusion simulations in histology-based micro-geometries. We envision creating the ultimate diffusion MRI instrument capable of capturing the complex multi-scale organization of the living human brain - from the microscopic scale needed to probe cellular geometry, heterogeneity and plasticity, to the mesoscopic scale for quantifying the distinctions in cortical structure and connectivity that define cyto- and myeloarchitectonic boundaries, to improvements in estimates of macroscopic connectivity.

Investigating cardiac stimulation limits of MRI gradient coils using electromagnetic and electrophysiological simulations in human and canine body models.

PURPOSE: Cardiac stimulation (CS) limits to gradient coil switching speed are difficult to measure in humans; instead, current regulatory guidelines (IEC 60601-2-33) are based on animal experiments and electric field-to-dB/dt conversion factors computed for a simple, homogeneous body model. We propose improvement to this methodology by using more detailed CS modeling based on realistic body models and electrophysiological models of excitable cardiac fibers. METHODS: We compute electric fields induced by a solenoid, coplanar loops, and a commercial gradient coil in two human body models and a canine model. The canine simulations mimic previously published experiments. We generate realistic fiber topologies for the cardiac Purkinje and ventricular muscle fiber networks using rule-based algorithms, and evaluate CS thresholds using validated electrodynamic models of these fibers. RESULTS: We were able to reproduce the average measured canine CS thresholds within 5%. In all simulations, the Purkinje fibers were stimulated before the ventricular fibers, and therefore set the effective CS threshold. For the investigated gradient coil, simulated CS thresholds for the x-, y-, and z-axis were at least one order of magnitude greater than the International Electrotechnical Commission limit. CONCLUSION: We demonstrate an approach to simulate gradient-induced CS using a combination of electromagnetic and electrophysiological modeling. Pending additional validation, these simulations could guide the assessment of CS limits to MRI gradient coil switching speed. Such an approach may lead to less conservative, but still safe, operation limits, enabling the use of the maximum gradient amplitude versus slew rate parameter space of recent, powerful gradient systems.

Prediction of peripheral nerve stimulation thresholds of MRI gradient coils using coupled electromagnetic and neurodynamic simulations.

PURPOSE: As gradient performance increases, peripheral nerve stimulation (PNS) is becoming a significant constraint for fast MRI. Despite its impact, PNS is not directly included in the coil design process. Instead, the PNS characteristics of a gradient are assessed on healthy subjects after prototype construction. We attempt to develop a tool to inform coil design by predicting the PNS thresholds and activation locations in the human body using electromagnetic field simulations coupled to a neurodynamic model. We validate the approach by comparing simulated and experimentally determined thresholds for 3 gradient coils. METHODS: We first compute the electric field induced by the switching fields within a detailed electromagnetic body model, which includes a detailed atlas of peripheral nerves. We then calculate potential changes along the nerves and evaluate their response using a neurodynamic model. Both a male and female body model are used to study 2 body gradients and 1 head gradient. RESULTS: There was good agreement between the average simulated thresholds of the male and female models with the experimental average (normalized root-mean-square error: <10% and <5% in most cases). The simulation could also interrogate thresholds above those accessible by the experimental setup and allowed identification of the site of stimulation. CONCLUSIONS: Our simulation framework allows accurate prediction of gradient coil PNS thresholds and provides detailed information on location and "next nerve" thresholds that are not available experimentally. As such, we hope that PNS simulations can have a potential role in the design phase of high performance MRI gradient coils.

Fast three-dimensional inner volume excitations using parallel transmission and optimized k-space trajectories.

PURPOSE: To design short parallel transmission (pTx) pulses for excitation of arbitrary three-dimensional (3D) magnetization patterns. METHODS: We propose a joint optimization of the pTx radiofrequency (RF) and gradient waveforms for excitation of arbitrary 3D magnetization patterns. Our optimization of the gradient waveforms is based on the parameterization of k-space trajectories (3D shells, stack-of-spirals, and cross) using a small number of shape parameters that are well-suited for optimization. The resulting trajectories are smooth and sample k-space efficiently with few turns while using the gradient system at maximum performance. Within each iteration of the k-space trajectory optimization, we solve a small tip angle least-squares RF pulse design problem. Our RF pulse optimization framework was evaluated both in Bloch simulations and experiments on a 7T scanner with eight transmit channels. RESULTS: Using an optimized 3D cross (shells) trajectory, we were able to excite a cube shape (brain shape) with 3.4% (6.2%) normalized root-mean-square error in less than 5 ms using eight pTx channels and a clinical gradient system (Gmax = 40 mT/m, Smax = 150 T/m/s). This compared with 4.7% (41.2%) error for the unoptimized 3D cross (shells) trajectory. Incorporation of B0 robustness in the pulse design significantly altered the k-space trajectory solutions. CONCLUSION: Our joint gradient and RF optimization approach yields excellent excitation of 3D cube and brain shapes in less than 5 ms, which can be used for reduced field of view imaging and fat suppression in spectroscopy by excitation of the brain only. Magn Reson Med 76:1170-1182, 2016. © 2015 Wiley Periodicals, Inc.

Influence of peripheral axon geometry and local anatomy on magnetostimulation chronaxie.

Objective.Rapid switching of magnetic resonance imaging (MRI) gradient fields induces electric fields that can cause peripheral nerve stimulation (PNS) and so accurate characterization of PNS is required to maintain patient safety and comfort while maximizing MRI performance. The minimum magnetic gradient amplitude that causes stimulation, the PNS threshold, depends on intrinsic axon properties and the spatial and temporal properties of the induced electric field. The PNS strength-duration curve is widely used to characterize simulation thresholds for periodic waveforms and is parameterized by the chronaxie and rheobase. Safety limits to avoid unwanted PNS in MRI rely on a single chronaxie value to characterize the response of all nerves. However, experimental magnetostimulation peripheral nerve chronaxie values vary by an order of magnitude. Given the diverse range of chronaxies observed and the importance of this number in MRI safety models, we seek a deeper understanding of the mechanisms contributing to chronaxie variability.Approach.We use a coupled electromagnetic-neurodynamic PNS model to assess geometric sources of chronaxie variability. We study the impact of the position of the stimulating magnetic field coil relative to the body, along with the effect of local anatomical features and nerve trajectories on the driving function and the resulting chronaxie.Main results.We find realistic variation of local axon and tissue geometry can modulate a given axon's chronaxie by up to two-fold. Our results identify the temporal rate of charge redistribution as the underlying determinant of the chronaxie.Significance.This charge distribution is a function of both intrinsic axon properties and the spatial stimulus along the nerve; thus, examination of the local tissue topology, which shapes the electric fields, as well as the nerve trajectory, are critical for better understanding chronaxie variations and defining more biologically informed MRI safety guidelines.

Optimization of MRI Gradient Coils With Explicit Peripheral Nerve Stimulation Constraints.

Peripheral Nerve Stimulation (PNS) limits the acquisition rate of Magnetic Resonance Imaging data for fast sequences employing powerful gradient systems. The PNS characteristics are currently assessed after the coil design phase in experimental stimulation studies using constructed coil prototypes. This makes it difficult to find design modifications that can reduce PNS. Here, we demonstrate a direct approach for incorporation of PNS effects into the coil optimization process. Knowledge about the interactions between the applied magnetic fields and peripheral nerves allows the optimizer to identify coil solutions that minimize PNS while satisfying the traditional engineering constraints. We compare the simulated thresholds of PNS-optimized body and head gradients to conventional designs, and find an up to 2-fold reduction in PNS propensity with moderate penalties in coil inductance and field linearity, potentially doubling the image encoding performance that can be safely used in humans. The same framework may be useful in designing and operating magneto- and electro-stimulation devices.

Optimizing selective stimulation of peripheral nerves with arrays of coils or surface electrodes using a linear peripheral nerve stimulation metric.

OBJECTIVE: We present a PNS oracle, which solves these computation time and linearity problems and is, therefore, well-suited for fast optimization of voltage distributions in contact electrode arrays and current drive patterns in non-contact magnetic coil arrays. APPROACH: The PNS oracle metric for a nerve fiber is computed from an electric field map using only linear operations (projection, differentiation, convolution, scaling). Due to its linearity, this PNS metric can be precomputed for a set of coil or electrode segments, allowing rapid PNS prediction and comparison of any possible coil or electrode stimulation configuration constructed from this set. The PNS oracle is closely related to the classical activating function and modified driving functions but is adjusted to better correlate with full neurodynamic modeling of myelinated mammalian nerves. MAIN RESULTS: We validated the PNS oracle in three MRI gradient coils and two body models and found good correlation between the PNS oracle and the full neurodynamic modeling approach (R 2 > 0.995). Finally, we demonstrated its potential utility by optimizing the driving currents and voltages of arrays of 108 magnetic coils or 108 contact electrodes to selectively stimulate target nerves in the lower leg. SIGNIFICANCE: Peripheral nerve stimulation (PNS) by electromagnetic fields can be accurately simulated using coupled electromagnetic and neurodynamic modeling. Such simulations are slow and non-linear in the electric field, which makes it difficult to iteratively optimize coil and electrode configurations or drive patterns aiming to avoid PNS or to initiate it for therapeutic purposes.

The 'virtual DBS population': five realistic computational models of deep brain stimulation patients for electromagnetic MR safety studies.

We design, develop, and disseminate a 'virtual population' of five realistic computational models of deep brain stimulation (DBS) patients for electromagnetic (EM) analysis. We found five DBS patients in our institution' research patient database who received high quality post-DBS surgery computer tomography (CT) examinations of the head and neck. Three patients have a single implanted pulse generator (IPG) and the two others have two IPGs (one for each lead). Moreover, one patient has two abandoned leads on each side of the head. For each patient, we combined the head and neck volumes into a 'virtual CT', from which we extracted the full-length DBS path including the IPG, extension cables, and leads. We corrected topology errors in this path, such as self-intersections, using a previously published optimization procedure. We segmented the virtual CT volume into bones, internal air, and soft tissue classes and created two-manifold, watertight surface meshes of these distributions. In addition, we added a segmented model of the brain (grey matter, white matter, eyes and cerebrospinal fluid) to one of the model (nickname Freddie) that was derived from a T1-weighted MR image obtained prior to the DBS implantation. We simulated the EM fields and specific absorption rate (SAR) induced at 3 Tesla by a quadrature birdcage body coil in each of the five patient models using a co-simulation strategy. We found that inter-subject peak SAR variability across models was independent of the target averaging mass and equal to ~45%. In our simulations of the full brain segmentation and six simplified versions of the Freddie model, the error associated with incorrect dielectric property assignment around the DBS electrodes was greater than the error associated with modeling the whole model as a single tissue class. Our DBS patient models are freely available on our lab website (Webpage of the Martinos Center Phantom Resource 2018 https://phantoms.martinos.org/Main_Page).

Evaluating the effects of receive-only arrays in specific absorption rate simulations at 3 and 7 T.

Specific absorption rate (SAR) simulations are performed for most clinical and research transmit coil configurations. Such simulations allow the determination of limits in transmit power for patient safety. Different human models and coil configurations have been previously investigated using these simulations. However, only a few works have accounted for the effect of the receive (Rx) arrays in the SAR calculations and they have used very specialized setups or simplified detuning modeling of the Rx elements. In this work, we performed electromagnetic simulations using a clinical alike setup for whole-body scans at 3 T and head scans at 7 T. SAR simulations are performed for both setups with and without Rx arrays. A difference below 10% percent was found for max SAR. The maximum difference for the mean SAR values of the 3 T setups remained within 8% and within 15% of the 7 T setup.

Sensitivity analysis of neurodynamic and electromagnetic simulation parameters for robust prediction of peripheral nerve stimulation.

Peripheral nerve stimulation (PNS) has become an important limitation for fast MR imaging using the latest gradient hardware. We have recently developed a simulation framework to predict PNS thresholds and stimulation locations in the body for arbitrary coil geometries to inform the gradient coil optimization process. Our approach couples electromagnetic field simulations in realistic body models to a neurodynamic model of peripheral nerve fibers. In this work, we systematically analyze the impact of key parameters on the predicted PNS thresholds to assess the robustness of the simulation results. We analyze the sensitivity of the simulated thresholds to variations of the most important simulation parameters, including parameters of the electromagnetic field simulations (dielectric tissue properties, body model size, position, spatial resolution, and coil model discretization) and parameters of the neurodynamic simulation (length of the simulated nerves, position of the nerve model relative to the extracellular potential, temporal resolution of the nerve membrane dynamics). We found that for the investigated setup, the subject-dependent parameters (e.g. tissue properties or body size) can affect PNS prediction by up to ~26% when varied in a natural range. This is in accordance with the standard deviation of ~30% reported in human subject studies. Parameters related to numerical aspects can cause significant simulation errors (>30%), if not chosen cautiously. However, these perturbations can be controlled to yield errors below 5% for all investigated parameters without an excessive increase in computation time. Our sensitivity analysis shows that patient-specific parameter fluctuations yield PNS threshold variations similar to the variations observed in experimental PNS studies. This may become useful to estimate population-average PNS thresholds and understand their standard deviation. Our analysis indicates that the simulated PNS thresholds are numerically robust, which is important for ranking different MRI gradient coil designs or assessing different PNS mitigation strategies.

Evaluation of stacked resonators to enhance the performance of a surface receive-only array for prostate MRI at 3 Tesla.

Prostate MRI is an important tool to diagnose and characterize cancer. High local sensitivity and good parallel imaging performance are of paramount importance for diagnostic quality and efficiency. The purpose of this work was to evaluate stacked resonators as part of a surface receiver array for prostate MRI at 3 Tesla. A base array of 6-channels consisting of a flexible anterior and a rigid posterior part were built each with three loop coils. A pair of stacked resonators was added concentrically to the center loops (anterior and posterior) of the base array. The evaluated stacked resonators were butterflies, composites and dipoles which yielded a total of three 8-channel arrays. The arrays were compared using noise correlations and single-channel signal-to-noise ratio maps in a phantom. Combined signal-to-noise ratio maps and parallel imaging performances were measured and compared in vivo in 6 healthy volunteers. The results were compared to the base and a commercial array. The SNR values in the prostate yielded by all the arrays were not statistically different using fully sampled k-space. However, significant differences were found in the parallel imaging performance of the arrays. More specifically, up to 88% geometric factor reduction was found compared to the commercial array and up to 83% reduction compared to the base array using butterfly coils. Thus, signal-to-noise ratio improvements were observed with stacked resonators when using parallel imaging. The use of stacked elements, in particular butterfly coils, can improve the performance of a base array consisting solely of single loops when using parallel imaging. We expect prostate MRI at 3 Tesla to benefit from using combinations of single loops and stacked resonators.