Modeling and measurement of lead tip heating and resonant length for implanted, insulated wires.

PURPOSE: To study implant lead tip heating because of the RF power deposition by developing mathematical models and comparing them with measurements acquired at 1.5 T and 3 T, especially to predict resonant length. THEORY AND METHODS: A simple exponential model and an adapted transmission line model for the electric field transfer function were developed. A set of wavenumbers, including that calculated from insulated antenna theory (King wavenumber) and that of the embedding medium were considered. Experiments on insulated, capped wires of varying lengths were performed to determine maximum temperature rise under RF exposure. The results are compared with model predictions from analytical expressions derived under the assumption of a constant electric field, and with those numerically calculated from spatially varying, simulated electric fields from body coil transmission. Simple expressions for the resonant length bounded between one-quarter and one-half wavelength are developed based on the roots of transcendental equations. RESULTS: The King wavenumber for both models more closely matched the experimental data with a maximum root mean square error of 9.81°C at 1.5 T and 5.71°C at 3 T compared to other wavenumbers with a maximum root mean square error of 27.52°C at 1.5 T and 22.01°C for 3 T. Resonant length was more accurately predicted compared to values solely based on the embedding medium. CONCLUSION: Analytical expressions were developed for implanted lead heating and resonant lengths under specific assumptions. The value of the wavenumber has a strong effect on the model predictions. Our work could be used to better manage implanted device lead tip heating.

Retrospective correction of second-order concomitant fields in 3D axial stack-of-spirals imaging on a high-performance gradient system.

PURPOSE: MRI using 3D stack-of-spirals (SoS) readout on a high-performance gradient system is subject to strong second-order, spatially varying concomitant fields, which can lead to signal dropout and blurring artifacts that become more significant at locations farther from the gradient isocenter. A method for compensating for second-order concomitant fields in 3D axial SoS image reconstruction is described. METHODS: We retrospectively correct for second-order concomitant field-induced phase error in the 3D SoS data by slice-dependent k-space phase compensation based on the nominal spiral readout trajectories. The effectiveness of the method was demonstrated in phantom and healthy volunteer scans in which 3D pseudo-continuous arterial spin labeling imaging was performed with SoS fast spin-echo readout at 3 T. RESULTS: Substantial reduction in blurring was observed with the proposed method. In phantom scans, blurring was reduced by about 53% at 98 mm from the gradient isocenter. In the in vivo 3D pseudo-continuous arterial spin labeling scans, differences of up to 10% were observed at 78 mm from the isocenter, especially around the white-matter and gray-matter interfaces, between the corrected and uncorrected proton density images, perfusion-weighted images, and cerebral blood flow maps. CONCLUSIONS: The described retrospective correction method provides a means to correct erroneous phase accruals due to second-order concomitant fields in 3D axial stack-of-spirals imaging.

Generalized magnetostatic target field method for shielded magnetic field coils in a separable coordinate system.

A theoretical method is described to analytically calculate a pair of surface current densities, which produce a desired static magnetic field in one region of the space and zero magnetic field in another. The analysis is based on the known relationship between a surface current density and a stream function, the equivalence of stream functions and surface magnetic dipole density, and the scalar potential representation of the associated magnetic field in free space. From these relations, we formulate the magnetostatic problem, which is often treated as a vector field problem, as a scalar field problem in which a two-dimensional scalar field (stream function) is related to a three-dimensional one (magnetic scalar potential) via the differentiation of the electrostatic Green's function 1/|r-rs|. It is shown that, in a coordinate system in which a separated form of the Green's function exists (separable coordinate system), there exists a simple relationship between a harmonic component of a stream function and a harmonic component of the magnetic scalar potential. The method is applied to calculate idealized surface current patterns for actively shielded, linear gradient field coils in the Cartesian, cylindrical, and spherical coordinates.

Calibration of concomitant field offsets using phase contrast MRI for asymmetric gradient coils.

