Intra-scan RF power amplifier drift correction.

PURPOSE: The drift in radiofrequency (RF) power amplifiers (RFPAs) is assessed and several contributing factors are investigated. Two approaches for prospective correction of drift are proposed and their effectiveness is evaluated. METHODS: RFPA drift assessment encompasses both intra-pulse and inter-pulse drift analyses. Scan protocols with varying flip angle (FA), RF length, and pulse repetition time (TR) are used to gauge the influence of these parameters on drift. Directional couplers (DICOs) monitor the forward waveforms of the RFPA outputs. DICOs data is stored for evaluation, allowing calculation of correction factors to adjust RFPAs' transmit voltage. Two correction methods, predictive and run-time, are employed: predictive correction necessitates a calibration scan, while run-time correction calculates factors during the ongoing scan. RESULTS: RFPA drift is indeed influenced by the RF duty-cycle, and in the cases examined with a maximum duty-cycle of 66%, the potential drift is approximately 41% or 15%, depending on the specific RFPA revision. Notably, in low transmit voltage scenarios, FA has minimal impact on RFPA drift. The application of predictive and run-time drift correction techniques effectively reduces the average drift from 10.0% to less than 1%, resulting in enhanced MR signal stability. CONCLUSION: Utilizing DICO recordings and implementing a feedback mechanism enable the prospective correction of RFPA drift. Having a calibration scan, predictive correction can be utilized with fewer complexity; for enhanced performance, a run-time approach can be employed.

Mitigating susceptibility-induced distortions in high-resolution 3DEPI fMRI at 7T.

Geometric distortion is a major limiting factor for spatial specificity in high-resolution fMRI using EPI readouts and is exacerbated at higher field strengths due to increased B0 field inhomogeneity. Prominent correction schemes are based on B0 field-mapping or acquiring reverse phase-encoded (reversed-PE) data. However, to date, comparisons of these techniques in the context of fMRI have only been performed on 2DEPI data, either at lower field or lower resolution. In this study, we investigate distortion compensation in the context of sub-millimetre 3DEPI data at 7T. B0 field-mapping and reversed-PE distortion correction techniques were applied to both partial coverage BOLD-weighted and whole brain MT-weighted 3DEPI data with matched distortion. Qualitative assessment showed overall improvement in cortical alignment for both correction techniques in both 3DEPI fMRI and whole-brain MT-3DEPI datasets. The distortion-corrected MT-3DEPI images were quantitatively evaluated by comparing cortical alignment with an anatomical reference using dice coefficient (DC) and correlation ratio (CR) measures. These showed that B0 field-mapping and reversed-PE methods both improved correspondence between the MT-3DEPI and anatomical data, with more substantial improvements consistently obtained using the reversed-PE approach. Regional analyses demonstrated that the largest benefit of distortion correction, and in particular of the reversed-PE approach, occurred in frontal and temporal regions where susceptibility-induced distortions are known to be greatest, but had not led to complete signal dropout. In conclusion, distortion correction based on reversed-PE data has shown the greater capacity for achieving faithful alignment with anatomical data in the context of high-resolution fMRI at 7T using 3DEPI.

Prediction of motion induced magnetic fields for human brain MRI at 3 T.

OBJECTIVE: Maps of B0 field inhomogeneities are often used to improve MRI image quality, even in a retrospective fashion. These field inhomogeneities depend on the exact head position within the static field but acquiring field maps (FM) at every position is time consuming. Here we propose a forward simulation strategy to obtain B0 predictions at different head-positions. METHODS: FM were predicted by combining (1) a multi-class tissue model for estimation of tissue-induced fields, (2) a linear k-space model for capturing gradient imperfections, (3) a dipole estimation for quantifying lower-body perturbing fields (4) and a position-dependent tissue mask to model FM alterations caused by large motion effects. The performance of the combined simulation strategy was compared with an approach based on a rigid body transformation of the FM measured in the reference position to the new position. RESULTS: The transformed FM provided inconsistent results for large head movements (> 5° rotation, approximately), while the simulation strategy had a superior prediction accuracy for all positions. The simulated FM was used to optimize B0 shims with up to 22.2% improvement with respect to the transformed FM approach. CONCLUSION: The proposed simulation strategy is able to predict movement-induced B0 field inhomogeneities yielding more precise estimates of the ground truth field homogeneity than the transformed FM.

Correcting inter-scan motion artifacts in quantitative R1 mapping at 7T.

