An automatic L1-based regularization method for the analysis of FFC dispersion profiles with quadrupolar peaks

Fast Field-Cycling Nuclear Magnetic Resonance relaxometry is a non-destructive technique to investigate molecular dynamics and structure of systems having a wide range of applications such as environment, biology, and food. Besides a considerable amount of literature about modeling and application of such technique in specific areas, an algorithmic approach to the related parameter identification problem is still lacking. We believe that a robust algorithmic approach will allow a unified treatment of different samples in several application areas. In this paper, we model the parameters identification problem as a constrained $L_1$-regularized non-linear least squares problem. Following the approach proposed in [Analytical Chemistry 2021 93 (24)], the non-linear least squares term imposes data consistency by decomposing the acquired relaxation profiles into relaxation contributions associated with 1H-1H and 1H-14N dipole-dipole interactions. The data fitting and the L1-based regularization terms are balanced by the so-called regularization parameter. For the parameters identification, we propose an algorithm that computes, at each iteration, both the regularization parameter and the model parameters. In particular, the regularization parameter value is updated according to a Balancing Principle and the model parameters values are obtained by solving the corresponding $L_1$-regularized non-linear least squares problem by means of the non-linear Gauss-Seidel method. We analyse the convergence properties of the proposed algorithm and run extensive testing on synthetic and real data. A Matlab software, implementing the presented algorithm, is available upon request to the authors.