PURPOSE: Asymmetric gradient coils introduce zeroth- and first-order concomitant field terms, in addition to higher-order terms common to both asymmetric and symmetric gradients. Salient to compensation strategies is the accurate calibration of the concomitant field spatial offset parameters for asymmetric coils. A method that allows for one-time calibration of the offset parameters is described. THEORY AND METHODS: A modified phase contrast pulse sequence with single-sided bipolar flow encoding is proposed to calibrate the offsets for asymmetric, transverse gradient coils. By fitting the measured phase offsets to different gradient amplitudes, the spatial offsets were calculated by fitting the phase variation. This was used for calibrating real-time pre-emphasis compensation of the zeroth- and first-order concomitant fields. RESULTS: Image quality improvement with the proposed corrections was demonstrated in phantom and healthy volunteers with non-Cartesian and Cartesian trajectory acquisitions. Concomitant field compensation using the calibrated offsets resulted in a residual phase error <3% at the highest gradient amplitude and demonstrated substantial reduction of image blur and slice position/selection artifacts. CONCLUSIONS: The proposed implementation provides an accurate method for calibrating spatial offsets that can be used for real-time concomitant field compensation of zeroth and first-order terms, substantially reducing artifacts without retrospective correction or sequence specific waveform modifications.

In vivo direct imaging of neuronal activity at high temporospatial resolution.

There has been a long-standing demand for noninvasive neuroimaging methods that can detect neuronal activity at both high temporal and high spatial resolution. We present a two-dimensional fast line-scan approach that enables direct imaging of neuronal activity with millisecond precision while retaining the high spatial resolution of magnetic resonance imaging (MRI). This approach was demonstrated through in vivo mouse brain imaging at 9.4 tesla during electrical whisker-pad stimulation. In vivo spike recording and optogenetics confirmed the high correlation of the observed MRI signal with neural activity. It also captured the sequential and laminar-specific propagation of neuronal activity along the thalamocortical pathway. This high-resolution, direct imaging of neuronal activity will open up new avenues in brain science by providing a deeper understanding of the brain's functional organization, including the temporospatial dynamics of neural networks.

Systematic Dimensional Analysis of the Scaling Relationship for Gradient and Shim Coil Design Parameters.

PURPOSE: To demonstrate systematic, linear algebra-based, dimensional analysis to derive a scaling relationship among the design parameters of MRI gradient and harmonic shim coils. THEORY AND METHODS: The dimensions of five physical quantities relevant for gradient coil design (inductance, gradient amplitude, inner diameter [ d$$ d $$ ], current, and the permeability of free space) were decomposed into fundamental units, and their exponents were arranged into a dimensional matrix. The resulting set of homogenous equations was solved using standard linear algebraic methods. Inclusion of the number of turns as an additional unit yielded a 5 × 5 dimensional matrix with a unique, nontrivial solution. The analysis was extended to harmonic shim coils. The gradient coil scaling relationship was compared with data from 24 published gradient coil sets. RESULTS: Only when the unit of turns was included did the linear algebra-based analysis uniquely produce the known scaling relationship that gradient inductance is proportional to gradient efficiency squared times d5$$ {d}^5 $$ . By applying the same methodology to an lth order shim coil, a novel result is obtained: Shim inductance is proportional to its efficiency squared times d2l+3$$ {d}^{2l+3} $$ . The predicted power-law relationship between inductance-normalized gradient efficiency and the diameter accounted for > 92% of the efficiency variation of the surveyed gradient coils. A dimensionless parameter is proposed as an intrinsic figure-of-merit of gradient coil efficiency. CONCLUSION: Systematic application of linear algebra-based dimensional analysis can provide new insight in gradient and shim coil design by revealing fundamental scaling relations and helping to guide the design and comparison of coils with different diameters.

Therapeutic Quadrisected Annular Array for Improving Magnetic Resonance Compatibility.

