Enhanced and robust contrast in CEST MRI: Saturation pulse shape design via optimal control.

PURPOSE: To employ optimal control for the numerical design of Chemical Exchange Saturation Transfer (CEST) saturation pulses to maximize contrast and stability against B 0 $$ {\mathrm{B}}_0 $$ inhomogeneities. THEORY AND METHODS: We applied an optimal control framework for the design pulse shapes for CEST saturation pulse trains. The cost functional minimized both the pulse energy and the discrepancy between the corresponding CEST spectrum and the target spectrum based on a continuous radiofrequency (RF) pulse. The optimization is subject to hardware limitations. In measurements on a 7 T preclinical scanner, the optimal control pulses were compared to continuous-wave and Gaussian saturation methods. We conducted a comparison of the optimal control pulses with Gaussian, block pulse trains, and adiabatic spin-lock pulses. RESULTS: The optimal control pulse train demonstrated saturation levels comparable to continuous-wave saturation and surpassed Gaussian saturation by up to 50 % in phantom measurements. In phantom measurements at 3 T the optimized pulses not only showcased the highest CEST contrast, but also the highest stability against field inhomogeneities. In contrast, block pulse saturation resulted in severe artifacts. Dynamic Bloch-McConnell simulations were employed to identify the source of these artifacts, and underscore the B 0 $$ {\mathrm{B}}_0 $$ robustness of the optimized pulses. CONCLUSION: In this work, it was shown that a substantial improvement in pulsed saturation CEST imaging can be achieved by using Optimal Control design principles. It is possible to overcome the sensitivity of saturation to B0 inhomogeneities while achieving CEST contrast close to continuous wave saturation.

Robust dual-angle T 1 $$ {T}_1 $$ measurement in magnetization transfer spectroscopy by time-optimal control.

Magnetization transfer spectroscopy relies heavily on the robust determination of T 1 $$ {T}_1 $$ relaxation times of nuclei participating in metabolic exchange. Challenges arise due to the use of surface RF coils for transmission (high B 1 + $$ {B}_1^{+} $$ variation) and the broad resonance band of most X nuclei. These challenges are particularly pronounced when fast T 1 $$ {T}_1 $$ mapping methods, such as the dual-angle method, are employed. Consequently, in this work, we develop resonance offset and B 1 + $$ {B}_1^{+} $$ robust excitation RF pulses for 31P magnetization transfer spectroscopy at 7T through ensemble-based time-optimal control. In our approach, we introduce a cost functional for designing robust pulses, incorporating the full Bloch equations as constraints, which are solved using symmetric operator splitting techniques. The optimal control design of the RF pulses developed demonstrates improved accuracy, desired phase properties, and reduced RF power when applied to dual-angle T 1 $$ {T}_1 $$ mapping, thereby improving the precision of exchange-rate measurements, as demonstrated in a preclinical in vivo study quantifying brain creatine kinase activity.

Advanced design of MRI inversion pulses for inhomogeneous field conditions by optimal control.

Non-selective inversion pulses find widespread use in MRI applications, where requirements on them are increasingly demanding. With the use of high and ultra-high field strength systems, robustness to Δ B 0 and B 1 + inhomogeneities, while tackling SAR and hardware limitations, has rapidly become important. In this work, we propose a time-optimal control framework for the optimization of Δ B 0 - and B 1 + -robust inversion pulses. Robustness is addressed by means of ensemble formulations, while allowing inclusion of hardware and energy limitations. The framework is flexible and performs excellently for various optimization goals. The optimization results are analyzed extensively in numerical experiments. Furthermore, they are validated, and compared with adiabatic RF pulses, in various phantom and in vivo measurements on a 3 T MRI system.

Accuracy and performance analysis for Bloch and Bloch-McConnell simulation methods.

PURPOSE: To introduce new solution methods for the Bloch and Bloch-McConnell equations and compare them quantitatively to different known approaches. THEORY AND METHODS: A new exact solution per time step is derived by means of eigenvalues and generalized eigenvectors. Fast numerical solution methods based on asymmetric and symmetric operator splitting, which are already known for the Bloch equations, are extended to the Bloch-McConnell equations. Those methods are compared to other numerical methods including spin domain, one-step and multi-step methods, and matrix exponential. Error metrics are introduced based on the exact solution method, which allows to assess the accuracy of each solution method quantitatively for arbitrary example data. RESULTS: Accuracy and performance properties for nine different solution methods are analyzed and compared in extensive numerical experiments including various examples for non-selective and slice-selective MR imaging applications. The accuracy of the methods heavily varies, in particular for short relaxation times and long pulse durations. CONCLUSION: In absence of relaxation effects, the numerical results confirm the rotation matrices approach as accurate and computationally efficient Bloch solution method. Otherwise, as well as for the Bloch-McConnell equations, symmetric operator splitting methods are recommended due to their excellent numerical accuracy paired with efficient run time.

Personalization of electro-mechanical models of the pressure-overloaded left ventricle: fitting of Windkessel-type afterload models.

