Over recent years, compact and low-field time-domain nuclear magnetic resonance (TD-NMR) instruments have become increasingly available, expanding their use in the characterization of biomaterials across food, plant, and agro-industrial research. In this context, the inverse Laplace transform (ILT) has emerged as a powerful mathematical approach for extracting relaxation time distributions from TD-NMR signals. However, despite its widespread use, ILT analysis is often restricted to proprietary software or requires advanced expertise in numerical methods, limiting its accessibility to non-specialist users. In this work, we present an open-access WebApp for performing ILT analysis of TD-NMR signals in a transparent and user-friendly manner. The implemented algorithm is based on non-negative least squares combined with Tikhonov regularization and singular value decomposition, allowing robust inversion of ill-posed relaxation data. The platform supports the main TD-NMR experiments used in practice, including Carr–Purcell–Meiboom–Gill (CPMG), inversion recovery, and saturation recovery pulse sequences, and is compatible with data from instruments of any manufacturer. In addition to describing the mathematical formulation and implementation of the algorithm, a concise methodological discussion of ILT in the context of TD-NMR is provided. The performance of the WebApp is evaluated using both simulated datasets and representative experimental signals, demonstrating that the obtained relaxation time distributions are consistent with those produced by established ILT approaches. By lowering the barrier to advanced signal processing, the proposed WebApp represents a useful open scientific tool for research and teaching in magnetic resonance applications.
Conselho Nacional de Desenvolvimento Científico e Tecnológico306611/2022-8Fundação de Amparo à Pesquisa do Estado de São Paulo2025/09284-9Introduction
The two best-known nuclear magnetic resonance (NMR) techniques in academia and industry are magnetic resonance imaging (MRI) and high-resolution NMR (HR-NMR). The former is widely applied in radiology for non-invasive medical diagnostics, whereas the latter is a well-established analytical technique used in chemical analysis laboratories for the precise characterization of molecular structures, playing a fundamental role in the pharmaceutical industry .
Time-domain NMR (TD-NMR), also referred to as low-field NMR (LF-NMR) or low-resolution NMR (LR-NMR) , is based on the same physical principles as MRI and HR-NMR but employs simpler instrumentation based on permanent magnets, without the need for large superconducting systems or cryogenic liquids. As a result, TD-NMR sensors and spectrometers have reduced size and cost, enabling the development of portable and benchtop instruments. These characteristics make TD-NMR a versatile technique for industrial and applied research, where it has increasingly been employed as a powerful analytical sensor . Recent advances in electronics and permanent magnet technology have further contributed to the reduction of production and maintenance costs, fostering the emergence of several companies dedicated to compact TD-NMR instrumentation, such as Bruker, Magritek, Oxford, Nanalysis, SpinLock, FIT, and Niumag .
For materials of biological origin, such as plants, foods, oils, seeds, and other agro-industrial products, magnetic resonance offers a non-destructive and nucleus-selective analytical approach, enabling the development of highly specific sensing methods. Among the available nuclei, 1H detection is the most commonly employed due to its high sensitivity, although the use of other nuclei has steadily increased in low-field NMR applications as instrumentation has continued to improve .
A wide variety of NMR methods can be implemented by applying different radio-frequency (rf) pulse sequences, resulting in distinct signal responses that encode diverse physicochemical information . The analysis of these signals is typically performed using mathematical approaches available in spectrometer software, including the Fourier transform (FT); exponential fitting; the inverse Laplace transform (ILT); and, more recently, machine learning techniques . Despite these developments, advanced signal-processing algorithms – particularly ILT – are often restricted to proprietary software or require specialized expertise, motivating many users to seek alternative implementations in environments such as C++, Python, MATLAB, or OriginLab .
The objective of this paper is to present an open-access WebApp (https://nmr-ilt.esalq.usp.br/, last access: 15 April 2026) for the application of the inverse Laplace transform (ILT) to time-domain NMR (TD-NMR) signals, together with a concise methodological discussion of the underlying principles. The implemented algorithm is based on a non-negative least squares approach combined with Tikhonov regularization and singular value decomposition (SVD), enabling robust inversion of ill-posed relaxation data. In contrast to our previous ILT implementation reported in the literature , which focused on local computational routines and specific applications, the present work provides a platform-independent, freely accessible web-based framework suitable for TD-NMR data from instruments of any manufacturer. The platform supports the main TD-NMR experiments used in practice, including Carr–Purcell–Meiboom–Gill (CPMG), inversion recovery, and saturation recovery pulse sequences. The main features of the proposed approach are validated using simulated datasets and representative experimental signals, demonstrating its consistency with established ILT methodologies and its potential as an open scientific tool for research and teaching in magnetic resonance applications.
Magnetic resonance background and motivation for ILT
Magnetic resonance techniques can be broadly distinguished by magnetic field strength and homogeneity, which directly determine signal characteristics and the appropriate data analysis strategy. High-field NMR systems employ highly homogeneous magnetic fields, enabling frequency-domain analysis with high spectral resolution through the Fourier transform (FT). In contrast, time-domain and low-field NMR (TD-NMR/LF-NMR) systems are typically based on permanent magnets, resulting in compact, cost-effective, and robust instrumentation but often with limited magnetic field homogeneity.
In many TD-NMR and benchtop systems, magnetic field inhomogeneities prevent reliable frequency-domain resolution, making relaxometry-based approaches more suitable. Pulse sequences such as Carr–Purcell–Meiboom–Gill (CPMG), inversion recovery (IR), and related methods are therefore widely used to probe longitudinal (T1) and transverse (T2) relaxation processes. The signals acquired using these sequences consist of superpositions of exponentially decaying components, whose amplitudes and characteristic times reflect molecular dynamics and local interactions within the sample.
From a signal-processing perspective, relaxometry data represent an inverse problem in which the measured signal is modeled as a weighted sum of exponential decays. Unlike the Fourier transform, which provides a direct and well-posed transformation, the determination of relaxation time distributions from TD-NMR signals constitutes an ill-posed inverse problem. In such cases, small perturbations in the input data – such as noise or experimental uncertainty – can lead to large variations in the recovered solution, necessitating the use of constrained optimization and regularization strategies.
Within this context, the inverse Laplace transform (ILT) has emerged as a powerful framework for analyzing TD-NMR relaxometry data. Practical implementations of the ILT rely on numerical algorithms that balance stability, resolution, and robustness through appropriate regularization and non-negativity constraints. These methodological aspects are central to the present work and motivate the development of a transparent, accessible, and reproducible platform for ILT-based analysis of TD-NMR signals, as described in the following section.
Inverse Laplace transform
In time-domain and low-field NMR experiments, signal analysis is predominantly performed in the time domain, particularly through relaxometry-based methods. In heterogeneous systems, the measured signal generally reflects a superposition of contributions arising from distinct spin populations characterized by different relaxation times, leading to multi-exponential signal behavior.
The present work does not aim to review relaxometry techniques in detail, which are extensively discussed elsewhere , but instead focuses on the mathematical processing of multi-exponential signals acquired using pulse sequences such as Carr–Purcell–Meiboom–Gill (CPMG), inversion recovery (IR), saturation recovery (SR), and related methods.
In simple systems, relaxation data can often be adequately described by mono-exponential models. However, in more complex samples, the measured signal is better represented as a superposition of multiple exponential components, each associated with a distinct relaxation time. In such cases, conventional multi-exponential fitting becomes unstable and strongly dependent on the assumed number of components.
In these situations, the analysis of relaxation data is commonly formulated as an inverse problem and addressed using inverse Laplace transform (ILT) methodologies, which aim to recover continuous distributions of relaxation times from experimental signals .
Figure schematically illustrates the objective of ILT processing, in which a time-domain relaxometry signal is transformed into a relaxation time distribution. The positions and areas of the resulting peaks reflect the relative contributions of different spin populations present in the sample.
Comparative schematic illustration of time-domain NMR signal processing using the Fourier transform (FT) and the inverse Laplace transform (ILT). In the FT approach (left panel), which is the most widely known and routinely used method in NMR spectroscopy, the free induction decay (FID) consists of a mixture of oscillatory, exponentially damped sinusoids in the time domain. After Fourier transformation, these components are separated in the frequency domain, producing a spectrum whose peak positions, line widths, and amplitudes provide chemical and structural information about the sample. In contrast, the ILT approach (right panel) is applied to non-oscillatory signals. On the left side of this panel, a Carr–Purcell–Meiboom–Gill (CPMG)-type signal is represented by the echo amplitudes as a function of time, showing purely monotonic exponential decays rather than oscillations. By processing this signal with the ILT, the result is not a frequency spectrum but a distribution of relaxation times (T2), shown on the right, also called a relaxogram. The positions and areas of the T2 peaks reflect physical and compositional properties of the analyzed material, such as pore size distribution, water and oil content, and meat-to-fat ratio. Together, the figure highlights the conceptual analogy between both transforms – each mapping time-domain signals into characteristic distributions – while emphasizing the fundamental difference between oscillatory (FT) and non-oscillatory (ILT) data.
In this way, the ILT process reveals the contribution of different spin groups, where the two wide peaks illustrated in the T2 distribution can be a better representation of a real sample once real samples are not straightly homogeneous with singular values of amplitude and relaxation times (narrow peaks), but they are actually inhomogeneous samples, with a dispersion of relaxation times.
The distribution obtained in the processing with the ILT can arise from a growing or decaying exponential damped curve, as with CPMG (T2) or inversion recovery (T1), saturation recovery (T1), diffusion (D), among other experiments, where what changes in the processing is the kernel equation used in the transformation, as will be discussed in the next section.
There are several software platforms to perform this processing, some with free access and others only being supplied with licenses by the spectrometer manufacturing companies, each with different processing characteristics. The most famous ILT algorithms are CONTIN (CONTINuous distribution) , UpenWin (Uniform PENalty) , BRD (Butler–Reeds–Dawson) , TSVD (truncated singular value decomposition) , NNLS (non-negative least squares) , PDCO (primal–dual interior method) , among bi-dimensional ILT and other methods.
In the next section, the basic concepts of the ILT and a simple-to-use WebApp developed for performing the transformation of TD-NMR signals are presented, freely available for research and teaching. A mathematical description of the implemented processing is then provided, followed by results with simulated and experimental data and, finally, the procedures for using the program. A video tutorial demonstrating the code and its usage is available on the WebApp page or can be requested from the authors by email.
Mathematical description of the ILT
There are several mathematical approaches to perform the inversion of these TD-NMR signals, and these techniques are generally referred to as inverse Laplace transform by the NMR scientific community. It is worth noting that the method discussed here is different from the procedure of the same name found in mathematical physics books . Despite the confusion with the name, the origin of the uses of this term by the NMR community is related to an analogy with the procedure of the inverse Fourier transform. Further details of this history can be found in the work of Fordham and colleagues .
In the context of TD-NMR, for most of the pulse sequences used, the acquired signals show exponential and/or Gaussian growth or decay. The inversion problem to be solved here is to determine all of the amplitudes and relaxation times present in the data. Therefore, what we call ILT is a particular case of a more general problem of inverting Fredholm integral equations, which can be written in the following form:
f(x)=∫abϕ(t)K(x,t)dt,
where K(x,t) is the kernel function of the system, and the function ϕ(t) represents the function to be determined from the experimental data curve f(x).
In the case of an experimental signal from a CPMG (Carr–Purcell–Meiboom–Gill) pulse sequence, the model function F(t) is given by F(t)=c(0)exp(-t/T2), where c(0) is the initial amplitude of the experimental signal, t is the time variable of the acquired signal, and T2 is the relaxation time. For a sample composed of multiple relaxation times, we have the following:
F(t)=∑i=1kci(0)exp(-t/T2i),
where k denotes the number of components present in the sample.
In real complex samples, it is physically expected not to have only a singular value of relaxation time but rather to have a continuous distribution of relaxation times. Thus, we can express the NMR signal as given by the following integral:
F(t)=∫g(T2)exp(-t/T2)dT2,
where exp(-t/T2) is the kernel function, and g(T2) represents the distribution of relaxation times T2, providing the amplitudes of each infinitesimal component dT2.
The ILT algorithm will determine the distribution of amplitudes g(T2) from the experimental data c(tn); however, experimental NMR signals c(tn) are a set of n data points acquired over the discrete time points tn, with the presence of experimental noise ϵ(tn); thereby
c(tn)=∑i=1Ng(T2i)K(tn,T2i)+ϵ(tn),
where K is the kernel equation of this summation. Note that N denotes the total points of the T2 dimension, which need to be defined in a window of expected relaxation times, for example, from 0.001 to 10 s using N=100 or 200 points. For a CPMG experiment, the kernel should be defined as a composition of decreasing exponential, as shown in Eq. (). The inversion of this problem can be obtained through the minimization of the mean square errors, given by
X2=‖c(tn)-F(tn)‖2,
where F(tn) is the model function that best describes the experimental data c(tn). This minimization problem is not trivial as it involves a mathematically ill-posed problem .
Several regularization methods can be used to obtain solutions for this Fredholm integral problem, each presenting different characteristics of stability and reliability. Physically, in our TD-NMR experiments, it is reasonable to expect that our solution g(T2) is a continuous distribution of relaxation times without the existence of delta peaks, discontinuous points, or negative values of amplitude in the distribution (Fig. ). These constraints can be introduced into the program's algorithm through a term that smooths the distribution and enforces non-negativity, as proposed by the NNLS (non-negative least squares) method. This regularization of the solution enforces more stable and continuous distributions, and the regularization parameter (α) is responsible for the smoothing of this solution.
Different values of the regularization parameter α can yield optimal curve fits, resulting in distributions with certain variations. Determining which distribution best represents the physical system is not trivial and, therefore, involves some subjectivity as the fit strongly depends on the signal-to-noise ratio of the data. Various methods have been used and proposed to optimize the choice of the regularization parameter α, such as the L-curve method . In the present implementation, no automatic optimization strategy for the regularization parameter α is employed. Instead, α is selected by the user based on prior knowledge of the system and the signal-to-noise characteristics of the data, an approach commonly adopted in practical ILT analyses. Future versions of the WebApp will incorporate automated strategies for α selection.
Therefore, the algorithm of the WebApp uses the Tikhonov regularization method with singular value decomposition while assuming the non-negativity of the response, also known as NNLS (non-negative least squares), widely used by the NMR community.
ILT-NMR WebApp
An open-access WebApp was developed to enable inverse Laplace transform (ILT) analysis of TD-NMR data without the need for local software installation. The application is deployed using Docker and hosted on secure servers of the University of São Paulo (USP) within the InterNuvem infrastructure, ensuring stable and reliable access for the scientific community. The computational core of the WebApp was implemented in Python , making use of established numerical and visualization libraries, including NumPy and Plotly. The implemented algorithm represents a platform-independent adaptation of our previously reported ILT routine developed for OriginLab software .
The WebApp supports ILT processing of relaxation data acquired using Carr–Purcell–Meiboom–Gill (CPMG), inversion recovery (IR), and saturation recovery (SR) pulse sequences. For each experiment type, the corresponding kernel function K(t,T) is defined as follows:
CPMG: K(t,T)=[exp(-t/T)],IR: K(t,T)=[1-2exp(-t/T)],SR: K(t,T)=[1-exp(-t/T)],
where t denotes the experimental time axis, and T represents the relaxation time domain of the resulting ILT distribution.
Figure summarizes the data-processing workflow implemented in the WebApp. Input signals are provided as text-based files containing the experimental time axis and one or more corresponding signal amplitudes. After specifying the experiment type and processing parameters, the ILT inversion is performed, and the resulting relaxation time distributions are returned in graphical and numerical formats, allowing direct inspection and further analysis. The uploaded data are temporarily processed on the server and are automatically discarded after the analysis, without being stored or logged.
WebApp pipeline.
The WebApp interface was designed to provide flexible control over key processing parameters, including the definition of the relaxation time window, discretization of the T axis, baseline correction, signal normalization, and the choice of the regularization parameter α. While the present implementation requires manual selection of α, this approach allows the user to incorporate prior knowledge of the sample and signal-to-noise characteristics into the analysis. Automated strategies for α optimization will be incorporated into future versions of the platform. To facilitate reproducibility and broader use, the WebApp accepts multiple signals sharing a common time axis, enabling batch processing and comparative analysis of TD-NMR experiments.
Figure shows the Processing page of the user interface, available at: https://nmr-ilt.esalq.usp.br/ (last access: 15 April 2026). The full navigation bar on the left includes Home (main page with general information about the ILT process), Processing (the core section for data processing), Citation, and Contact (to send a message to the authors).
WebApp simplified processing interface.
In the Processing bar, it is required that information be provided regarding the input parameters and information:
Kernel type. Information regarding what type of data the user is processing must be provided. There are three options, namely CPMG, IR, and SR kernels, as shown in Eq. ().
Remove start points. If desired, some initial points of the data can be removed. It is useful when delaying with spectrometer dead-time effects or spurious data points.
Alpha. The regularization parameter α is always necessary to be defined, and it is related to the resolution of the spectrum. Usually, the larger the α value, the lower spectral resolution; however, the forcing resolution can provide unreliable spectra in data with a low signal-to-noise ratio. In this WebApp, its value is defined by the user. Typical values used are 0.01, 0.1, 1, and 10.
Starting T value. Set the minimum value of the relaxation times expected, in the timescale defined.
Ending T value. Set the maximum value of the relaxation times expected, in the timescale defined.
Number of points. Define the total number of points in the resulting distribution spectrum. Typical values are 100 or 200.
Timescale. It must be defined whether the time axis is in seconds, milliseconds, or microseconds. After processing, the relaxation time T axis in the ILT spectra will be of the same order of magnitude.
Force to zero/center. Spurious peaks can arise from a small displacement of the baseline of the input signals (offset), and this function can be used to perform this pre-processing step.
Normalization. This function will divide the data in each column by the first value of the signal, thus normalizing the signals.
Perform ILT. Once the processing is initiated, in a few seconds (or minutes if many columns are provided), the results are presented in the figures. An option to download all of the images and data points is also available.
The signal input needed to have a minimum of two columns: first, the time axis and, second, the amplitude signal c(t). The results of ILT processing are presented in panels (b) to (e), with (b) the relaxation time distribution, (d) the comparison of the fitting and input signal, and (e) its residues. All data results can be downloaded in an xlsx file.
Figure presents an example of the input data format and the corresponding ILT processing results, including the relaxation time distribution, fitting, residuals, and exported data. The input data must be provided in a .txt or .dat file, where the first column represents the time axis, and the second contains the signal intensity. Multiple signals can also be processed simultaneously by placing them in the subsequent columns, provided they share the same time axis (first column). In this case, all signals will be processed with the same input parameters, which is convenient for comparison and analysis of dynamic process experiments. The first row may optionally be used to define the axis labels.
After the ILT processing, resulting figures are generated in the WebPage. Figure b shows the relaxation time distribution obtained, with the x axis in logarithmic scale, (d) shows the comparison between the fitting and the input signal, and (e) its respective residues. All data and figures can be downloaded in a zip file.
Three example files with experimental data are available on the Home page. The ILT WebApp can be accessed at https://nmr-ilt.esalq.usp.br (last access: 15 April 2026), and a video tutorial demonstrating its use is available at https://www.youtube.com/watch?v=_n_sN_G3Cnk (last access: 15 April 2026).
Results and discussions
To evaluate the performance and robustness of the implemented ILT algorithm, a series of tests were conducted using simulated and experimental TD-NMR signals acquired under different conditions for CPMG, inversion recovery (IR), and saturation recovery (SR) experiments. The results presented below focus on the influence of regularization parameters, the fidelity of the recovered relaxation time distributions, and the consistency of the results with known ground-truth distributions and literature data.
Influence of the regularization parameter <inline-formula><mml:math id="M61" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula>
Figure illustrates the effect of the regularization parameter α on the relaxation time distribution obtained from the ILT of a bi-exponential CPMG signal with relaxation times of 100 and 500 ms. As expected, increasing α leads to broader distributions with reduced peak amplitudes while preserving the total spectral area. This behavior reflects the trade-off between resolution and stability that is characteristic of regularized inverse problems.
Effect of the regularization parameter (α) on the resulting relaxation distribution spectrum. The larger the value of α, the broader the peaks become.
The choice of α is particularly critical for signals with a limited signal-to-noise ratio, where excessive regularization may oversmooth physically relevant features, whereas insufficient regularization can amplify noise-induced artifacts. In practice, the selection of α is guided by balancing these two effects, aiming to minimize spurious oscillations while preserving meaningful spectral features.
A common strategy is to select the smallest α value that yields a stable and physically interpretable distribution, without the appearance of noise-driven peaks. In this context, visual inspection of the resulting distribution, combined with knowledge of the expected number of components and their approximate relaxation times, can provide a practical and effective criterion.
Although automated optimization strategies such as the L-curve method have been proposed in the literature , the present implementation allows manual selection of α, enabling the incorporation of prior knowledge about the sample and experimental conditions. This flexibility is particularly useful in heterogeneous systems, where the definition of an optimal α may depend on the specific analysis objective, such as resolution of closely spaced components or suppression of noise artifacts.
Validation using simulated and experimental signals
The algorithm was first validated using simulated TD-NMR signals generated from predefined relaxation time distributions. Broad log-Gaussian distributions were used to emulate realistic heterogeneous systems, and synthetic CPMG and IR signals were obtained by numerical integration of these distributions, followed by the addition of Gaussian white noise at different rms levels. Figure a–d present representative simulated signals with increasing noise levels, corresponding to rms values of 0 %, 1 %, 2 %, and 5 %, respectively. For each case, the corresponding ILT distribution obtained from the simulated signal is shown below, alongside the original distribution used to generate the signal, allowing a direct comparison between the recovered and true distributions.
Simulated CPMG signals were generated to evaluate the performance of the ILT WebApp. Panels (a)–(d) show signals with increasing rms noise levels, obtained by numerical integration of a predefined relaxation time distribution (shown below as dist, red line). The original distribution consists of two peaks centered at relaxation times of 100 and 500 ms. After processing the time-domain data using the ILT algorithm, the recovered distributions (ILT, black line) are presented below each signal. Several similar simulated datasets were analyzed, and the results were found to be consistent with those obtained using established ILT software.
From 0 % to 2 % rms, the recovered distributions show good agreement with the ground-truth distributions used in the simulations, accurately reproducing peak positions, widths, and relative amplitudes. The relative deviation of peak centers remained within a few percent for all tested cases, consistently with previously reported results for established ILT implementations. As the noise level increases from 2 % to 5 % rms, the recovered ILT distributions begin to exhibit spurious peaks, reflecting the inherent ill-posed nature of the ILT problem in the presence of noise, as widely documented.
Comparison of CPMG signals (left) and the corresponding relaxation time distributions (right) of maize seeds under different conditions: (a) dry seeds with moisture content below 10 %; (b) seeds after 24 h of water absorption; and (c) seeds after several days, when germination had begun, as illustrated in the inset above. These results are consistent with Song et al. (2022)'s analysis of maize seeds .
Beyond simulations, the algorithm was applied to experimental TD-NMR data acquired from maize seed samples at different hydration and germination stages (Fig. ). The 1H TD-NMR measurements were performed on an 11.3 MHz spectrometer (0.27 T for 1H resonance frequency), SLK-200 (SpinLock, Argentina), using a 30 mm probe at 30 °C. The Carr–Purcell–Meiboom–Gill (CPMG) pulse sequence was used, with 90 and 180° pulses of 9.0 and 18.5 µs, respectively, and an echo time of 200 µs, with a total of 4000 echoes, a recycle delay of 3 s, and 32 scans. For ILT processing, α=1 and 100 points were used.
The three ILT spectra of Fig. are from a set of nine maize seeds measured simultaneously: first, in (a), measured when they are dry with less than 10 % moisture; in (b), measured after the set of maize seeds absorbed water by 24 h; and, in (c), after some days, when they started germination, as showed in the insert illustration of a maize seed. It can be noted in (a) that the main peak, at first at T2=1 ms, shifted to 3 ms in (b) and also got wider. Also, the two smaller peaks in (a), around 40 and 150 ms, shift and widen in (b). These changes reflect the process of water absorption and the increase in heterogeneity inside the maize seeds. These processes can be monitored in more detail and studied to interpret what each peak represents inside the maize seed. These peaks of T2 relaxation times in seeds typically represent signals from different portions of water and/or oil mobility due the porous structures.
In Fig. c, the ILT spectrum of germinated maize seeds is presented, as shown in the inset illustration of a maize seed. These resulting relaxation time distributions reveal systematic shifts and broadening of the T2 peaks as water uptake and metabolic activity increase, reflecting changes in molecular mobility and structural heterogeneity within the seeds. These trends are consistent with previously reported TD-NMR studies of maize seeds and validate the applicability of the proposed approach to real agro-biological systems .
Taken together, these results demonstrate that the implemented ILT algorithm reliably recovers relaxation time distributions from both simulated and experimental TD-NMR data, producing results that are consistent with established methods. Moreover, the same core algorithm has been successfully applied by our research group in a wide range of previous studies, including food analysis , meat aging and color evaluation , physiological disorders in poultry and fruits , seed and plant science , enzymatic activity in cassava roots , and the detection of counterfeit spirits , further supporting its robustness and practical applicability.
Tips and tricks
To ensure accurate analysis of the ILT processing, certain practical aspects need to be emphasized, as illustrated in Fig. :
Signals with exponential decay and/or growth should be fully acquired during the experimental stage, with a sufficiently long baseline. Discontinuities in the initial points or a step in the final baseline, as in illustration (i), can lead to artificial (spurious) peaks in the resulting ILT spectrum.
Artificial peaks may arise in the ILT spectrum from slight displacements of the baseline (amplitude offset). To address this, the WebApp provides the option of a pre-processing function, namely “Force to zero/center”, for this correction.
The ILT spectrum is strongly influenced by the signal-to-noise ratio of the acquired data. Noisy experimental signals can result in unreliable distributions.
Note that the ILT spectrum will have the same time unit as the input signal. We recommend using times in seconds or milliseconds for the input signal as relaxation times of biologic samples typically fall within the range of 0.001 and 10 s.
Note that the ILT spectrum window (starting T value and ending T value) and the number of points (100) must be selected carefully. If the window is too short or too long, it can generate distributions with artificial and distorted peaks.
For studies comparing several signals, it is advisable to consistently use the same parameter values in both the experimental acquisition stage and the ILT processing. For the ILT processing, these parameters include the number of points, the α parameter, and the values of the starting T value and ending T value.
Illustration with tips and tricks that must be observed when processing signals with the ILT WebApp. Panel (i) shows the problem of discontinuity experimental data signal, panel (ii) shows offset amplitude, panel (iii) shows the signal-to-noise ratio, panel (iv) shows the time axis, panel (v) shows the selection of an expected window of relaxation times, and panel (vi) shows the uses of the same parameters values to compare several signals. Figure adapted from reference .
By paying attention to these practical considerations, more accurate and meaningful ILT analyses can be achieved.
Conclusion
This work presented the fundamental principles of the inverse Laplace transform (ILT) applied to time-domain NMR (TD-NMR) data and introduced an open-access WebApp designed to perform this analysis in a transparent and platform-independent manner. The implemented approach is based on non-negative least squares combined with Tikhonov regularization, providing a robust framework for the inversion of ill-posed relaxometry problems commonly encountered in TD-NMR experiments.
Validation using both simulated datasets with known ground-truth distributions and representative experimental signals demonstrated that the proposed implementation reliably recovers relaxation time distributions, yielding results consistent with those obtained using established ILT software. These findings confirm the methodological soundness of the algorithm and support its use as a complementary tool for TD-NMR data analysis in research and educational contexts.
The WebApp is freely accessible online at https://nmr-ilt.esalq.usp.br/ (last access: 15 April 2026), offering an accessible environment for ILT analysis without the need for proprietary software. Future developments will focus on extending the available functionalities and incorporating automated strategies for regularization parameter selection.
Code and data availability
Code is not publicly available because it is implemented within a web-based platform that is freely accessible to users. The source code is not distributed in order to ensure appropriate use within the intended platform and to avoid unauthorized commercial reuse. Data are available from the authors upon request.
Sample availability
Samples are available from the authors upon request.
Video supplement
A video tutorial demonstrating the use of the WebApp is available at https://youtu.be/_n_sN_G3Cnk (Moraes, 2026).
Author contributions
T. B. Moraes: conceptualization, methodology, validation, data curation, writing (review and editing), funding acquisition, supervision, project administration. G. V. Von Atzingen: software, data processing, formal analysis, writing (original draft). L. P. Mazzero: software, web platform development, implementation. W. S. Mendes: software, web platform development, implementation. M. B. Zacharias: investigation, experimental data acquisition, validation, writing (original draft). M. C. B. Cardinali: software, data processing, methodology, writing (review and editing).
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. The authors bear the ultimate responsibility for providing appropriate place names. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Acknowledgements
The authors would like to thank the Universidade de São Paulo (USP) and Instituto Federal de Educação, Ciência e Tecnologia de São Paulo (IFSP), campus de Piracicaba/SP for supporting the development of this project. This work was financially supported by the following grants from Brazilian agencies: Conselho Nacional de Desenvolvimento Científico e Tecnológico CNPq grant (no. 306611/2022-8) and Fundação de Amparo à Pesquisa do Estado de São Paulo FAPESP (grant no. 2025/09284-9).
Financial support
This research has been supported by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (grant no. 306611/2022-8) and the Fundação de Amparo à Pesquisa do Estado de São Paulo (grant no. 2025/09284-9).
Review statement
This paper was edited by Geoffrey Bodenhausen and reviewed by three anonymous referees.
References Arfken, G. B. and Weber, H. J.: Mathematical Methods for Physicists, 4th edn., Academic Press, San Diego, ISBN: 978-0120598151, 1995.Azachi, M. H. and Wiesman, Z.: Correlating Molecular and Textural Properties of Raw Soy-Based and Beef Burgers Using TD-NMR and TPA, J. Food Sci., 91, e70839, 10.1111/1750-3841.70839, 2026. Berman, P., Levi, O., Parmet, Y., Saunders, M., and Wiesman, Z.: Laplace Inversion of Low-Resolution NMR Relaxometry Data Using Sparse Representation Methods, Concept. Magn. Reson. A, 42, 72–88, 2013.Bizzani, M., Flores, D. W. M., Moraes, T. B., Colnago, L. A., Ferreira, M. D., and Spoto, M. H. F.: Non-invasive detection of internal flesh breakdown in intact Palmer mangoes using time-domain nuclear magnetic resonance relaxometry, Microchem. J., 158, 105208, 10.1016/j.microc.2020.105208, 2020.Blümich, B.: Introduction to compact NMR: A review of methods, TrAC-Trend. Anal. Chem., 83, 2–11, 10.1016/j.trac.2015.12.012, 2016. Blümich, B., Haber-Pohlmeier, S., and Zia, W.: Compact NMR, De Gruyter, ISBN: 9783110266283, 2014.Borgia, G., Brown, R., and Fantazzini, P.: Uniform-Penalty Inversion of Multiexponential Decay Data, J. Magn. Reson., 132, 65–77, 10.1006/jmre.1998.1387, 1998.Bortolotti, V., Brizi, L., Fantazzini, P., Landi, G., and Zama, F.: Upen2DTool: A Uniform PENalty Matlab tool for inversion of 2D NMR relaxation data, SoftwareX, 10, 100302, 10.1016/j.softx.2019.100302, 2019.Bruker Corporation: https://www.bruker.com (last access: 15 April 2026), 2023. Butler, J. P., Reeds, J. A., and Dawson, S. V.: Estimating Solutions of First Kind Integral Equations with Nonnegative Constraints and Optimal Smoothing, SIAM J. Numer. Anal., 18, 381–397, 1981.Caputo, P., Loise, V., Ashimova, S., Teltayev, B., Vaiana, R., and Oliviero Rossi, C.: Inverse Laplace Transform (ILT) NMR: A powerful tool to differentiate a real rejuvenator and a softener of aged bitumen, Colloid. Surface. A, 574, 154–161, 10.1016/j.colsurfa.2019.04.080, 2019.Cheng, S., Wang, X., Li, R., Yang, H., Wang, H., Wang, H., and Tan, M.: Influence of multiple freeze-thaw cycles on quality characteristics of beef semimembranous muscle: With emphasis on water status and distribution by LF-NMR and MRI, Meat Sci., 147, 44–52, 10.1016/j.meatsci.2018.08.020, 2019.Cônsolo, N. R., Samuelsson, L. M., Barbosa, L. C., Monaretto, T., Moraes, T. B., Buarque, V. L., Higuera-Padilla, A. R., Colnago, L. A., Silva, S. L., Reis, M. M., Fonseca, A. C., da S. Araújo, C. S., de S. Leite, B. G., Roque, F. A., and Araújo, L. F.: Characterization of chicken muscle disorders through metabolomics, pathway analysis, and water relaxometry: a pilot study, Poultry Sci., 99, 6247–6257, 10.1016/j.psj.2020.06.066, 2020.Cônsolo, N. R., Silva, J., Buarque, V. L., Samuelsson, L. M., Miller, P., Maclean, P. H., Moraes, T. B., Barbosa, L. C., Higuera-Padilla, A., Colnago, L. A., Saran Netto, A., Gerrard, D. E., and Silva, S. L.: Using TD-NMR relaxometry and 1D 1H NMR spectroscopy to evaluate aging of Nellore beef, Meat Sci., 181, 108606, 10.1016/j.meatsci.2021.108606, 2021.da Silva Ferreira, J., Moraes, T. B., Colnago, L. A., and Pereira, F. M. V.: Enzymatic Activity Prediction Using Time-Domain Nuclear Magnetic Resonance (TD-NMR) and Multivariate Analysis: A Case Study Using Cassava Roots, Appl. Magn. Reson., 49, 653–664, 10.1007/s00723-018-0995-0, 2018.Day, I. J.: On the inversion of diffusion NMR data: Tikhonov regularization and optimal choice of the regularization parameter, J. Magn. Reson., 211, 178–185, 10.1016/j.jmr.2011.05.014, 2011.de Lima, J. T., Zawadzki, S. F., Soares, F. L. F., Salome, K. S., Barison, A., de Moraes, T. B., and D'Oca, C. D. R. M.: The forensic ability of TD-NMR in detecting counterfeit spirits by analyzing bottle caps, Microchem. J., 191, 108896, 10.1016/j.microc.2023.108896, 2023.Engelsen, S. B. and van den Berg, F. W. J.: Quantitative Analysis of Time Domain NMR Relaxation Data, Springer International Publishing, Cham, 1–19, 10.1007/978-3-319-28275-6_21-1, 2017.FIT - Fine Instrument Technology: https://fitinstrument.com/ (last access: 15 April 2026), 2023. Fordham, E., Venkataramanan, L., Mitchell, J., and Valori, A.: What are, and what are not, Inverse Laplace Transforms, The Open-Access Journal for the Basic Principles of Diffusion Theory, Experiment and Application, 29, 1–8, 2017.Hansen, P. C.: Analysis of Discrete Ill-Posed Problems by Means of the L-Curve, SIAM Rev., 34, 561–580, 10.1137/1034115, 1992.Harris, C. R., Millman, K. J., van der Walt, S. J., Gommers, R., Virtanen, P., Cournapeau, D., Wieser, E., Taylor, J., Berg, S., Smith, N. J., Kern, R., Picus, M., Hoyer, S., van Kerkwijk, M. H., Brett, M., Haldane, A., Fernández del Río, J., Wiebe, M., Peterson, P., Gérard-Marchant, P., Sheppard, K., Reddy, T., Weckesser, W., Abbasi, H., Gohlke, C., and Oliphant, T. E.: Array programming with NumPy, Nature, 585, 357–362, 10.1038/s41586-020-2649-2, 2020.Istratov, A. A. and Vyvenko, O. F.: Exponential analysis in physical phenomena, Rev. Sci. Instrum., 70, 1233–1257, 10.1063/1.1149581, 1999. Jacobsen, N. E.: NMR Spectroscopy Explained, John Wiley & Sons, Ltd, New Jersey, ISBN: 978-0-471-73096-5, 2007.Lawson, C. L. and Hanson, R. J.: Solving least squares problems, in: Classics in Applied Mathematics, Chap. 26, 10.1137/1.9781611971217, 1976.Magritek: https://magritek.com (last access: 15 April 2026), 2023.Marassi, A. G., de Araújo-Ferreira, A. G., Oliveira, E. L., Vidoto, E. L. G., Amorim, A. D. F., Trevizan, W. A., and Bonagamba, T. J.: Transverse relaxation measurements for moving samples in the presence of strong magnetic field gradients, J. Magn. Reson., 354, 107522, 10.1016/j.jmr.2023.107522, 2023.Maus, A., Hertlein, C., and Saalwächter, K.: A Robust Proton NMR Method to Investigate Hard/Soft Ratios, Crystallinity, and Component Mobility in Polymers, Macromol. Chem. Phys., 207, 1150–1158, 10.1002/macp.200600169, 2006.Mitchell, J., Gladden, L., Chandrasekera, T., and Fordham, E.: Low-field permanent magnets for industrial process and quality control, Prog. Nucl. Mag. Res. Sp., 76, 1–60, 10.1016/j.pnmrs.2013.09.001, 2014.Monaretto, T., Moraes, T. B., and Colnago, L. A.: Recent 1D and 2D TD–NMR Pulse Sequences for Plant Science, Plants, 10, 10.3390/plants10050833, 2021.Moraes, T. B.: Transformada Inversa de Laplace para analise de sinais de Ressonancia Magnetica Nuclear de Baixo Campo, Quim. Nova, 44, 1020–1027, 10.21577/0100-4042.20170751, 2021.Moraes, T. B.: An open-access WebApp for Inverse Laplace Transform analysis of TD-NMR signals, YouTube, https://www.youtube.com/watch?v=_n_sN_G3Cnk (last access: 15 April 2026), 2026.Moraes, T. B. and Colnago, L. A.: Noninvasive analyses of food products using low-field time-domain NMR: a review of relaxometry methods, Braz. J. Phys., 10.1007/s13538-022-01055-1, 2022.Moreira, L. F. P. P., Ferrari, A. C., Moraes, T. B., Reis, R. A., Colnago, L. A., and Pereira, F. M. V.: Prediction of beef color using time-domain nuclear magnetic resonance (TD-NMR) relaxometry data and multivariate analyses, Magn. Reson. Chem., 54, 800–804, 10.1002/mrc.4456, 2016.Nanalysis: https://www.nanalysis.com/ (last access: 15 April 2026), 2023.Niumag: https://www.nmranalyzer.com/ (last access: 15 April 2026), 2023.Oxford: https://nmr.oxinst.com/ (last access: 15 April 2026), 2023.Provencher, S. W.: CONTIN: A general purpose constrained regularization program for inverting noisy linear algebraic and integral equations, Comput. Phys. Commun., 27, 229–242, 10.1016/0010-4655(82)90174-6, 1982.Song, P., Yue, X., Gu, Y., and Yang, T.: Assessment of maize seed vigor under saline-alkali and drought stress based on low field nuclear magnetic resonance, Biosyst. Eng., 220, 135–145, 10.1016/j.biosystemseng.2022.05.018, 2022.Song, Y.-Q.: Magnetic Resonance of Porous Media (MRPM): A perspective, J. Magn. Reson., 229, 12–24, 10.1016/j.jmr.2012.11.010, 2013.SpinLock: https://www.spinlock.com.ar/ (last access: 15 April 2026), 2023. Tikhonov, A. N.: On the solution of ill-posed problems and the method of regularization, Dokl. Akad. Nauk SSSR+, 151, 501–504, 1963.van Duynhoven, J., Voda, A., Witek, M., and As, H. V.: Time-Domain NMR Applied to Food Products, in: Annual Reports on NMR Spectroscopy, vol. 69, Academic Press, 145–197, 10.1016/S0066-4103(10)69003-5, 2010.Van Rossum, G. and Drake Jr., F. L.: Python reference manual, Centrum voor Wiskunde en Informatica Amsterdam, ISBN: 978-1441412690, 1995. Venkataramanan, L., Song, Y.-Q., and Hurlimann, M.: Solving Fredholm integrals of the first kind with tensor product structure in 2 and 2.5 dimensions, IEEE T. Signal Proces., 50, 1017–1026, 10.1109/78.995059, 2002.Williams, P.: 200 and More NMR Experiments: A Practical Course, 3rd Edition Edited by S. Berger and S. Braun (University of Leipzig). Wiley-VCH, Weinheim, Germany. 2004. xv + 838 pp. 17 × 24 cm. $89.95. ISBN 3-527-31067-3, J. Nat. Prod., 68, 632–632, 10.1021/np058231c, 2005.Zalesskiy, S. S., Danieli, E., Blümich, B., and Ananikov, V. P.: Miniaturization of NMR Systems: Desktop Spectrometers, Microcoil Spectroscopy, and “NMR on a Chip” for Chemistry, Biochemistry, and Industry, Chem. Rev., 114, 5641–5694, 10.1021/cr400063g, 2014.