OBJECTIVE: Focused ultrasound has been applied in brain therapeutics. Although focusing ultrasonic beams on multiple arbitrary regions under the guidance of magnetic resonance imaging(MRI) is needed for precise treatments, current therapeutic transducers with large pitch sizes have been optimized to focus on deep brain regions. While annular arrays can adjust the beam foci from cortical to deep regions, their circular shape may generate eddy current-induced magnetic flux during MRI. In this study, a quadrisected annular array is proposed to address these limitations. METHODS: Conventional and quadrisected annular arrays with three elements were implemented by loading the electrode patterns onto an 850 kHz 1-3 composite PZT disc, with a diameter of 31 mm, including three rings. MR compatibilities were demonstrated by imaging an MRI phantom with pulse sequences for B0 and B1 mapping and spin-echo imaging. Acoustic beam profiles, with and without a macaque monkey skull, were measured. A quadrisected transducer was also used to open the blood-brain barrier(BBB). RESULTS: The flip angle distortion improved by 20% in spin-echo MR imaging. The acoustic beam distortions shifting the focal point from 36 to 41mm and elongating the focal zone from 10 to 15 mm could be recovered to nearly the original values. BBB openings in the hippocampus and basal region were also demonstrated. CONCLUSION: The MR compatibility was improved by the increased resistance of the electrodes in the quadrisected array maintaining dynamic focusing capabilities. SIGNIFICANCE: The quadrisected annular design can be a fundamental structure for a larger MR-compatible segmented array transducer generating multiple acoustic foci.

Rapid calculation of static magnetic field perturbation generated by magnetized objects in arbitrary orientations.

PURPOSE: Most previous work on the calculation of susceptibility-induced static magnetic field (B0 ) inhomogeneity has considered strictly unidirectional magnetic fields. Here, we present the theory and implementation of a computational method to rapidly calculate static magnetic field vectors produced by an arbitrary distribution of voxelated magnetization vectors. THEORY AND METHODS: Two existing B0 calculation methods were systematically extended to include arbitrary orientations of the magnetization and the magnetic field; they are (1) Fourier-domain convolution with k-space-discretized (KD) dipolar field, and (2) generalized susceptibility voxel convolution (gSVC). The methods were tested on an analytical ellipsoid model and a tilted human head model, as well as against experimentally measured B0 fields induced by a stainless-steel implant located in an inhomogeneous region of a clinical 3T MRI magnet. RESULTS: Both methods were capable of correctly calculating B0 fields inside a magnetized ellipsoid in all tested orientations. The KD method generally required a larger grid and longer computation time to achieve accuracy comparable to gSVC. Measured B0 fields due to the implant showed a good match with the gSVC-calculated fields that accounted for the spatial variation of the applied magnetic field including the radial components. CONCLUSION: Our method can provide a reliable and efficient computational tool to calculate B0 perturbation by magnetized objects under a variety of circumstances, including those with inhomogeneous magnetizing fields, anisotropic susceptibility, and a rotated coordinate system.

Extra-basal ganglia iron content and non-motor symptoms in drug-naïve, early Parkinson's disease.

BACKGROUND: Although iron dyshomeostasis is associated with Parkinson's disease (PD) pathogenesis, the relationship between iron deposition and non-motor involvement in PD is not fully understood. In this study, we investigated basal ganglia and extra-basal ganglia system iron contents and their correlation with non-motor symptoms in drug-naïve, early-stage PD patients. METHODS: We enrolled 14 drug-naïve, early-stage PD patients and 12 age/sex-matched normal controls. All participants underwent brain magnetic resonance imaging to obtain the effective transverse relaxation rate (R2*) and quantitative susceptibility mapping (QSM). Deep brain structures, including the nucleus accumbens, caudate nucleus, putamen, globus pallidus, thalamus, hippocampus, and amygdala, were delineated using the FSL-FIRST; the substantia nigra, red nucleus, and dentate nucleus were segmented manually. Inter-group differences in R2* and QSM values, as well as their association with clinical parameters of PD, were investigated. RESULTS: Substantia nigra and putamen R2* values were significantly higher in PD patients than in normal controls, despite no significant difference in QSM values. Regarding the non-motor symptom scales, PD sleep scale score negatively correlated with R2* values in the red nucleus and right amygdala, Scales for Outcomes in Parkinson's disease-Autonomic scores were positively correlated with R2* values in the right amygdala and left hippocampus, and cardiovascular sub-score of Non-Motor Symptoms Scale for PD was positively associated with the QSM value in the left hippocampus. CONCLUSION: In this study, iron content in the extra-basal ganglia system was significantly correlated with non-motor symptoms, especially sleep problems and dysautonomia, even in early-stage PD.