PURPOSE: Inter-scan motion is a substantial source of error in R1 estimation methods based on multiple volumes, for example, variable flip angle (VFA), and can be expected to increase at 7T where B1 fields are more inhomogeneous. The established correction scheme does not translate to 7T since it requires a body coil reference. Here we introduce two alternatives that outperform the established method. Since they compute relative sensitivities they do not require body coil images. THEORY: The proposed methods use coil-combined magnitude images to obtain the relative coil sensitivities. The first method efficiently computes the relative sensitivities via a simple ratio; the second by fitting a more sophisticated generative model. METHODS: R1 maps were computed using the VFA approach. Multiple datasets were acquired at 3T and 7T, with and without motion between the acquisition of the VFA volumes. R1 maps were constructed without correction, with the proposed corrections, and (at 3T) with the previously established correction scheme. The effect of the greater inhomogeneity in the transmit field at 7T was also explored by acquiring B1+ maps at each position. RESULTS: At 3T, the proposed methods outperform the baseline method. Inter-scan motion artifacts were also reduced at 7T. However, at 7T reproducibility only converged on that of the no motion condition if position-specific transmit field effects were also incorporated. CONCLUSION: The proposed methods simplify inter-scan motion correction of R1 maps and are applicable at both 3T and 7T, where a body coil is typically not available. The open-source code for all methods is made publicly available.

Developing formalin-based fixative agents for post mortem brain MRI at 9.4 T.