Computer models of left ventricular (LV) electro-mechanics (EM) show promise as a tool for assessing the impact of increased afterload upon LV performance. However, the identification of unique afterload model parameters and the personalization of EM LV models remains challenging due to significant clinical input uncertainties. Here, we personalized a virtual cohort of N = 17 EM LV models under pressure overload conditions. A global-local optimizer was developed to uniquely identify parameters of a three-element Windkessel (Wk3) afterload model. The sensitivity of Wk3 parameters to input uncertainty and of the EM LV model to Wk3 parameter uncertainty was analysed. The optimizer uniquely identified Wk3 parameters, and outputs of the personalized EM LV models showed close agreement with clinical data in all cases. Sensitivity analysis revealed a strong dependence of Wk3 parameters on input uncertainty. However, this had limited impact on outputs of EM LV models. A unique identification of Wk3 parameters from clinical data appears feasible, but it is sensitive to input uncertainty, thus depending on accurate invasive measurements. By contrast, the EM LV model outputs were less sensitive, with errors of less than 8.14% for input data errors of 10%, which is within the bounds of clinical data uncertainty. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Time optimal control-based RF pulse design under gradient imperfections.

PURPOSE: This study incorporates a gradient system imperfection model into an optimal control framework for radio frequency (RF) pulse design. THEORY AND METHODS: The joint design of minimum-time RF and slice selective gradient shapes is posed as an optimal control problem. Hardware limitations such as maximal amplitudes for RF and slice selective gradient or its slew rate are included as hard constraints to assure practical applicability of the optimized waveforms. In order to guarantee the performance of the optimized waveform with possible gradient system disturbances such as limited system bandwidth and eddy currents, a measured gradient impulse response function (GIRF) for a specific system is integrated into the optimization. RESULTS: The method generates optimized RF and pre-distorted slice selective gradient shapes for refocusing that are able to fully compensate the modeled imperfections of the gradient system under investigation. The results nearly regenerate the optimal results of an idealized gradient system. The numerical Bloch simulations are validated by phantom and in-vivo experiments on 2 3T scanners. CONCLUSIONS: The presented design approach demonstrates the successful correction of gradient system imperfections within an optimal control framework for RF pulse design.

A time domain signal equation for multi-echo spin-echo sequences with arbitrary excitation and refocusing angle and phase.

Accurate T2 mapping using multi-echo spin-echo is usually impaired by non-ideal refocusing due to B1+ inhomogeneities and slice profile effects. Incomplete refocusing gives rise to stimulated echo and so called "T1-mixing" and consequently a non-exponential signal decay. Here we present a time domain formula that incorporates all relaxation and pulse parameters and enables accurate and realistic modelling of the magnetization decay curve. By pulse parameters here we specifically mean the actual refocusing angle and axis, and phase angle of both the excitation and refocusing pulse. The method used for derivation comprises the so called Generating functions approach with subsequent back-transformation to the time domain. The proposed approach was validated by simulations using realistic RF pulse shapes as well as by comparison to phantom measurements. Excellent agreement between simulations and measurements underpin the validity of the presented approach. Conclusively, we here present a complete time domain formula ready to use for accurate T2 mapping with multi-echo spin-echo sequences.

Simultaneous multislice refocusing via time optimal control.

PURPOSE: Joint design of minimum duration RF pulses and slice-selective gradient shapes for MRI via time optimal control with strict physical constraints, and its application to simultaneous multislice imaging. THEORY AND METHODS: The minimization of the pulse duration is cast as a time optimal control problem with inequality constraints describing the refocusing quality and physical constraints. It is solved with a bilevel method, where the pulse length is minimized in the upper level, and the constraints are satisfied in the lower level. To address the inherent nonconvexity of the optimization problem, the upper level is enhanced with new heuristics for finding a near global optimizer based on a second optimization problem. RESULTS: A large set of optimized examples shows an average temporal reduction of 87.1% for double diffusion and 74% for turbo spin echo pulses compared to power independent number of slices pulses. The optimized results are validated on a 3T scanner with phantom measurements. CONCLUSION: The presented design method computes minimum duration RF pulse and slice-selective gradient shapes subject to physical constraints. The shorter pulse duration can be used to decrease the effective echo time in existing echo-planar imaging or echo spacing in turbo spin echo sequences.

Magnetic Resonance RF Pulse Design by Optimal Control With Physical Constraints.

Optimal control approaches have proved useful in designing RF pulses for large tip-angle applications. A typical challenge for optimal control design is the inclusion of constraints resulting from physiological or technical limitations that assure the realizability of the optimized pulses. In this paper, we show how to treat such inequality constraints, in particular, amplitude constraints on the B1 field, the slice-selective gradient, and its slew rate, as well as constraints on the slice profile accuracy. For the latter, a pointwise profile error and additional phase constraints are prescribed. Here, a penalization method is introduced that corresponds to a higher order tracking instead of the common quadratic tracking. The order is driven to infinity in the course of the optimization. We jointly optimize for the RF and slice-selective gradient waveform. The amplitude constraints on these control variables are treated efficiently by semismooth Newton or quasi-Newton methods. The method is flexible, adapting to many optimization goals. As an application, we reduce the power of refocusing pulses, which is important for spin echo-based applications with a short echo spacing. Here, the optimization method is tested in numerical experiments for reducing the pulse power of simultaneous multislice refocusing pulses. The results are validated by phantom and in-vivo experiments.

Efficient high-resolution RF pulse design applied to simultaneous multi-slice excitation.

RF pulse design via optimal control is typically based on gradient and quasi-Newton approaches and therefore suffers from slow convergence. We present a flexible and highly efficient method that uses exact second-order information within a globally convergent trust-region CG-Newton method to yield an improved convergence rate. The approach is applied to the design of RF pulses for single- and simultaneous multi-slice (SMS) excitation and validated using phantom and in vivo experiments on a 3T scanner using a modified gradient echo sequence.