Line scan-based rapid magnetic resonance imaging of repetitive motion.

Two-dimensional (2D) line scan-based dynamic magnetic resonance imaging (MRI) is examined as a means to capture the interior of objects under repetitive motion with high spatiotemporal resolutions. The method was demonstrated in a 9.4-T animal MRI scanner where line-by-line segmented k-space acquisition enabled recording movements of an agarose phantom and quail eggs in different conditions-raw and cooked. A custom MR-compatible actuator which utilized the Lorentz force on its wire loops in the scanner's main magnetic field effectively induced the required periodic movements of the objects inside the magnet. The line-by-line k-space segmentation was achieved by acquiring a single k-space line for every frame in a motion period before acquisition of another line with a different phase-encode gradient in the succeeding motion period. The reconstructed time-course images accurately represented the objects' displacements with temporal resolutions up to 5.5 ms. The proposed method can drastically increase the temporal resolution of MRI for imaging rapid periodic motion of objects while preserving adequate spatial resolution for internal details when their movements are driven by a reliable motion-inducing mechanism.

An equal-TE ultrafast 3D gradient-echo imaging method with high tolerance to magnetic susceptibility artifacts: Application to BOLD functional MRI.

PURPOSE: To develop an ultrafast 3D gradient echo-based MRI method with constant TE and high tolerance to B0 inhomogeneity, dubbed ERASE (equal-TE rapid acquisition with sequential excitation), and to introduce its use in BOLD functional MRI (fMRI). THEORY AND METHODS: Essential features of ERASE, including spin behavior, were characterized, and a comparison study was conducted with conventional EPI. To demonstrate high tolerance to B0 inhomogeneity, in vivo imaging of the mouse brain with a fiber-optic implant was performed at 9.4 T, and human brain imaging (including the orbitofrontal cortex) was performed at 3 T and 7 T. To evaluate the performance of ERASE in BOLD-fMRI, the characteristics of SNR and temporal SNR were analyzed for in vivo rat brains at 9.4 T in comparison with multislice gradient-echo EPI. Percent signal changes and t-scores are also presented. RESULTS: For both mouse brain and human brain imaging, ERASE exhibited a high tolerance to magnetic susceptibility artifacts, showing much lower distortion and signal dropout, especially in the regions involving large magnetic susceptibility effects. For BOLD-fMRI, ERASE provided higher temporal SNR and t-scores than EPI, but exhibited similar percent signal changes in in vivo rat brains at 9.4 T. CONCLUSION: When compared with conventional EPI, ERASE is much less sensitive, not only to EPI-related artifacts such as Nyquist ghosting, but also to B0 inhomogeneity including magnetic susceptibility effects. It is promising for use in BOLD-fMRI, providing higher temporal SNR and t-scores with constant TE when compared with EPI, although further optimization is needed for human fMRI.

Radio-Frequency Vector Magnetic Field Mapping in Magnetic Resonance Imaging.

