Model-Based Reconstruction for MPI

A prerequisite for system matrix-based image reconstruction in Magnetic Particle Imaging (MPI) is the acquisition of a system matrix that describes the mapping between the MPI tracer and the measured signal. A common method for its acquisition is a time-consuming calibration procedure, during which the scanner is blocked for other uses. For this reason, a model-based approach that allows the system matrix to be obtained by simulation has great appeal. However, an accurate model that describes the magnetization behavior of the tracer, allows the identification of its parameters, and is computationally feasible has not yet been found.

In an ongoing collaboration with Hannes Albers and Tobias Kluth from the University of Bremen we investigate and refine a magnetization model based on the Néel rotation for the magnetic moments of the particles (see Kluth et al., 2019 and Albers et al., 2022). On the one hand, the identification of the parameters of the model is in focus, since these are unknown a priori, on the other hand, measured 2D MPI system matrices are describe with much higher accuracy than the current MPI models. Moreover, we are also interested in the limitations current MPI models in the context of fluid dynamics (see Möddel et al., 2023).

A comparison of the frequency components between a measured and model-based system matrix shows significant differences, as the simple and widely used equilibrium model is not able to fully capture the complex magnetization dynamics.

Project Publications

[164777]
Title: Modeling the magnetization dynamics for large ensembles of immobilized magnetic nanoparticles in multi-dimensional magnetic particle imaging.
Written by: H. Albers, T. Knopp, M. Möddel, M. Boberg, and T. Kluth
in: <em>Journal of Magnetism and Magnetic Materials</em>. February (2022).
Volume: <strong>543</strong>. Number:
on pages: 168534
Chapter:
Editor:
Publisher:
Series:
Address:
Edition:
ISBN:
how published:
Organization:
School:
Institution:
Type:
DOI: 10.1016/j.jmmm.2021.168534
URL: https://arxiv.org/abs/2106.08040
ARXIVID:
PMID:

[www]

Note: article, model-based

Abstract: Magnetic nanoparticles (MNPs) play an important role in biomedical applications including imaging modalities such as magnetic resonance imaging (MRI) and magnetic particle imaging (MPI). The latter one exploits the non-linear magnetization response of a large ensemble of magnetic nanoparticles to magnetic fields which allows determining the spatial distribution of the MNP concentration from measured voltage signals. The image-to-voltage mapping is linear and described by a system matrix. Currently, modeling the voltage signals of large ensembles of MNPs in an MPI environment is not yet accurately possible, especially for liquid tracers in multi-dimensional magnetic excitation fields. As an immediate consequence, the system matrix is still obtained in a time consuming calibration procedure. While the ferrofluidic case can be seen as the typical setting, more recently immobilized and potentially oriented MNPs have received considerable attention. By aligning the particles magnetic easy axis during immobilization one can encode the angle of the particle’s magnetic easy axis into the magnetization response providing a relevant benchmark system for model-based approaches. In this work we address the modeling problem for immobilized and oriented MNPs in the context of MPI. We investigate a model-based approach where the magnetization response is simulated by a Néel rotation model for the particle’s magnetic moments and the ensemble magnetization is obtained by solving a Fokker–Planck equation approach. Since the parameters of the model are a-priori unknown, we investigate different methods for performing a parameter identification and discuss two different models: One where a single function vector is used from the space spanned by the model parameters and another where a superposition of function vectors is considered. We show that our model can much more accurately reproduce the orientation dependent signal response when compared to the equilibrium model, which marks the current state-of-the-art for model-based system matrix simulations in MPI.