PURPOSE: To develop fixative agents for high-field MRI with suitable dielectric properties and measure MR properties in immersion-fixed brain tissue. METHODS: Dielectric properties of formalin-based agents were assessed (100 MHz-4.5 GHz), and four candidate fixatives with/without polyvinylpyrrolidone (PVP) and different salt concentrations were formulated. B1 field and MR properties (T1 , R2∗ , R2 , R2' , and magnetic susceptibility [QSM]) were observed in white and gray matter of pig brain samples during 0.5-35 days of immersion fixation. The kinetics were fitted using exponential functions. The immersion time required to reach maximum R2∗ values at different tissue depths was used to estimate the Medawar coefficient for fixative penetration. The effect of replacing the fixatives with Fluoroinert and phosphate-buffered saline as embedding media was also evaluated. RESULTS: The dielectric properties of formalin were nonlinearly modified by increasing amounts of additives. With 5% PVP and 0.04% NaCl, the dielectric properties and B1 field reflected in vivo conditions. The highest B1 values were found in white matter with PVP and varied significantly with tissue depth and embedding media, but not with immersion time. The MR properties depended on PVP yielding lower T1 , higher R2∗ , more paramagnetic QSM values, and a lower Medawar coefficient (0.9 mm/h ; without PVP: 1.5). Regardless of fixative, switching to phosphate-buffered saline as embedder caused a paramagnetic shift in QSM and decreased R2∗ that progressed during 1 month of storage, whereas no differences were found with Fluorinert. CONCLUSION: In vivo-like B1 fields can be achieved in formalin fixatives using PVP and a low salt concentration, yielding lower T1 , higher R2∗ , and more paramagnetic QSM than without additives. The kinetics of R2∗ allowed estimation of fixative tissue penetration.

Transcranial direct current stimulation with functional magnetic resonance imaging: a detailed validation and operational guide.

Introduction: Transcranial direct current stimulation (tDCS) is a non-invasive brain stimulation technique used to modulate human brain and behavioural function in both research and clinical interventions. The combination of functional magnetic resonance imaging (fMRI) with tDCS enables researchers to directly test causal contributions of stimulated brain regions, answering questions about the physiology and neural mechanisms underlying behaviour. Despite the promise of the technique, advances have been hampered by technical challenges and methodological variability between studies, confounding comparability/replicability. Methods: Here tDCS-fMRI at 3T was developed for a series of experiments investigating language recovery after stroke. To validate the method, one healthy volunteer completed an fMRI paradigm with three conditions: (i) No-tDCS, (ii) Sham-tDCS, (iii) 2mA Anodal-tDCS. MR data were analysed in SPM12 with region-of-interest (ROI) analyses of the two electrodes and reference sites. Results: Quality assessment indicated no visible signal dropouts or distortions introduced by the tDCS equipment. After modelling scanner drift, motion-related variance, and temporal autocorrelation, we found no field inhomogeneity in functional sensitivity metrics across conditions in grey matter and in the three ROIs. Discussion: Key safety factors and risk mitigation strategies that must be taken into consideration when integrating tDCS into an fMRI environment are outlined. To obtain reliable results, we provide practical solutions to technical challenges and complications of the method. It is hoped that sharing these data and SOP will promote methodological replication in future studies, enhancing the quality of tDCS-fMRI application, and improve the reliability of scientific results in this field. Conclusions: The method and data provided here provide a technically safe, reliable tDCS-fMRI procedure to obtain high quality MR data. The detailed framework of the Standard Operation Procedure SOP ( https://doi.org/10.5281/zenodo.4606564) systematically reports the technical and procedural elements of our tDCS-fMRI approach, which we hope can be adopted and prove useful in future studies.

A 32-channel multi-coil setup optimized for human brain shimming at 9.4T.

PURPOSE: A multi-coil shim setup is designed and optimized for human brain shimming. Here, the size and position of a set of square coils are optimized to improve the shim performance without increasing the number of local coils. Utilizing such a setup is especially beneficial at ultrahigh fields where B0 inhomogeneity in the human brain is more severe. METHODS: The optimization started with a symmetric arrangement of 32 independent coils. Three parameters per coil were optimized in parallel, including angular and axial positions on a cylinder surface and size of the coil, which were constrained by cylinder size, construction consideration, and amplifiers specifications. B0 maps were acquired at 9.4T in 8 healthy volunteers for use as training data. The global and dynamic shimming performance of the optimized multi-coil were compared in simulations and measurements to a symmetric design and to the scanner's second-order shim setup, respectively. RESULTS: The optimized multi-coil performs better by 14.7% based on standard deviation (SD) improvement with constrained global shimming in comparison to the symmetric positioning of the coils. Global shimming performance was comparable with a symmetric 65-channel multi-coil and full fifth-order spherical harmonic shim coils. On average, an SD of 48.4 and 31.9 Hz was achieved for in vivo measurements after global and dynamic slice-wise shimming, respectively. CONCLUSIONS: An optimized multi-coil shim setup was designed and constructed for human whole-brain shimming. Similar performance of the multi-coils with many channels can be achieved with a fewer number of channels when the coils are optimally arranged around the target.

Impact of prospective motion correction, distortion correction methods and large vein bias on the spatial accuracy of cortical laminar fMRI at 9.4 Tesla.

Functional imaging with sub-millimeter spatial resolution is a basic requirement for assessing functional MRI (fMRI) responses across different cortical depths and is used extensively in the emerging field of laminar fMRI. Such studies seek to investigate the detailed functional organization of the brain and may develop to a new powerful tool for human neuroscience. However, several studies have shown that measurement of laminar fMRI responses can be biased by the image acquisition and data processing strategies. In this work, measurements with three different gradient-echo EPI BOLD fMRI protocols with a voxel size down to 650 â€‹µm isotropic were performed at 9.4 â€‹T. We estimated how prospective motion correction can help to improve spatial accuracy by reducing the number of spatial resampling steps in postprocessing. In addition, we demonstrate key requirements for accurate geometric distortion correction to ensure that distortion correction maps are properly aligned to the functional data and that strong variations of distortions near large veins can lead to signal overlays which cannot be corrected for during postprocessing. Furthermore, this study illustrates the spatial extent of bias induced by pial and other larger veins in laminar BOLD experiments. Since these issues under investigation affect studies performed with more conventional spatial resolutions, the methods applied in this work may also help to improve the understanding of the BOLD signal more broadly.

Design of a shim coil array matched to the human brain anatomy.

PURPOSE: The purpose of this study is to introduce a novel design method of a shim coil array specifically optimized for whole brain shimming and to compare the performance of the resulting coils to conventional spherical harmonic shimming. METHODS: The proposed design approach is based on the stream function method and singular value decomposition. Eighty-four field maps from 12 volunteers measured in seven different head positions were used during the design process. The cross validation technique was applied to find an optimal number of coil elements in the array. Additional 42 field maps from 6 further volunteers were used for an independent validation. A bootstrapping technique was used to estimate the required population size to achieve a stable coil design. RESULTS: Shimming using 12 and 24 coil elements outperforms fourth- and fifth-order spherical harmonic shimming for all measured field maps, respectively. Coil elements show novel coil layouts compared to the conventional spherical harmonic coils and existing multi-coils. Both leave-one-out and independent validation demonstrate the generalization ability of the designed arrays. The bootstrapping analysis predicts that field maps from approximately 140 subjects need to be acquired to arrive at a stable design. CONCLUSIONS: The results demonstrate the validity of the proposed method to design a shim coil array matched to the human brain anatomy, which naturally satisfies the laws of electrodynamics. The design method may also be applied to develop new shim coil arrays matched to other human organs.

Dynamic B0 shimming of the motor cortex and cerebellum with a multicoil shim setup for BOLD fMRI at 9.4T.

PURPOSE: We assessed how improved static magnetic field (B0 ) homogeneity with a dynamic multicoil shimming can influence the blood oxygen level dependent (BOLD) contrast to noise when echo planar imaging (EPI) sequence is used for a motor task functional MRI study. We showed that a multicoil shim setup can be a proper choice for dynamic shimming of 2 spatially distant areas with different inhomogeneity distributions. METHODS: A 16-channel multicoil shim setup is used to provide improved B0 homogeneity by dynamic slice-wise shimming. The performance of dynamic B0 shimming was investigated in 2 distinct brain regions, the motor cortex and the cerebellum, in the same experiment during a finger-tapping task. Temporal SNR (tSNR), geometric distortion of the EPIs, and results of an analysis with a general linear model before and after shimming with the multicoil were compared. RESULTS: Reduced B0 deviation by 30% and 52% in the cerebellum and motor cortex, respectively, resulted in higher tSNR and a reduction of distortions in the EPI. Statistical analysis applied to the EPIs showed higher t values and increased number of voxels above significance threshold when shimming with the multicoil setup. CONCLUSIONS: Improved B0 homogeneity leads to higher tSNR and enhances the detection of BOLD signal.

Spread-spectrum magnetic resonance imaging.

PURPOSE: A novel method for the acceleration of MRI acquisition is proposed that relies on the local modulation of magnetic fields. These local modulations provide additional spatial information for image reconstruction that is used to accelerate image acquisition. METHODS: In experiments and simulations, eight local coils connected to current amplifiers were used for rapid local magnetic field variation. Acquired and simulated data were reconstructed to quantify reconstruction errors as a function of the acceleration factor and applied modulation frequency and strength. RESULTS: Experimental results demonstrate a possible acceleration factor of 2 to 4. Simulations demonstrate the challenges and limits of this method in terms of required magnetic field modulation strengths and frequencies. A normalized mean squared error of below 10% can be achieved for acceleration factors of up to 8 using modulation field strengths comparable to the readout gradient strength at modulation frequencies in the range of 5 to 20 kHz. CONCLUSION: Spread-spectrum MRI represents a new approach to accelerate image acquisition, and it can be independently combined with traditional parallel imaging techniques based on local receive coil sensitivities.

Comparison of prospective head motion correction with NMR field probes and an optical tracking system.

PURPOSE: The aim of this study was to compare prospective head motion correction and motion tracking abilities of two tracking systems: Active NMR field probes and a Moiré phase tracking camera system using an optical marker. METHODS: Both tracking systems were used simultaneously on human subjects. The prospective head motion correction was compared in an MP2RAGE and a gradient echo sequence. In addition, the motion tracking trajectories for three subjects were compared against each other and their correlation and deviations were analyzed. RESULTS: With both tracking systems motion artifacts were visibly reduced. The precision of the field probe system was on the order of 50 µm for translations and 0.03° for rotations while the camera's was approximately 5 µm and 0.007°. The comparison of the measured trajectories showed close correlation and an average absolute deviation below 500 µm and 0.5°. CONCLUSION: This study presents the first in vivo comparison between NMR field probes and Moiré phase tracking. For the gradient echo images, the field probes had a similar motion correction performance as the optical tracking system. For the MP2RAGE measurement, however, the camera yielded better results. Still, both tracking systems substantially decreased image artifacts in the presence of subject motion. Thus, the motion tracking modality should be chosen according to the specific requirements of the experiment while considering the desired image resolution, refresh rate, and head coil constraints.

Dynamic B0 shimming of the human brain at 9.4 T with a 16-channel multi-coil shim setup.

PURPOSE: A 16-channel multi-coil shimming setup was developed to mitigate severe B0 field perturbations at ultrahigh field and improve data quality for human brain imaging and spectroscopy. METHODS: The shimming setup consisted of 16 circular B0 coils that were positioned symmetrically on a cylinder with a diameter of 370 mm. The latter was large enough to house a shielded 18/32-channel RF transceiver array. The shim performance was assessed via simulations and phantom as well as in vivo measurements at 9.4 T. The global and dynamic shimming performance of the multi-coil setup was compared with the built-in scanner shim system for EPI and single voxel spectroscopy. RESULTS: The presence of the multi-coil shim did not influence the performance of the RF coil. The performance of the proposed setup was similar to a full third-order spherical harmonic shim system in the case of global static and dynamic slice-wise shimming. Dynamic slice-wise shimming with the multi-coil setup outperformed global static shimming with the scanner's second-order spherical-harmonic shim. The multi-coil setup allowed mitigating geometric distortions for EPI. The combination of the multi-coil shim setup with the zeroth and first-order shim of the scanner further reduced the standard deviation of the B0 field in the brain by 12% compared with the case in which multi-coil was used exclusively. CONCLUSION: The combination of a multi-coil setup and the linear shim channels of the scanner provides a straightforward solution for implementing dynamic slice-wise shimming without requiring an additional pre-emphasis setup.