A method is presented to measure the radio-frequency (RF) vector magnetic field inside an object using magnetic resonance imaging (MRI). Conventional " [Formula: see text] mapping" in MRI can measure the proton co-rotating component ( [Formula: see text] of the RF field produced by a transmit coil. Here we show that by repeating [Formula: see text] mapping on the same object and coil at multiple (8) specific orientations with respect to the main magnet, the magnitudes and relative phases of all (x, y, z) Cartesian components of the RF field can be determined unambiguously. We demonstrate the method on a circularly polarized volume coil and a loop coil tuned at 123.25 MHz in a 3 Tesla MRI scanner, with liquid phantoms. The volume coil measurement showed the axial component of the RF field, which is normally unmeasurable in MRI, away from the center of the coil. The measured RF vector field maps of both coils compared favorably with numerical simulation, with volumetric normalized root-mean-square difference in the range of 7~20%. While the proposed method cannot be applied to human imaging at present, applications to phantoms and small animals could provide a useful experimental tool to validate RF simulation and verify certain assumptions in [Formula: see text] map-based electrical properties tomography (EPT).

Magnetic susceptibility imaging of human habenula at 3 T.

The habenula plays an important role in brain reward circuitry and psychiatric conditions. While much work has been done on the function and structure of the habenula in animal models, in vivo imaging studies of the human habenula have been relatively scarce due to its small size, deep brain location, and lack of clear biomarkers for its heterogeneous substructure. In this paper, we report high-resolution (0.5 × 0.5 × 0.8 mm3) MRI of the human habenula with quantitative susceptibility mapping (QSM) at 3 T. By analyzing 48 scan datasets collected from 21 healthy subjects, we found that magnetic susceptibility contrast is highly non-uniform within the habenula and across the subjects. In particular, we observed high prevalence of elevated susceptibility in the posterior subregion of the habenula. Correlation analysis between the susceptibility and the effective transverse relaxation rate (R2*) indicated that localized susceptibility enhancement in the habenula is more associated with increased paramagnetic (such as iron) rather than decreased diamagnetic (such as myelin) sources. Our results suggest that high-resolution QSM could make a potentially useful tool for substructure-resolved in vivo habenula imaging, and provide a groundwork for the future development of magnetic susceptibility as a quantitative biomarker for human habenula studies.

Strategies for rapid reconstruction in 3D MRI with radial data acquisition: 3D fast Fourier transform vs two-step 2D filtered back-projection.

For 3D radial data reconstruction in magnetic resonance imaging (MRI), fast Fourier transform via gridding (gFFT) is widely used for its fast processing and flexibility. In comparison, conventional 3D filtered back projection (cFBP), while more robust against common radial k-space centering errors, suffers from long computation times and is less frequently used. In this study, we revisit another back-projection reconstruction strategy, namely two-step 2D filtered back-projection (tsFBP), as an alternative 3D radial MRI reconstruction method that combines computational efficiency and certain error tolerance. In order to compare the three methods (gFFT, cFBP, and tsFBP), theoretical analysis was performed to evaluate the number of computational steps involved in each method. Actual reconstruction times were also measured and compared using 3D radial-MRI data of a phantom and a human brain. Additionally, the sensitivity of tsFBP to artifacts caused by radial k-space centering errors was compared with the other methods. Compared to cFBP, tsFBP dramatically improved the reconstruction speed while retaining the benefit of tolerance to the radial k-space errors. Our study therefore suggests that tsFBP can be a promising alternative to the conventional back projection method for 3D radial MRI reconstruction.

Feasibility of head-tilted brain scan to reduce susceptibility-induced signal loss in the prefrontal cortex in gradient echo-based imaging.

Susceptibility-induced static field (B0) inhomogeneity near the nasal cavity degrades high-field MRI image quality. Many studies have addressed this problem by hardware- or sequence-based methods to improve local B0 shimming or minimize the impact of inhomogeneity. Here, we investigate the feasibility of the head-tilted brain scan as an easily accessible way to reduce B0 inhomogeneity and associated gradient echo signal loss in the prefrontal cortex (PFC). We exploit the fact that the region of intense local B0 gradient can be steered away from the PFC by head reorientation with respect to the main magnetic field. We found that the required chin-up head tilting by a substantial angle (> 30°) can be readily achieved for a group of healthy subjects when their back was raised by about 10  cm. Eleven subjects were scanned at 3T, using a standard 20 channel head-neck coil, for whole-head B0 mapping and gradient-echo EPI-based functional MRI (fMRI) performing a reward-punishment task in normal and tilted head orientations. Additionally, multi-echo gradient echo and resting-state fMRI scans were performed on six subjects in both orientations. Head-tilted sessions, which lasted for at least 20 min, were well-tolerated by all subjects and demonstrated a marked reduction of localized signal loss in the gradient echo-based images and EPI images in the PFC compared to normal orientation scans. Imaging in tilted orientation reduced the group-averaged B0 standard deviation and peak B0 gradient in the orbital gyrus beyond what was possible with simulated 3rd order shimming. The behavioral performance in the head-tilted fMRI scans indicated that the subjects were able to perform a cognitive task with little difficulty, and the tilted fMRI scans successfully produced a robust whole-brain functional activation map consistent with the literature. Our study proposes that the back-raised, head-tilted imaging can benefit the shimming of the prefrontal brain regions while being compatible with moderate-length neuroimaging scans on healthy, cooperating subjects.

Experimental setup for bulk susceptibility effect-minimized, multi-orientation MRI of ex vivo tissue samples.

PURPOSE: Many conventional ex vivo magnetic resonance imaging (MRI) setups utilize cylindrical or other nonspherical tissue containers which can cause static field (B0 ) inhomogeneity affecting the accuracy of the measurements in an orientation-dependent manner. In this work we demonstrate an experimental method to obtain MRI of ex vivo tissue samples held in a spherical container in order to minimize bulk susceptibility-induced B0 inhomogeneity in arbitrary orientations. METHODS: B0 inhomogeneity caused by tissue-air susceptibility mismatch can be theoretically eliminated if the surface of susceptibility discontinuity is spherical. This situation can be approximated by putting a tissue sample in a spherical shell filled with materials with tissue-like magnetic susceptibility. We achieved this on an intact monkey brain by (a) holding the brain with a three-dimensional (3D)-printed holder with tissue-like (within 0.5 ppm) susceptibility, and (b) enclosing the brain and the holder in an acrylic spherical shell filled with diamagnetic liquid. Furthermore, the sphere and the radio-frequency coil for MRI were mounted on a 3D-printed frame designed to reduce B0 inhomogeneity contributions. The sphere could be rotated freely without disturbing the RF coil to facilitate multi-orientation imaging. We verified our setup by mapping B0 in the monkey brain at 13 different orientations in a human 7T scanner, and measuring orientation-dependent R 2 ∗ relaxation rates in the white and gray matters of the brain. The results were then compared with a setup where the brain was held inside a cylindrical container. RESULTS: In all orientations, the B0 standard deviation in the brain in the spherical setup (converted to Larmor frequency offset) was less than about 10 Hz. This corresponds to two-sigma deviation of B0 of <0.07 ppm. The B0 gradient was <9 Hz/mm in 95 % of the brain voxels in all orientations. In high-resolution imaging with e.g. voxel size <0.4 mm, this corresponds to voxel line broadening of <4 Hz (0.013 ppm). R 2 ∗ in the corpus callosum showed distinctly different orientation dependence compared to the gray matter. The B0 uniformity and R 2 ∗ reliability were much reduced in the cylindrical container setup. CONCLUSIONS: We have demonstrated an experimental method to effectively minimize bulk susceptibility-induced B0 perturbation in multi-orientation ex vivo MRI. The method promises to benefit a range of tissue orientation-dependent MR property studies.

Fabrication of a spherical inclusion phantom for validation of magnetic resonance-based magnetic susceptibility imaging.

Fabrication of a spherical multi-compartment MRI phantom is demonstrated that can be used to validate magnetic resonance (MR)-based susceptibility imaging reconstruction. The phantom consists of a 10 cm diameter gelatin sphere that encloses multiple smaller gelatin spheres doped with different concentrations of paramagnetic contrast agents. Compared to previous multi-compartment phantoms with cylindrical geometry, the phantom provides the following benefits: (1) no compartmental barrier materials are used that can introduce signal voids and spurious phase; (2) compartmental geometry is reproducible; (3) spherical susceptibility boundaries possess a ground-truth analytical phase solution for easy experimental validation; (4) spherical geometry of the overall phantom eliminates background phase due to air-phantom boundary in any scan orientation. The susceptibility of individual compartments can be controlled independently by doping. During fabrication, formalin cross-linking and water-proof surface coating effectively blocked water diffusion between the compartments to preserve the phantom's integrity. The spherical shapes were realized by molding the inner gel compartments in acrylic spherical shells, 3 cm in diameter, and constructing the whole phantom inside a larger acrylic shell. From gradient echo images obtained at 3T, we verified that the phantom produced phase images in agreement with the theoretical prediction. Factors that limit the agreement include: air bubbles trapped at the gel interfaces, imperfect magnet shimming, and the susceptibility of external materials such as the phantom support hardware. The phantom images were used to validate publicly available codes for quantitative susceptibility mapping. We believe that the proposed phantom can provide a useful testbed for validation of MR phase imaging and MR-based magnetic susceptibility reconstruction.

Electrical properties tomography: Available contrast and reconstruction capabilities.

MR-based electrical properties tomography converts the MRI transmit/receive RF field measurements to tissue electrical property maps through dedicated reconstruction algorithms. Recent reports showed that despite limitations, electrical properties tomography holds promise for generating additional contrast for tumor detection and patient-specific modeling of tissue-RF field interactions. This review summarizes the available tissue electrical property contrasts and compares them with the capabilities of the most commonly used electrical properties tomography reconstruction method. Future directions and prospects of clinical translation are discussed.

Lightweight, compact, and high-performance 3T MR system for imaging the brain and extremities.

PURPOSE: To build and evaluate a small-footprint, lightweight, high-performance 3T MRI scanner for advanced brain imaging with image quality that is equal to or better than conventional whole-body clinical 3T MRI scanners, while achieving substantial reductions in installation costs. METHODS: A conduction-cooled magnet was developed that uses less than 12 liters of liquid helium in a gas-charged sealed system, and standard NbTi wire, and weighs approximately 2000 kg. A 42-cm inner-diameter gradient coil with asymmetric transverse axes was developed to provide patient access for head and extremity exams, while minimizing magnet-gradient interactions that adversely affect image quality. The gradient coil was designed to achieve simultaneous operation of 80-mT/m peak gradient amplitude at a slew rate of 700 T/m/s on each gradient axis using readily available 1-MVA gradient drivers. RESULTS: In a comparison of anatomical imaging in 16 patients using T2 -weighted 3D fluid-attenuated inversion recovery (FLAIR) between the compact 3T and whole-body 3T, image quality was assessed as equivalent to or better across several metrics. The ability to fully use a high slew rate of 700 T/m/s simultaneously with 80-mT/m maximum gradient amplitude resulted in improvements in image quality across EPI, DWI, and anatomical imaging of the brain. CONCLUSIONS: The compact 3T MRI system has been in continuous operation at the Mayo Clinic since March 2016. To date, over 200 patient studies have been completed, including 96 comparison studies with a clinical 3T whole-body MRI. The increased gradient performance has reliably resulted in consistently improved image quality.

Rapid, theoretically artifact-free calculation of static magnetic field induced by voxelated susceptibility distribution in an arbitrary volume of interest.

PURPOSE: To demonstrate a computationally efficient and theoretically artifact-free method to calculate static field (B0 ) inhomogeneity in a volume of interest induced by an arbitrary voxelated susceptibility distribution. METHODS: Our method computes B0 by circular convolution between a zero-filled susceptibility matrix and a shifted, voxel-integrated dipolar field kernel on a grid of size NS +NT - 1 in each dimension, where NS and NT are the sizes of the susceptibility source and B0 target grids, respectively. The computational resource requirement is independent of source-target separation. The method, called generalized susceptibility voxel convolution, is demonstrated on three susceptibility models: an ellipsoid, MR-compatible screws, and a dynamic human heartbeat model. RESULTS: B0 in an ellipsoid calculated by generalized susceptibility voxel convolution matched an analytical solution nearly exactly. The method also calculated screw-induced B0 in agreement with experimental data. Dynamic simulation demonstrated its computational efficiency for repeated B0 calculations on time-varying susceptibility. On the contrary, conventional and alias-subtracted k-space-discretized Fourier convolution methods showed nonnegligible aliasing and Gibbs ringing artifacts in the tested models. CONCLUSION: Generalized susceptibility voxel convolution can be a fast and reliable way to compute susceptibility-induced B0 when the susceptibility source is not colocated with the B0 target volume of interest, as in modeling B0 variations from motion and foreign objects.