Numerical solutions to the time-dependent Bloch equations revisited.

Area of Science:

• Magnetic Resonance Imaging (MRI)
• Biophysics
• Computational Methods

Background:

• The Bloch equations are fundamental to understanding nuclear magnetic resonance (NMR) phenomena.
• Solving time-dependent Bloch equations is computationally intensive, especially for complex models like chemical exchange saturation transfer (CEST) and amide proton transfer (APT).
• Existing numerical methods, such as the Runge-Kutta-Fehlberg (RKF) method, can be time-consuming.

Purpose of the Study:

• To present a novel, simple, and rapid method for solving the time-dependent Bloch equations.
• To validate the proposed method against analytical solutions and established numerical techniques.
• To apply the method to analyze complex MRI contrast mechanisms like CEST and APT.

Main Methods:

• Reduced the time-dependent Bloch equations to a homogeneous linear differential equation.
• Derived a simple matrix operation for solving the differential equation.
• Compared results with analytical solutions for constant radiofrequency irradiation.
• Applied the method to the two-pool exchange model for CEST/APT MRI and compared with the RKF method.

Main Results:

• The matrix-based method demonstrated good agreement with analytical solutions.
• Calculated Z-spectra and asymmetry spectra for CEST/APT MRI using the new method showed excellent agreement with RKF results.
• The proposed method was significantly faster than the RKF method.

Conclusions:

• The developed matrix-based method is a valid and efficient approach for solving time-dependent Bloch equations.
• This method offers a substantial speed advantage over the RKF method for MRI applications.
• It will be valuable for analyzing CEST/APT contrast mechanisms and optimizing CEST/APT MRI protocols.