Dataset Information

Interslice current change constrained B₀ shim optimization for accurate high-order dynamic shim updating with strongly reduced eddy currents.

ABSTRACT:

Purpose

To overcome existing challenges in dynamic B₀ shimming by implementing a shim optimization algorithm which limits shim current amplitudes and their temporal variation through the application of constraints and regularization terms.

Theory and methods

Spherical harmonic dynamic B₀ shimming is complicated by eddy currents, ill-posed optimizations, and the need for strong power supplies. Based on the fact that eddy current amplitudes are proportional to the magnitude of the shim current changes, and assuming a smoothness of the B₀ inhomogeneity variation in the slice direction, a novel algorithm was implemented to reduce eddy current generation by limiting interslice shim current changes. Shim degeneracy issues and resulting high current amplitudes are additionally addressed by penalizing high solution norms. Applicability of the proposed algorithm was validated in simulations and in phantom and in vivo measurements.

Results

High-order dynamic shimming simulations and measurements have shown that absolute shim current amplitudes and their temporal variation can be substantially reduced with negligible loss in achievable B₀ homogeneity. Whereas conventional dynamic shim updating optimizations improve the B₀ homogeneity, on average, by a factor of 2.1 over second-order static solutions, our proposed routine reached a factor of 2.0, while simultaneously providing a 14-fold reduction of the average maximum shim current changes.

Conclusions

The proposed algorithm substantially reduces the shim amplitudes and their temporal variation, while only marginally affecting the achievable B₀ homogeneity. As a result, it has the potential to mitigate the remaining challenges in dynamic B₀ shimming and help in making its application more readily available.

SUBMITTER: Schwerter M

PROVIDER: S-EPMC8168441 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:PurposeThe purpose of this study was to measure gantry angle-related eddy currents in a 0.35-T MRI-Linac and determine if B0 (zeroth order) eddy currents are the primary cause of gantry angle-dependent imaging isocenter shifts vs other potential causes like B0 inhomogeneities and gradient (first order) eddy currents. For conventional Cartesian acquisitions, B0 eddy currents can cause imaging isocenter shifts along both phase encode and readout directions. Gradient eddy currents can cause spatial distortion along both the phase encode and readout directions. Center frequency offsets can cause imaging isocenter shifts along the readout direction that vary with readout gradient polarity.MethodsMRI-related eddy currents and imaging isocenter shifts were measured on a 0.35-T MRI-Linac at gantry angles from 0° to 330° in increments of 30° . All measurements were made after gradient shimming and center frequency tuning at each planned gantry angle. Eddy current and field homogeneity measurements were conducted using a 24-cm diameter spherical phantom. Gradient and B0 eddy currents were calculated from the free induction decays (FIDs) resulting from selective excitation of slices located ±5 cm from isocenter. B0 eddy currents were also calculated from FIDs acquired with nonselective excitation and compared with B0 eddy current values derived using selective excitation. B0 inhomogeneities and center frequency offsets were measured by acquiring FIDs with nonselective excitation. Imaging isocenter shifts were measured using a 33x33x10.5 cm3 uniformity linearity (grid) phantom and a 3D true fast imaging with steady-state precession (TrueFISP) sequence used in MRI-guided radiation therapy. Eddy currents were compared to vendor specifications and correlated with the imaging isocenter shifts. Measurements were conducted before and after the MRI-Linac's waveguide was replaced with an updated design to reduce eddy currents.ResultsB0 eddy currents were highly correlated (r = 0.986, P << 0.001) for measurements made with vs without selective excitation. Transverse (X and Y) axis B0 eddy currents before and after the waveguide upgrade were out of specification (specification: ≤0.1 μT m/mT for delays < 10 ms) for most of the measured gantry angles. Gradient eddy currents before and after the upgrade were within specifications for the measured gantry angles (≤0.1% for delays < 10 ms). B0 eddy currents and imaging isocenter shifts were highly correlated (r = 0.965, P << 0.001). After the Linac waveguide upgrade, root mean square (RMS) peak B0 and gradient eddy currents dropped 45% and 11%, respectively, for delays <10 ms, while imaging isocenter shifts dropped 53%. Isocenter shifts were observed in both phase encode and readout directions. Center frequency offsets were <26 Hz while B0 inhomogeneities were <33 Hz full width at half maximum (FWHM).ConclusionsImaging isocenter shifts measured in a 0.35-T MRI-Linac were highly correlated with B0 eddy currents. The eddy currents and imaging isocenter shifts decreased after the MRI-Linac's waveguide was replaced.

Project description:We derive a low-scaling G0W0 algorithm for molecules using pair atomic density fitting (PADF) and an imaginary time representation of the Green's function and describe its implementation in the Slater type orbital (STO)-based Amsterdam density functional (ADF) electronic structure code. We demonstrate the scalability of our algorithm on a series of water clusters with up to 432 atoms and 7776 basis functions and observe asymptotic quadratic scaling with realistic threshold qualities controlling distance effects and basis sets of triple-ζ (TZ) plus double polarization quality. Also owing to a very small prefactor, a G0W0 calculation for the largest of these clusters takes only 240 CPU hours with these settings. We assess the accuracy of our algorithm for HOMO and LUMO energies in the GW100 database. With errors of 0.24 eV for HOMO energies on the quadruple-ζ level, our implementation is less accurate than canonical all-electron implementations using the larger def2-QZVP GTO-type basis set. Apart from basis set errors, this is related to the well-known shortcomings of the GW space-time method using analytical continuation techniques as well as to numerical issues of the PADF approach of accurately representing diffuse atomic orbital (AO) products. We speculate that these difficulties might be overcome by using optimized auxiliary fit sets with more diffuse functions of higher angular momenta. Despite these shortcomings, for subsets of medium and large molecules from the GW5000 database, the error of our approach using basis sets of TZ and augmented double-ζ (DZ) quality is decreasing with system size. On the augmented DZ level, we reproduce canonical, complete basis set limit extrapolated reference values with an accuracy of 80 meV on average for a set of 20 large organic molecules. We anticipate our algorithm, in its current form, to be very useful in the study of single-particle properties of large organic systems such as chromophores and acceptor molecules.